An iterative algorithm for estimating the Moore-Penrose generalized inverse is developed. The main motive for the construction of algorithm is simultaneous usage of Penrose equations (2) and (4). Convergence properties of the introduced method are considered as well as their first-order and the second-order error terms. Numerical experience is also presented. AMS Subj. Class.: 15A09.