Computing eigenelements of Sturm-Liouville problems by using Daubechies wavelets

Abstract

This work is our first step to get multiresolution approximation of eigenelements of Sturm-Liouville problems within bounded domain of varied nature. The formula for obtaining elements of representation of Sturm-Liouville operator involving polynomial coefficients in wavelet basis of Daubechies family have been derived in a form which can be readily used for their computations by a simple computer program. Estimates of errors for both the eigenvalues and eigenfunctions are also presented here. The proposed wavelet-Galerkin scheme based on scale functions and wavelets of Daubechies family having three or four vanishing moments of their wavelets has been applied to get approximate eigenelements of regular and singular Sturm-Liouville problems within bounded domain and compared with the exact or approximate results whenever available. From our study it appears that the proposed method is efficient and rapidly convergent in comparison to other approximation schemes based on variational method in Haar basis or finite difference methods studied by Bujurke et al. [39].