Richard Slevinsky

Philippe Gaudreau, Richard Slevinsky, Hassan Safouhi

Sturm-Liouville problems are abundant in the numerical treatment of scientific and engineering problems. In the present contribution, we present an efficient and highly accurate method for computing eigenvalues of singular Sturm-Liouville boundary value problems. The proposed method uses the double exponential formula coupled with Sinc collocation method. This method produces a symmetric positive-definite generalized eigenvalue system and has exponential convergence rate. Numerical examples are presented and comparisons with single exponential Sinc collocation method clearly illustrate the advantage of using the double exponential formula.

Philippe Gaudreau, Richard Slevinsky, Hassan Safouhi

A quantum anharmonic oscillator is defined by the Hamiltonian ${\cal H}= -\frac{ {\rm d^{2}}}{{\rm d}x^{2}} + V(x)$, where the potential is given by $V(x) = \sum_{i=1}^{m} c_{i} x^{2i}$ with $c_{m}>0$. Using the Sinc collocation method combined with the double exponential transformation, we develop a method to efficiently compute highly accurate approximations of energy eigenvalues for anharmonic oscillators. Convergence properties of the proposed method are presented. Using the principle of minimal sensitivity, we introduce an alternate expression for the mesh size for the Sinc collocation method which improves considerably the accuracy in computing eigenvalues for potentials with multiple wells. We apply our method to a number of potentials including potentials with multiple wells. The numerical results section clearly illustrates the high efficiency and accuracy of the proposed method. All our codes are written using the programming language Julia and are available upon request.