A Three-Stage Multiderivative Explicit Runge-Kutta Method

In recent years, the derivation of Runge-Kutta methods with higher derivatives has been on the increase. In this paper, we present a new class of three stage Runge-Kutta method with first and second derivatives. The consistency and stability of the method is analyzed. Numerical examples with excellent results are shown to verify the accuracy of the proposed method compared with some existing methods.


Introduction
The derivation of Runge-Kutta schemes involving higher derivatives is now on the increase.Traditionally, given an initial value problem (IVP), classical explicit Runge-Kutta methods are derived with the intention of performing multiple evaluations of   f y in each internal stage for a given accuracy.Recently, Akanbi [1,2] derived multi-derivative explicit Runge-Kutta method involving up to second derivative.Goeken and Johnson [3] also derived explicit Runge-Kutta schemes of stages up to four with the first derivative of   f y .However, the new scheme is derived with the notion of incorporating higher order derivatives of   f y up to the second derivative.The cost of internal stage evaluations is reduced greatly and there is an appreciable improvement on the attainable order of accuracy of the method.

Derivation of the Proposed Scheme
The general form of a single step method for solving the Initial Value Problem (IVP) is defined as where and for the autonomous case of (1), in which The proposed scheme of this paper is of the form where   Expanding 2 and 3 in Taylor's series and substituting the result into (5), the coefficients of the powers of are then compared with that of (3) to obtain the following system of equations: x y h  is obtained using the Taylor's series expansion of an arbitrary function: Solving the above system of equations, we have the set of solutions in Table 1.
The above solution set gives rise to a family of 3-stage multi-derivative explicit Runge-Kutta schemes.The proposed scheme denoted by 3sMERK above is thus given by  

Existence and Uniqueness of Solution
The properties of the incremental function  3sMERK , ; n n x y h  of the newly derived scheme are in general, very crucial to its stability and convergence characteristics [1,[4][5][6][7][8].
Let  ,  f x y , where , be defined and continuous for all x y in the region defined by , where a, b are finite, and let there exist a constant L such that  .The requirement (7) is known as the Lipschitz's condition, and the constant is a Lipschitz's constants [6,7,[9][10][11].We shall assume that the hypothesis of this theorem is satisfied by the IVP (1).The following lemma will be useful for establishing the aforementioned characteristics.
n be a set of real numbers.If there exist finite constants  and  such that Then, from (8) i j  implies that 1 .

Accuracy and Stability
Usually, during the implementation of a computational scheme, errors are generated.The magnitude of the error determines how accurate and stable a scheme is.For instance, if the magnitude of the error is sufficiently small, the computational results would be accurate.However, if the magnitude of the error becomes so large, it can make the method unstable.The sources of error for these schemes and their principal error functions are discussed in Butcher [5,6], Fatunla [7] and Lambert [9,10].The following theorem guarantees the stability of the 3sMERK methods.
Applying triangle inequality and (13), we have If we assume , and , then Lemma 3.2 implies that . This implies the stability of the 3sMERK scheme.

Numerical Experiments
The proposed 3sMERK scheme ( 6) is applied to the two IVPs below and the results obtained are compared with the standard 3-stage methods of Runge-Kutta (Heun's) [7,9,10] and that of Goeken and Johnson [3] stated in ( 16) and ( 17) respectively.  .

Conclusions
The results generated by the proposed scheme in this paper when applied to the problems above, evidently proved the extent of accuracy of the scheme.Tables 2  and 3 above show the absolute error associated with the schemes for the test problems with the variation of the step length.The computations above clearly show the accuracy of the method.The standard Heun's (third order) method grows faster in error than the method of Goeken and the newly derived scheme.However, 3sMERK performed best among the three methods.
Based on the two problems solved above, it follows that the scheme is quite efficient.We therefore conclude that the 3sMERK method proposed is reliable, stable and with high accuracy.

)
holds every     , , , x y x y D  , then for any n   , there exist a unique solution   y x of the problem (1), where   y x is continuous and differentiable for all   , x y D

Table 1 . Examples of three-stage MERK methods.
Theorem 3.3.Suppose the IVP (1) satisfies the hypotheses of Theorem .1, then the new 3sMERK algorithm is stable.3