Numerical Computation of Spectral Solutions for Sturm-Liouville Eigenvalue Problems

This paper focuses on the study of Sturm-Liouville eigenvalue problems. In the classical Chebyshev collocation method, the Sturm-Liouville problem is discretized to a generalized eigenvalue problem where the functions represent interpolants in suitably rescaled Chebyshev points. We are concerned with the computation of high-order eigenvalues of Sturm-Liouville problems using an effective method of discretization based on the Chebfun software algorithms with domain truncation. We solve some numerical Sturm-Liouville eigenvalue problems and demonstrate the efficiency of computations.