A Neumann series of Bessel functions representation for solutions of Sturm-Liouville equations
This provides a computationally efficient representation for solving Sturm-Liouville problems, particularly beneficial for applications requiring solutions over wide spectral ranges.
The authors derive a Neumann series of Bessel functions representation for solutions of Sturm-Liouville equations, where the truncation error is independent of the spectral parameter ω for real ω, enabling efficient computation over large ω intervals. Error estimates and an algorithm for solving related problems are provided.
A Neumann series of Bessel functions (NSBF) representation for solutions of Sturm-Liouville equations and for their derivatives is obtained. The representation possesses an attractive feature for applications: for all real values of the spectral parameter $ω$ the difference between the exact solution and the approximate one (the truncated NSBF) depends on $N$ (the truncation parameter) and the coefficients of the equation and does not depend on $ω$. A similar result is valid when $ω\in\mathbb{C}$ belongs to a strip $|Imω|<C$. This feature makes the NSBF representation especially useful for applications requiring computation of solutions for large intervals of $ω$. Error and decay rate estimates are obtained. An algorithm for solving initial value, boundary value or spectral problems for the Sturm-Liouville equation is developed and illustrated on a test problem.