[en] A new algorithm is presented for solving nonlinear second-order coupled equations of the form y'' = f(r, y). This method consists of a predictor, a corrector and a modifier so that it does not require iteration or matrix inversion. The method retains both the advantages of exponentially fitted two-step methods and of the predictor-corrector Numerov method. (orig.)