Module 810 in TOMSSLEIGN2: Computation of eigenvalues and eigenfunctions of Sturm-Liouville problems, regular and singular self-adjoint, with both separated and coupled boundary conditions. Approximates the continuous spectrum in the singular case. Upgrades and extends the earlier SLEIGN software (TOMS Algorithm 700). (See P. Bailey et al., ACM TOMS 27 (2001) pp. 143-192).