A four-step predictor-corrector method with a higher derivative is presented for the direct solution of second-order initial value problems (IVPs). In deriving the algorithm, power series and Bernstein polynomial in combine was used as an approximate function. The approximate solution is interpolated while its second and third derivatives are collocated at the given interval. The main predictor was developed using the same approach and is of the equal order as the corrector. The developed method (Corrector) was found to be consistent, zero stable, and convergence. Numerical examples are tested to demostrate the accuracy of the method.
 Omolehin, J.O., Ibiejugba, M.A., Alabi, M.O., and Evans, D.J, A new class of adams-bashforth schemes for ODEs. Int. Journal of Comput. Math. 80 (2010) 629-638. DOI:10.1080/0020716021000023051
 Kayode, S.J, A Zero Stable Optimal Order Method for Direct Solution of Second Order Differential Equations. Journal of Mathematics Statistics, Vol. 6(3), (2010) 367-371
 Adeyeye, O. and Omar, Z, Maximal Order Block Method For the solution of Second Order Ordinary Differential Equations. IAENG International Journal of Applied Mathematics Vol. 46 (4), (2016) IJAM-46-4-03
Financed by the National Centre for Research and Development under grant No. SP/I/1/77065/10 by the strategic scientific research and experimental development program:
SYNAT - “Interdisciplinary System for Interactive Scientific and Scientific-Technical Information”.