A Fortran program which implements this algorithm is provided and some test results are. This approach can be used to compute scattering phase shifts in an efficient and reliable manner. This algorithm uses a high-order, variable step Runge-Kutta like method in the region where the potential term dominates, and an exponential or Bessel fitted method in the asymptotic region. An efficient algorithm is described for the accurate numerical integration of the one-dimensional Schrödinger equation.
