Use a linearly constrained linear least squares solver.
lsqlin in MATLAB https://www.mathworks.com/help/optim/ug/lsqlin.html
lsei in R https://rdrr.io/rforge/limSolve/man/lsei.html
The easiest way might be using CVX https://cvxr.com/cvx under MATLAB. CVXR https://cvxr.rbind.io/ under R, or CVXPY https://www.cvxpy.org/ under Python.
Here is the code for CVX:
A*x >= 0
which will transform the problem into a Second Order Cone Problem, send it to a solver, and transform the solver results back to the original problem as entered. You can include the factor of 1/2 (harmless) and square the norm, which doesn't affect the solution but needlessly makes the problem solution less numerically robust.
Edit: Extra details as requested in chat:
CVX calls a numerical optimization solver to solve the optimization problem. The solver enforces the specified constraints (within solver tolerance).
As mentioned above, CVX actually transforms this into an SOCP (Second Order Cone Problem) by converting the problem into epigraph formulation. It does this by introducing a new variable, t, and in effect moving the original objective to the constraints. Thus produce the problem.
norm((A*x-b) <= t
A*x >= 0
There might also be a slight rearrangement of the constraint
A*x >= 0. CVX calls a Second Order Cone solver optimization solver such as SeDuMi, SDPT3, Gurobi, or Mosek to solve this problem. It then transforms the results back to the original problem formulation as entered by the user.