numerical matrix methods.xmcd

This worksheet is based on an article in AJP using numerical matrix methods to estimate energies and r-component wave functions for central potentials.

The Mathcad results appear internally consistent, but do not agree at all with the results of the authors who used Mathematica.  Anyone see where I have made an error?