The close-coupling approximation method involves the numerical solution of systems of coupled second order ordinary differential equations. The solutions can display instability which is made apparent by dependence of the resonance energy on H (step-size). This instability has been examined and corrected. The comparative efficiency, time-wise and storage-wise, of a number of algorithms for the integration of the system of equations is presented.