Spectral corrections for Sturm-Liouville problems (Q5946116)

The authors discuss the numerical integration of linear 1-D Sturm-Liouville problems in finite and infinite intervals. A shooting method is implemented both for the eigenvalues and the eigenfunctions. Newton iterations, with their typical advantages and risks, are used to calculate an approximation of an eigenvalue. A general formula for the global discretisation error is proposed together with an estimate for the eigenvalue error reported to work well for higher order eigenvalues. Solutions to a number of test examples support the findings.