A problem in computation of pseudoinverses
In this paper we describe implementation of three methods for computing generalized inverses in the programming packages MATHEMATICA and DELPHI. Grevile's method, the limit representation and the determinantal representation of generalized inverses are considered. A major problem arising in the implementation of these methods is very rapid increment of arithmetic operations. We describe a few different algorithms to avoid this common problem.