ABSTRACT

We describe SLDRIVER, an interactive Fortran package for experimenting with the numerical solution of Sturm-Liouville problems. SLDRIVER incorporates into one program four of the main Sturm-Liouville solver codes, namely SL02F, SLEDGE, SLEIGN and SLEIGN2. A Fortran implementation of the 60 test problems in the author’s 1993 text “Numerical Solution of Sturm-Liouville problems” is included in SLDRIVER. This makes it easy to compare the performance of different solvers on various problems. Other test problem sets can be incorporated simply by including a different subroutine at link time.

A general boundary condition can be specified not only at a regular endpoint but also at a limit-circle (nonoscillatory or oscillatory) endpoint for each problem in the set, allowing the recently developed facilities of SLEIGN2 to be exercised. Eigenfunctions can be computed as well as eigenvalues.

The program records computed results in the form of MATLAB m-files. This makes it simple to plot eigenfunctions, or do other analyses, using the powerful facilities of MATLAB.