A Neumann series of Bessel functions representation for solutions of SturmLiouville equations
Abstract
A Neumann series of Bessel functions (NSBF) representation for solutions of SturmLiouville 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 (the truncation parameter) and the coefficients of the equation and does not depend on . A similar result is valid when belongs to a strip . 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 SturmLiouville equation is developed and illustrated on a test problem.
1 Introduction
In the recent work [13] a new representation for solutions of the onedimensional Schrödinger equation
(1.1) 
was obtained in terms of socalled (see, e.g., [24], [25] and [2]) Neumann series of Bessel functions (NSBF). The representation possesses an attractive feature for applications: for all the difference between the exact solution and the approximate one (the truncated NSBF) depends on (the truncation parameter) and and does not depend on . A similar result is valid when belongs to a strip . This feature makes the NSBF representation especially useful for applications requiring computation of solutions of (1.1) for large intervals of . For example, as was shown in [13], the NSBF representation allows one to compute hundreds or if necessary even thousands of eigendata with a nondeteriorating for large and remarkable accuracy. In [15] the NSBF representation was extended onto perturbed Bessel equations.
In the present work we derive an NSBF representation for solutions of the SturmLiouville equation
(1.2) 
The coefficients are assumed to admit the application of the Liouville transformation (see, e.g., [6], [26]). The main result consists in an NSBF representation for solutions of (1.2) and for their derivatives, preserving the same attractive feature described above. Error and decay rate estimates are obtained. An algorithm for solving initial value, boundary value or spectral problems for (1.2) is developed and illustrated on a test problem.
Besides this Introduction the paper contains the following sections. Section 2 presents some well known facts concerning the Liouville transformation together with recent results from [11] showing how the system of formal powers associated with (1.2) is transformed by the Liouville transformation. In Section 3 we recall some relevant results from [13] which are used in Section 4 and Section 5 for obtaining the main results of this work, the NSBF representation for solutions of (1.2) as well as for their derivatives. In Section 6 convenient for computation formulas for the coefficients of the NSBF representations are derived. In Section 7 a computational algorithm based on the NSBF representation is formulated and discussed. Section 8 presents some illustrations of its numerical performance on a test problem admitting an exact solution. In Appendix A we prove error and decay rate estimates of the NSBF representation.
2 Preliminaries on the Liouville transformation
2.1 Definition of the Liouville transformation
Consider the SturmLiouville differential equation
(2.1) 
with being a finite interval. Let , and be such that , and for all . Define the mapping by
(2.2) 
where . Denote . Then (2.1) is related to the onedimensional Schrödinger differential equation
(2.3) 
by the Liouville transformation of the variables and into and defined as (see, e.g., [6], [26])
and the coefficient is given by the relation
(2.4) 
The Liouville transformation can be considered as an operator acting according to the rule
Let us introduce the following notations for the differential expressions
The following proposition summarizes the main properties of the operator .
Proposition 2.1.

The inverse operator is defined by .

The uniform norms of the operators and are and .

