We describe a new algorithm to compute the eigenvalues of singular Sturm-Liouville problems with separated self-adjoint boundary conditions for both the limit-circle nonoscillatory and oscillatory cases. Also described is a numerical code implementing this algorithm and how it compares with SLEIGN. The latter is the only effective general purpose software available for the computation of the ei...
