Approximations to the Moore-Penrose generalized inverse are obtained via iteration in the application of regularization. Uniform error bounds are obtained for linear operators with closed range. For operators with arbitrary range pointwise error estimates are derived assuming certain smoothness conditions on the data. The stability of the iteration is considered and error bounds are obtained for“noisy”data.