Many of the numerical methods for obtaining the solution of the first order ordinary differential equation
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  . This is evident as Runge-Kutta type methods, Obrechkoff methods and general linear multistep methods usually produce very poor solutions around singularity points    . The authors in  were the first to develop quadrature formulas based on rational interpolating functions. On the other hand, the rational interpolation schemes proposed in , was seen to be effective in the neighbourhood of the singularity and even beyond as reiterated in . The work of authors in  was modified by  with a replacement of the general rational function Luke et al. replaced the general rational function of  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     . The works of the authors in       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
To construct an explicit single-step method with (2) for (1), it requires that satisfies the following:
Substituting for expressions and simplifying (3) yields
Eliminating the undetermined coefficients , , , , and in (4) results in
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
is given as
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
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.
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.
To get the stability behaviour of the constructed scheme, the scheme is implemented on the standard test problem
and the stability polynomial , where is obtained. The stability function of (5) is obtained as
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
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  and the derivative-free methods proposed in  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      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.
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  as given in Figure 4.
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.
 Luke, Y.L., Fair, W. and Wimp, J. (1975) Predictor Corrector Formulas Based on Rational Interpolants. Computers & Mathematics with Applications, 1, 3-12.
 Lambert, J.D. and Shaw, B. (1965) On the Numerical Solution of y’ = f(x,y) by a Class of Formulae Based on Rational y’ Approximation. Mathematics of Computation, 19, 456-462.
 Fatunla, S.O. (1986) Numerical Treatment of Singular Initial Value Problems. An International Journal of Computers and Mathematics with Applications, 128, 1109-1115.
 Fatunla, S.O. (1982) Nonlinear Multistep Methods for Initial Value Problems. An International Journal of Computers and Mathematics with Applications, 8, 231-239.
 Tasneem, A., Asif, A.S. and Sania, Q. (2018) Development of a Nonlinear Hybrid Numerical Method. Advances in Differential Equations and Control Process, 19, 275-285.
 Ying, T.Y., Omar, Z. and Mansor, K.H. (2014) Modified Exponential-Rational Methods for the Numerical Solution of First Order Initial Value Problems. Sains Malaysiana, 43, 1951-1959.