Received 23 October 2015; accepted 22 February 2016; published 25 February 2016
The theory of ordinary homogeneous linear differential equations of the second order, containing a large parameter, is well established  -  . The aim of this paper is to investigate detailed analytical solutions of equations of the form;
where is and is a real parameter. We shall investigate the behaviour of solutions of this differential equation, and the stability of the origin as. Without loss of generality, we take First, we make the following remarks:
a) Any second order linear ODE of the form; can be reduced to by a suitable transformation.
b) Furthermore, any equation of the form is conservative. We shall demonstrate this shortly. This will help us in our asymptotic stability analysis.
c) In Equation (1.1) if we take; then, we have the well known Sturm-Liouville problem;
where, is positive and of class and.
Introducing the new variables;
If we suppress the variable t for the moment, it then follows that;
Since, then, and after transforming the interval into, with further algebraic manipulations, the ode (1.2) becomes;
is a continuous function of. It can be shown that the solutions of (1.3) satisfy the Volterra integral equation;
where and are arbitrary. and have the same value, and the same derivate, at. The solution to the integral Equation (1.5) can be obtained by successive approximation in the form;
Assuming that the function is bounded, i.e., there exists a constant M such that, then, one can show by induction that;
In the case of a finite interval, it follows that (1.7) is uniformly convergent for, and is also an asymptotic expansion of as. Unfortunately, the is very difficult to compute. Other approximations for large may be obtained from formal solutions, and these are usually divergent.
2. Formal Solutions
Let us now consider the general ode;
If is a formal power series in with coefficients which depend on x, then two linearly independent solutions of (2.1) may also be represented by a formal power series in. However if the formal expansion of in powers of contains positive powers of, then the formal expansion of x will be a Laurent series. We shall discover that in the case that, as a function of, has a pole at, we can still construct formal solutions.
In (2.1), we shall assume that is of the form;
where the are independent of, and Furthermore, we assume that does not vanish in the interval over which t varies. We shall adopt a first formal solution of the form;
Substituting (2.3) into (2.1), with the convention that for and for and also for All summations may then be assumed over all the integers, and we obtain
Picking out the coefficients of we obtain;
This first condition arises when Setting in (2.4) we obtain;
It then follows that;
Consequent upon these relations, we may restrict our summation to in the first sum in Equation (2.4). Now for in (2.4) we get;
and when we replace n by in (2.4) we obtain
It is now obvious that Equation (2.3) satisfies (2.1), provided that and satisfy (2.5) to (2.8).
In these equations, empty sums (i.e. those with upper limit
we may choose a branch of, and then (2.5) determines up to an additive constant. Moreover,
, and hence (2.6) determines recurrently, up to an additive constant in each. Equation (2.7) determines up to a constant factor, and (2.8) determines recurrently, up to an additive con-
stant multiple of in each. Corresponding to the two branches of, we obtain two formal solutions of the form (2.3).
3. Another Formal Solution
A second type of formal solution is given by
Substituting (2.9) into (2.1) we get;
Equating coefficients of,
There are two linearly independent formal solutions of this type. The obvious connection between these two types of formal solutions can be seen from the fact that equations (2.10) and (2.11) are identical with (2.5) and
(2.6), and is the formal expansion of.
In the foregoing, we have assumed that as a function of, has a pole of even order at. If the pole is of odd order, then no solution of the form (2.3) or (2.9) exists, and instead of powers of, we must expand in powers of.
3.2. Asymptotic Solutions
We shall now demonstrate that under certain assumptions, the differential Equation (2.1) possesses a fundamental system of solutions which are represented asymptotically by the formal solutions obtained in preceding section. It actually does not matter whether we compare solution of (2.1) with
where the and satisfy (2.5) to (2.8), or with
where the satisfy (2.10) to (2.12), for the q’s and’s can be so chosen that the ratio of these two expressions is.
We now fix a positive integer N, and set;
with, and for each j, the satisfy (2.10) to (2.12). These coefficients are completely determined by, and certain derivatives of these functions, and we shall say that the are sufficiently often differentiable if all the derivatives entering the determination of the exist and are continuous functions of t. We allow t to vary over a bounded and closed interval, and over a sectorial domain. We have the following theorem.
Let S and I be as defined above then for each fixed is a continuous function of t over I; If
Uniformly in t and, as in S, where the are sufficiently often differentiable in I, and
, then the differential equation
possesses a fundamental system of solutions, and, such that
Top establish the existence and asymptotic property of, we substitute
in Equation (3.4) to get
uniformly in t and in S, by (3.2) and (2.10) to (2.12). Equation (3.7) may be written as
By two successive integrations, and a suitable choice of the constants of integration, we obtain;
Since is an increasing function, we have and .
The existence of follows immediately from the theory of Volterra integral equations, or can be established by successive approximations. From (3.8) and (3.9), we have, uniformly in t and in S. Furthermore is differentiable, i.e.
This proves (3.5) for. The proof for is very much similar, except that b rather than a, must be chosen as fixed limit of integration in the integral equation.
The methods of the last two sections can be applied to prove the asymptotic formulae for the Bessel functions  , viz;
Equation (1) holds for, uniformly in if
Equations (2) and (3) hold for uniformly in if
We observe that the functions; are solutions of the differential equation
This equation is of the form (3.4) with all other vanishing identically. The points are singular points of (3.10) and is a so called transition point at which the condition (3.3) is violated for any value.
5. Stability Analysis
In Section 1.0, we claimed that any equation of the form is conservative. It turns out that such systems are characterized by closed curves in the phase plane. For the former equation, we only need to show that it possesses a Hamiltonian H, such that.
Let us begin by multiplying the equation by, i.e.,
Hence (3.11) becomes
Thus, the required Hamiltonian is.
The Bessel differential equation
can be recast in vector form as
Clearly the origin (0, 0) is the only critical point and the corresponding Hamiltonian is;
We use the above Hamiltonian to construct a Lyapunov function given by;
with and. We note that, furthermore;
Thus. Since, it follows that and hence the origin is asymptotically stable for all and.
6. Numerical Investigation of Asymptotic Solutions
In what follows, we employ the Runge-Kutta algorithm provided by MathCAD  software to obtain a numerical solution for large values of.
6.1. MathCAD Runge-Kutta Algorithm
We define the following for the MathCAD algorithm.
t0: = 0.2 t1: = 10 Solution interval endpoint
Initial condition vector
N: = 1500 Number of solution values on [t0, t1]
S: = rkfixed (ic, t0, t1, N, D) Runge-Kutta algorithm.
T: = S<0> Independent variable values.
X0: = S<1> First solution function values.
X1: = S<2> Second solution function values.
Remark: X0 represents solution values x satisfying the Bessel ODE, while X1 represents the derivative of X0 i.e.. S
Figure 1. Section of solution matrix S.
For solutions no longer exist as they become unbounded. From the graphs shown, it is clear that the given Bessel differential equation is very sensitive to the parameter, and as the effect is to increase the oscillations until the solutions become unstable and die out. Furthermore, the phase portrait depicted shows that the Bessel differential equation represents a conservative system. This is clearly evident from the closed curves. However for, the phase portrait no longer appears like a closed curve but more like an explosion from the centre.
In this work, we have studied asymptotic solutions of equations of the type, where is a large parameter. We have shown that equations of this form represent a conservative system, meaning that they possess a conserved quantity, namely the Hamiltonian which is computed. As a special example, we consider the
Bessel differential equation for which the behaviour of the solutions as well as
the stability of the origin is investigated numerically as the parameter grows indefinitely.
 MathSoft, Inc. (1997) MathCAD 7 User’s Guide.