A method for computing the generalized inverse of a matrix is described, which makes use of elementary orthogonal matrices and theGaussian elimination. The method also yields orthonormal bases for the ranges and the null spaces of the matrix and the generalized inverse. Modifications of the method for the solution of simultaneous linear equations are given. Compact storage schemes, in the case of sparse matrices, are also described.
展开▼