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 computations' efficiency.