Abstract
The modified Gram–Schmidt process of Rice [4] is adapted to give a numerically stable algorithm for computing the minimal norm least squares solution of the problem Ax = b. Numerical results of a FORTRAN implementation of this algorithm are also discussed