Computation of eigenvalues and solutions of regular Sturm–Liouville problems using Haar wavelets

The paper presents a novel method for the computation of eigenvalues and solutions of Sturm–Liouville eigenvalue problems (SLEPs) using truncated Haar wavelet series. This is an extension of the technique proposed by Hsiao to solve discretized version of variational problems via Haar wavelets. The proposed method aims to cover a wider class of problems, by applying it to historically important and a very useful class of boundary value problems, thereby enhancing its applicability. To demonstrate the effectiveness and efficiency of the method various celebrated Sturm–Liouville problems are analyzed for their eigenvalues and solutions. Also, eigensystems are investigated for their asymptotic and oscillatory behavior. The proposed scheme, unlike the conventional numerical schemes, such as Rayleigh quotient and Rayleigh–Ritz approximation, gives eigenpairs simultaneously and provides upper and lower estimates of the smallest eigenvalue, and it is found to have quadratic convergence with increase in resolution.