Module SLEIG in IMSLMSolves Sturm-Liouville problems of the form -(p*u'')''+q*u = lambda*r*u for eigenvalues, eigenfunctions and/or spectral density functions.
CALL SLEIG (CONS, COEFFN, ENDFIN, NUMEIG, INDEX, TEVLAB, TEVLRL, EVAL)| I1b3 | Eigenvalue (e.g., Sturm-Liouville) multipoint boundary value problems for ordinary differential equations |