Abstract
A method for computing the generalized inverse of a matrix is described, which makes use of the well-known Gauss–Jordan elimination method and the orthogonal triangularizations.