An Explicit Single-Step Nonlinear Numerical Method for First Order Initial Value Problems (IVPs)

Interest in the construction of efficient methods for solving initial value problems that have some peculiar properties with it or its solution is recently gaining wide popularity. Based on the assumption that the solution is representable by nonlinear trigonometric expressions, this work presents an explicit single-step nonlinear method for solving first order initial value problems whose solution possesses singularity. The stability and convergence properties of the constructed scheme are also presented. Implementation of the new method on some standard test problems compared with those discussed in the literature proved its accuracy and efficiency.


Introduction
Many of the numerical methods for obtaining the solution of the first order ordinary differential equation , , y f x y x x x X y x η ′ = ∈ = (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 ( ) u t is known to posse a singularity, then it is particularly inappropriate to represent ( ) y x , 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 Journal of Applied Mathematics and Physics 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 func- n Q x 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] [10]. The works of the authors in [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.

Construction of Method
In this work, we assumed that the theoretical solution ( ) y x of (1) can locally be represented by a rational interpolant ( ) To construct an explicit single-step method with (2) for (1), it requires that ( ) r x satisfies the following: Substituting for expressions and simplifying (3) yields 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.

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.

Local Truncation Error
Local Truncation Error: The local truncation error

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 ( )

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 the above indicates that the constructed schemes satisfy the consistency property.

Stability
To get the stability behaviour of the constructed scheme, the scheme is implemented on the standard test problem ( ) , 0 y y Re λ λ ′ = < (12) and the stability polynomial ( )  (13) and the region of absolute stability is seen in Figure 1.

Numerical Examples
The first problem considered in this work is the nonlinear initial value problem ( ) whose theoretical solution is given as 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.

Problem 2
The second test problem considered is given as ( ) The exact Solution is ( ) The logarithm of absolute errors for the solutions obtained is compared with other methods discussed in [12] as given in Figure 4.

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. Journal of Applied Mathematics and Physics