Development of a Numerical Scheme

In this paper, we developed a new numerical scheme which aimed to solve some initial value problems of ordinary differential equations. The full breakdown of this new numerical scheme derivation is presented. While in our subsequent research, we shall fully examine the characteristics of the scheme such as consistency, convergence and stability. Also, the implementation of this new numerical scheme shall be worked-on and comparison shall also be made with some existing methods.


Introduction
Many numerical analysts such as: S. O. Fatunla [1], E. A. Ibijola [2] [3], R. B. Ogunrinde [4] and even A. A. Obayomi [5] and so on, have developed schemes for the solution of some initial value problem of ordinary differential equations.The efficiency of all these contributed effort from this numerical analyst in numerical analysis had been measured and tested for their stability, accuracy, convergence and consistency properties.The accuracy properties of different methods are usually compared by considering the order of convergence as well as the truncation error coefficients of the various methods (C.F. Tischer, 1984).From literatures, this shows that so many methods which are suitable for solving some sets of initial value problems (ivps) in ordinary differential equations (ODEs) must have all the mentioned characteristics.
Ogunrinde, R. B. [4], developed a scheme in which standard finite difference schemes were developed.Similarly, Obayomi, A. A. [5] [6], also worked on some approximation techniques which was used to derive qualitatively stable non-standard finite difference schemes.
In this paper, a new numerical scheme was developed with the above mentioned characteristics in mind to solve some initial value problems of ordinary differential equations which was based on the local representation of the theoretical solution to initial value problem of the form: , where 0 a , 1 a , 2 a and b are real undetermined coefficients.

Derivation of the New Scheme
Suppose we have the initial value problem: Let us assume that the theoretical solution ( ) y x to (1) can be locally represented in the interval ( ) x x + , 0 n ≥ by the interpolating polynomial function:

Conclusions
We aim to develop a new numerical scheme which can favourably agree with the existing ones for solving some initial value problems of ordinary differential equations.Clearly, this paper has been able to show the development of the new numerical scheme as proposed.
In our subsequent research, we shall pay more attention on the implementation of this new scheme to solve some initial value problems (ivp) of the form (1) and also compare the results with the existing methods and thereafter we examine the characteristics properties such as the stability, convergence, accuracy and consistency )