Spectral theory for bounded linear operators is used to develop a general class of approximation methods for the Moore-Penrose generalized inverse of a closed, densely defined linear operator. Issues of convergence and stability are addressed and the methods are modified to provide a stable class of methods for evaluation of unbounded linear operators.