We consider the numerical solution of initial value problem of the form:
In practice, higher order ordinary differential equations of this form , is solved by reducing it to systems of first order differential equation of the form:
then an approximate method is applied to solve the resulting Equation (2) as widely discussed by Fatunla  and Lambert  and Spiegel  . The approach does not utilize additional information associated with the specific ordinary differential equation, and consequently, the oscillatory nature of the solution of the differential equation is always neglected. Thus, it would be more efficient to improve on the numerical method so that higher order ordinary differential equations could be solved without having to reduce to systems of first order as suggested by Chakravati and Worland  , Dahlquist  , sharp and Fine  , and Bun and Vasilsyer  . Actually, considerable attention has been devoted to solving ordinary differential equation of higher order directly without reduction for instance: methods of linear multistep method (LMM) were considered by Lambert and Watson  , Dormand and El-Mikkawy  , El-Mikkawy and El- Desouky  and Awoyemi     . Subsequently, LMM was independently proposed by Kayode  , Onumanyi et al.  and Adesanya et al.  in the predictor-corrector mode, based on collocation method. These authors proposed LMM with continuous coefficients where they adopted Taylor series algorithm to supply the starting values. Also, some notable scholars improve on the predictor-corrector method for solving ordinary differential equations of higher orders, for instance, Jator and Li  proposed five-step and four-step methods respectively in which they adopted a continuous LMM to obtain finite difference method. Moreover, Adesanya  adopted a method of collocation and interpolation to develop a continuous LMM which is evaluated at different grid points to give discret methods that generate independent solutions. Others that adopted block methods include Badmus and Yahaya  . One of the advantages of the method is that it provides direct solution of implicit Linear multistep method without developing separate predictors.
Although some of the aforementioned authors have made use of Taylor series, but little has been said with the use of Taylor series as a major method of implementation. So, Our idea is to use Taylor series algorithm to evaluate
and and calculate , by the
use partial derivative technique. Thus, two-step hybrid methods in the Taylor series mode are developed to solve second order ordinary differential equations directly.
In this section, power series is considered as an approximate solution to the general second order problems:
of the form:
The first and second derivative of (3) are respectively given as:
Combining (2) and (5), we generate the differential system
we develop the hybrid scheme using (3) and (5) as interpolation and collocation equations in this work.
Collocating (6) at selected grid and off-grid points, and interpolating (3) at selected grid and off-grid points, it results into a system of equations:
where, , solving Equations ((7) and (8)), , yield a method ex-
pressed in the form:
Writing these system of equations in matrix form:
Using Gaussian elimination method, the unknown coefficients can be obtained. Putting back into (3) gives (10):
The coefficients , are continuous coefficients obtained using the transformation ,
Then simplifying the continuous , , and taking their first derivatives, we have:
Then, putting gives:
its first derivative
with the order , error constant , and interval of absolute stability Implementation of the method using Taylor series algorithm to evaluate
where, and and,
Finding the partial derivative as follows
2.1. Analysis of the Properties of the Scheme
We shall consider the analysis of the basic properties of our methods which includes the order, the region of absolute stability and the zero stability of the methods.
2.2. Order of Accuracy of the Method
The local truncation error with k-step linear multistep m method which is in line with Lambert (1973), is taken to be linear difference operator defined by
Thus, expanding (21) as Taylor series about point x and comparing coefficients of , the scheme will be of order with error constant
where are the constant coefficients given as:
In line with  , k-step, linear multistep (21) has order p if and , where, is the error constant. Subjecting our schemes to equations 35, it is therefore established that linear multistep scheme is of order , relatively small error constant −0.0069941.
2.3. Consistency of the Scheme
A linear multistep method is consistent if the following conditions are satisfied:
1) The order .
2.4. Zero Stability of the Method
Equation (21) has its first characteristic polynomial to be:
The method is zero stable since they have roots twice.
2.5. Region of Absolute Stability of the Method
In order to establish the region of absolute stability, we apply the boundary locus method as in  . The method implies that
From scheme (21), we have:
Considering the values of for at intervals of gives the region of absolute stability to be
3. Numerical Experiments
We test the accuracy of the proposed scheme on some numerical problems, and the results are compared with existing methods.
Table 1. Results and errors for problem (1).
Table 2. Results and errors for problem (2).
Note: YEX = Yexact, YC = Ycomputed, ERRNew = Error in new method.
The numerical results of the problem is shown in Table 2.
A Linear Multistep method which implements a Taylor’s series algorithm is developed for the direct solution of general second order initial value problems of
ordinary differential equations without reduction to systems of first order differential equation. The derivatives of continuous scheme to any order were computed implementing Taylor’s series algorithm. The accuracy of the method was tested with two test problems, and results were compared with Awoyemi and Kayode  of order (8).
 Lambert, J.D. and Wastson, A. (1976) Symmetric Multistep Method for Periodic Initial Value Problem. Journal of the Institute of Mathematics Application, 18, 189- 202. https://doi.org/10.1093/imamat/18.2.189
 El-Mikkawy, M.E. and El-Desouky, R. (2003) A New Optimized Non FSAL Embedded Runge-Kutta-Nystrom Algorithms of Order 6 and 4 in Six Stages. Applied Mathematics and Computation, 145, 33-43.
 Awoyemi, D.O. and Kayode, S.J. (2005) A Maxima Order Collocation Method for Initial Value Problems of General Second Order Ordinary Differential Equation. Proceeding of the Conference Organized by the National Mathematical Center, Abuja.
 Awoyemi, D.O. (1999) A Class of Continuous Method for General Second Order Initial Value Problems in Ordinary Differential Equations. International Journal of Computer Mathematics, 72, 29-37.
 Awoyemi, D.O. (2001) A New Sixth Order Algorithms for General Second Order Ordinary Differential Equation. International Journal of Computer Mathematics, 77, 117-124.
 Awoyemi, D.O. (2005) An Algorithm Collocation Methods for Direct Solution of Special and General Fourth Order Initial Value Problems of Ordinary Differential Equation. Journal of Nigerian Association of Mathematical Physics, 6, 271-238.
 Onumanyi, P., Sirisons, U.W. and Dauda, Y. (2001) Toward Uniformly Accurate Continuous Finite Difference Approximation of Ordinary Differential Equation. Bayero Journal of Pure and Applied Sciences, 1, 5-8.
 Adesanya, A.O., Anake, T.A. and Udoh, M.O. (2008) Improved Continuous Method for Direct Solution of General Second Order Ordinary Differential Equation. Journal of Nigeria Association of Mathematical Physics, 13, 59-62.
 Jator, S.N. and Li, J. (2009) A Self Stationary Linear Multistep Method for a Direct Solution of General Second Order Initial Value Problem. International Journal of Computer Mathematics, 86, 817-1165.
 Badmus, A.M. and Yahaya, Y.A. (2009) An Accurate Uniform Order 6 Block Method for Direct Solution of General Second Order Ordinary Differential Equations. The Pacific Journal of Science and Technology, 10, 273-277.