An Explicit Single-Step Nonlinear Numerical Method for First Order Initial Value Problems (IVPs) ()
1. Introduction
Many of the numerical methods for obtaining the solution of the first order ordinary differential equation
(1)
are based on the assumption that the solution is locally representable by a polynomial. However, when a given initial value problem or its theoretical solution
is known to posse a singularity, then it is particularly inappropriate to represent
, in the neighbourhood of the singularity by a polynomial [1] [2]. This is evident as Runge-Kutta type methods, Obrechkoff methods and general linear multistep methods usually produce very poor solutions around singularity points [1] [3] [4] [5]. The authors in [4] were the first to develop quadrature formulas based on rational interpolating functions. On the other hand, the rational interpolation schemes proposed in [2], was seen to be effective in the neighbourhood of the singularity and even beyond as reiterated in [5]. The work of authors in [4] was modified by [2] with a replacement of the general rational function Luke et al. replaced the general rational function of [4] by
where
and
are respectively polynomial of degree m and n. The resulting schemes require analytic generation of first and higher order derivatives, hence, the major limitation of the schemes. Since rational functions are more appropriate for the representation of functions close to singularities than polynomials, the limitation is overcome by a local representation of the theoretical solution with a rational expression. Interestingly, this approach appears to be promising as several methods are now being constructed in this direction [6] [7] [8] [9] [10]. The works of the authors in [6] [7] [9] [10] [11] [12] showed that solution around singularity point are well approximated by this approach. In this work, an explicit single-step nonlinear method involving higher derivatives of the state function for solving (1) is presented. The local truncation error and absolute stability of the new method are also discussed.
2. Construction of Method
In this work, we assumed that the theoretical solution
of (1) can locally be represented by a rational interpolant
, of the form
(2)
To construct an explicit single-step method with (2) for (1), it requires that
satisfies the following:
(3)
Substituting for expressions and simplifying (3) yields
(4)
Eliminating the undetermined coefficients
,
,
,
,
and
in (4) results in
(5)
The resulting method (5) is explicit, self-starting and nonlinear. We shall refer to (5) as NLM4 which is the method proposed in this work. The new method NLM4 is suitable for solving initial value problems whose solution possesses singularities.
3. Local Truncation Error and Absolute Stability of Constructed Method
In this section, the local truncation error (lte) and the absolute stability properties of the new method proposed in this work are considered.
3.1. Local Truncation Error
Local Truncation Error: The local truncation error
at
of the general explicit one step method
(6)
is given as
(7)
where,
is the theoretical solution. Using the above definition, it follows that the local truncation error of the constructed one step method can be written as
(8)
3.2. Order of a Ordinary Differential Equation
Order of a Ordinary Differential Equation: A numerical method is said to be of order p if p is the largest integer for which
for every n and
. Following the above definition, the local truncation error of the method constructed in this work is obtained as the residual when
is replaced by
. Below is the local truncation error for the method constructed in this work.
(9)
3.3. Consistency
A scheme is said to be consistent if the difference equation of the integrating formula exactly approximates the differential equation it intends to solve as the step size approaches zero. In order to establish the consistency property of the constructed method, it is sufficient to show that
(10)
Now,
(11)
the above indicates that the constructed schemes satisfy the consistency property.
3.4. Stability
To get the stability behaviour of the constructed scheme, the scheme is implemented on the standard test problem
(12)
and the stability polynomial
, where
is obtained. The stability function of (5) is obtained as
(13)
and the region of absolute stability is seen in Figure 1.
Figure 1. Region of absolute stability of (5).
4. Numerical Examples
The first problem considered in this work is the nonlinear initial value problem
(14)
whose theoretical solution is given as
(15)
For this problem, the absolute errors of the results obtained by the method proposed in this work are first compared with those of Non-linear One-Step methods for initial value problems of [7] and the derivative-free methods proposed in [11] as shown in Figure 2. A comparison of the maximum absolute error obtained by the proposed methods against those produced by the methods of the authors in [4] [6] [7] [9] [13] is also presented in Figure 3.
Figure 2. Logarithm of absolute errors for the solutions of Problem 1 with step-size h = 0.01.
Figure 3. Log plot of maximum absolute errors for Problem 1 as a function of the step-size
.
Figure 4. Logarithm of absolute errors for the solutions of (16) with step-size h = 0.01.
Problem 2
The second test problem considered is given as
(16)
The exact Solution is
(17)
The logarithm of absolute errors for the solutions obtained is compared with other methods discussed in [12] as given in Figure 4.
5. Conclusion
The explicit single-step nonlinear method constructed in this work is consistent and absolutely stable. Its region of absolute stability is larger than those of the methods discussed in the literature. The method gave more accurate result on the standard test problems compared with other methods discussed. Hence, the method is suitable for solving problems whose solution possesses singularity.