The operator equality
is valid on .
2.2 Transformation of formal powers
Let be a nonvanishing solution of the equation
(2.5) 
such that
On the existence of a nonvanishing see Remark 2.5 below. Note that in general complex valued solutions are considered even when the coefficient is real valued. Denote .
Definition 2.2.
The families of functions and constructed according to the rules
are called systems of formal powers associated with equation (2.3).
Remark 2.3.
Analogously, let us introduce a system of formal powers corresponding to equation (2.1).
Let be a solution of the equation
(2.6) 
such that for all (see Remark 2.5). Then the following two families of auxiliary functions are well defined
Definition 2.4.
Remark 2.5.
The existence of a nonvanishing solution of (2.6) for complex valued and such that , and are continuous on was proved in [14, Remark 5] (see also [4]). Moreover, the only reason for the requirement of the absence of zeros of the functions and is to make sure that the formal powers be well defined. As was shown in [16] this is true even when and/or have zeros, though in that case corresponding formulas are slightly more complicated.
Theorem 2.6 ([11]).
In order that the equality be fulfilled, must be chosen so that
(2.7) 
We will assume this initial condition to be satisfied.
It is easy to see that for one obtains
(2.8) 
3 Solution of the onedimensional Schrödinger equation
Let , and , denote the solutions of (2.3) satisfying the following initial conditions in the origin
(3.1) 
where is an arbitrary complex number.
Theorem 3.1 ([13]).
The solutions and of the equation
admit the following representations
(3.2) 
and
(3.3) 
where stands for the spherical Bessel function of order , the functions are defined as follows
(3.4) 
where is the coefficient of in the Legendre polynomial of order . The series in (3.2) and (3.3) converge uniformly with respect to on and converge uniformly with respect to on any compact subset of the complex plane of the variable . Moreover, for the functions
(3.5) 
and
(3.6) 
the following uniform estimates hold
(3.7) 
for any , , and
(3.8) 
for any , belonging to the strip , . Here is a function satisfying as . These estimates are slight refinements for those presented in [13] using the ideas from [15]. Some estimates for are presented in Appendix A.
Inequalities (3.7) and (3.8) reveal an interesting feature of the representations (3.2) and (3.3): the accuracy of approximation of the exact solutions by the partial sums (3.5) and (3.6) does not depend on meanwhile it belongs to a strip in the complex plane. This feature was tested and confirmed numerically in [13]. An analogous result is valid for the derivatives of the solutions. We formulate it in the following statement.
Theorem 3.2 ([13]).
The derivatives of the solutions and with respect to admit the following representations
(3.9) 
and
(3.10) 
where are defined as follows
(3.11) 
The series in (3.9) and (3.10) converge uniformly with respect to on and converge uniformly with respect to on any compact subset of the complex plane of the variable .
Moreover, for the approximations
and
the following inequalities are valid
for any , , and
for any , belonging to the strip , .
Remark 3.3 ([13]).
The functions also can be constructed as solutions of the recurrent equations
(3.12) 
with the first functions given by and and initial conditions where .
The functions can be calculated from the equalities
(3.13) 
One can use as well.
4 Solution of the SturmLiouville equation
Theorem 4.1.
Let the functions , and satisfy the conditions from Subsection 2.1 and be a solution of (2.6) satisfying (2.7) and such that for all . Then two linearly independent solutions and of equation (2.1) for can be written in the form
(4.1) 
and
(4.2) 
with defined by (2.2), the coefficients being defined by the equalities
(4.3) 
where are from Definition 2.4. The solutions and satisfy the following initial conditions
(4.4)  
(4.5) 
The series in (4.1) and (4.2) converge uniformly with respect to on and converge uniformly with respect to on any compact subset of the complex plane of the variable .
Moreover, for the functions
(4.6) 
and
(4.7) 
the following estimates hold
(4.8) 
for any , , and
(4.9) 
for any , belonging to the strip , , where is a function from Theorem 3.1.
Proof.
Consider solutions (3.2) and (3.3) of (2.3). The functions
and
are then solutions of (2.1) with . In these representations all the magnitudes except the coefficients can be defined with no reference to equation (2.3). Let us show that the same observation is applicable to the functions
Note that
(4.10) 
for any function and any . Then
Due to Theorem 2.6 we obtain (4.3) and thus the solutions and of (2.1) have the form (4.1), (4.2).
Remark 4.2.
Inequalities (4.8) and (4.9) show that the representations (4.1) and (4.2) preserve the important property of the representations (3.2) and (3.3): the accuracy of approximation of the exact solutions by the partial sums (4.6) and (4.7) does not depend on meanwhile it belongs to a strip in the complex plane.
Remark 4.3.
Similarly to the coefficients (see Remark 3.3) the functions can be constructed as solutions of the recurrent equations
(4.12) 
obtained by applying to (3.12) and taking into account (4.10) and the equality . The first functions of the sequence are given by
(4.13) 
The initial conditions satisfied by have the form
(4.14) 
This sequence of equations for leads to a stable recursive integration procedure which is proposed in Section 6.
5 Representation of the derivatives of the solutions and
Due to (4.11) for and with the aid of (3.9), (3.10) and (4.1), (4.2) we obtain
(5.1) 
and
(5.2) 
where
and the functions and have the form (cf. [11, Remark 4.7])
(5.3) 
and
(5.4) 
Let us obtain a formula for calculating which would not involve magnitudes related to equation (2.3) but only those related to (2.1). Consider (3.11). Direct calculation gives us the relation
Hence, due to Theorem 2.6,
Thus,
(5.5) 
This is a direct formula for calculating in terms of the functions and . A recurrent formula for can be obtained from (3.13). We have
(5.6)  
(5.7) 
and
(5.8) 
One may also use
(5.9) 
6 Computation formulas for the coefficients
Besides the direct formulas (