High Order Continuous Extended Linear Multistep Methods for Approximating System of ODEs

  • I. M. Esuabana Department of Mathematics, University of Calabar, Nigeria
  • S. E. Ogunfeyitimi Department of Computer Sciences, Wellspring University, Benin City, Nigeria
Keywords: linear multistep formulas, extended Enright's methods, A-stable, ODEs


A class of high-order continuous extended linear multistep methods (HOCELMs) is proposed for solving systems of ordinary differential equations (ODEs). These continuous schemes are obtained through multistep collocation at various points to create a single block method with higher dimensions. This class of schemes consists of A-stable methods with a maximum order of $p\leq14$, capable of yielding moderately accurate results for equations with several eigenvalues of the Jacobians located close to the imaginary axis. The results obtained from numerical experiments indicate that these schemes show great promise and competitiveness when compared to existing methods in the literature.


Cash, J. R. (2003). Review paper. Efficient numerical methods for the solution of stiff initial-value problems and differential algebraic equations. Proc. R. Soc. London. Ser. Math. Phys. Eng. Sci., 459 (2032), 797-815. https://doi.org/10.1098/rspa.2003.1130

Dahlquist, G. (1963). A special stability problem for linear multistep methods. BIT, 3, 27-43. https://doi.org/10.1007/BF01963532

Curtis, C. F., & Hirschfelder, J. O. (1952). Integration of stiff equations. National Academy of Sciences, 38, 235-243. https://doi.org/10.1073/pnas.38.3.235

Bickart, T. A., & Rubin, W. B. (1974). Composite multistep methods and stiff stability. In R. A. Willoughby (Ed.), Stiff Differential Systems (pp. 293-307). Plenum Press, New York. https://doi.org/10.1007/978-1-4684-2100-2_2

Enright, W. H. (1974). Second derivative multistep methods for stiff ordinary differential equations. SIAM. J. Numer. Anal., 11, 321-331. https://doi.org/10.1137/0711029

Obrechkoff, N. (190). Neue Quadraturformeln. Abh Preuss Akad Wiss Math Nat K14.

Cash, J. R. (1981). Second derivative extended backward differentiation formula for the numerical integration of stiff system. SIAM. J. Numer. Anal., 18, 21-36. https://doi.org/10.1137/0718003

Jia-Xiang, X., & Jiao-Xun, K. (1998). A class of DBDF methods with the derivative modifying term. J Comput. Math., 6, 7-13.

Ngwane, F. F., & Jator, S. N. (2015). A family of trigonometrically fitted Enright second derivative methods for stiff and oscillatory initial-value problem. J. Appl. Math.. https://doi.org/10.1155/2015/343295

Moradi, A., Abdi, A., & Hojjati, G. (2022). Implicit-explicit second derivative general linear methods with strong stability preserving explicit part. Appl. Numer. Mat., 181, 23-45. https://doi.org/10.1016/j.apnum.2022.05.012

Gragg, W., & Stetter, H. J. (1964). Generalized multistep predictor-corrector methods. J. Assoc. Comput. Mach., 11, 188-209. https://doi.org/10.1145/321217.321223

Butcher, J. C. (1965). A modified multistep methods for the numerical integration of ordinary differential equations. J. ACM, 12, 125-135. https://doi.org/10.1145/321250.321261

Gear, C. W. (1965). Hybrid methods for initial value problems in ordinary differential equations. SIAM. J. Numer. Anal., 2, 69-86. https://doi.org/10.1137/0702006

Vigo-Aguiar, J., & Ramos, H. (2006). A new eighth-order A-stable method for solving differential systems arising in chemical reactions. J. Math. Chem., 40, 71-83. https://doi.org/10.1007/s10910-006-9121-x

Selva, M., Arevalo, C., & Fuherer, C. (2002). A collocation formulation of multistep methods for variable step-size extensions. Appl. Numer. Math., 42, 5-16. https://doi.org/10.1016/S0168-9274(01)00138-6

Okuonghae, R. I., & Ikhile, M. N. O. (2013). A class of hybrid linear multistep methods with A(a)-stability properties for stiff IVPs in ODEs. J. Numer. Math., 21, 157-172. https://doi.org/10.1515/jnum-2013-0006

Aiguobansimwin, I. B., & Okuonghae, R. I. (2019). A class of two-derivative two-step Runge-Kutta Methods for non-stiff ODEs. Appl Maths. Hindawi. https://doi.org/10.1155/2019/2459809

Esuabana, I. M., & Ekoro, S. E. (2017). Hybrid linear multistep methods with nested hybrid predictors for solving linear and non-linear initial value problems in ordinary differential equations. Mathematical Theory and Modeling, 11, 77-88.

Esuabana, I. M., & Ekoro, S. E. (2018). Derivation and implementation of new family of second derivative hybrid linear multistep methods for stiff ordinary differential equations. Global Journal of Mathematics, 2, 821-828.

Gupta, G. K. (1978). Implementation second-derivative methods using Nordsieck polynomial representation. Math. Comp., 332, 13-18. https://doi.org/10.1090/S0025-5718-1978-0478630-7

Golik, W. L., Amodio, P., & Mazzia, F. (1995). Variable-step boundary-value methods based on reverse Adams schemes and their grid distribution. Appl. Numer. Math., 18, 5-21. https://doi.org/10.1016/0168-9274(95)00044-U

Brugnano, L., & Trigiante, D. (1996). Convergence and stability of boundary-value methods for ordinary differential equations. J. Comput. Math., 66, 97-109. https://doi.org/10.1016/0377-0427(95)00166-2

Brugnano, L., & Trigiante, D. (1998). Solving differential problems by multistep initial and boundary-value methods. Gordon and Breach Science Publishers.

Aceto, L., & Trigiante, D. (2002). On the A-stable method in the GBDF class. Nonlinear Analysis Real World Appl., 3, 9-23. https://doi.org/10.1016/S1468-1218(01)00009-8

Ogunfeyitimi, S. E., & Ikhile, M. N. O. (2019). Second derivative generalized extended backward differentiation formulas for stiff problems. J. Korean Soc. Ind. Appl. Math., 23, 179-202.

Ogunfeyitimi, S. E., & Ikhile, M. Ν. Ο. (2020). Generalized second derivative linear multistep methods based on the methods of Enright. Int. J. Appl. and Comput. Math., 6, 76. https://doi.org/10.1007/s40819-020-00827-0

Akinfenwa, O. A., Jator, S. N., & Yao, N. M. (2013). Continuous block backward differentiation formula for solving stiff ordinary differential equation. J. Comput. Appl. Math. Appl., 65, 996-1005. https://doi.org/10.1016/j.camwa.2012.03.111

Jator, S. N., & Sahi, R. K. (2010). Boundary-value technique for initial value problems based on Adams-type second derivative methods. Int. J. Math. Educ. Sci. Educ., 1-8.

Ramos, H. (2017). An optimized two-step hybrid block method for solving first-order initial-value problems in ODEs Differential Geometry. Dynamical Systems, 19, 107-118.

Ramos, H., & Singh, G. (2017). A tenth order A-stable two-step hybrid block method for solving initial value problems of ODEs. Applied Mathematics and Computation, 310, 75-88. https://doi.org/10.1016/j.amc.2017.04.020

Ramos, H., & Popescu, P. (2018). How many k-step linear block methods exist and which of them is the most efficient and simplest one? Applied Mathematics and Computation, 316, 296-309. https://doi.org/10.1016/j.amc.2017.08.036

Butcher, J. C. (2003). Numerical methods for ordinary differential equations. Wiley, England. https://doi.org/10.1002/0470868279

Cash, J. R. (2000). Modified extended backward differentiation formulae for the numerical solution of stiff initial value problems in ODEs and DAEs. Journal of Computational and Applied Mathematics, 125, 117-130. https://doi.org/10.1016/S0377-0427(00)00463-5

Milner, W. E. (1953). Numerical solution of differential equations. John Wiley and Sons, New York.

Sarafyan, D. (1965). Multistep method for the numerical solution of ordinary differential equations made self-starting. Wisconsin Univ Madison Mathematics research centre Tech Rep 495.

Sommeijer, B., Couzy, W., & Houwen, P. (1989). A-stable parallel block methods. Report NM-R8919. Center for Mathematics and Computer Science, Amsterdam.

Axelson, A. O. H. (1969). A class of A-stable methods. BIT, 9, 185-197. https://doi.org/10.1007/BF01946812

Chartier, P. (1993). L-stable parallel one-block methods for ordinary differential equations. Technical Report 1650 INRIA.

Fatunla, S. (1990). Block methods for second-order ODEs. International Journal of Computer Mathematics, 14, 55-56. https://doi.org/10.1080/00207169108804026

Iavernaro, F., & Mazzia, F. (1999). Block-boundary value methods for the solution of ordinary differential equations. Journal of Scientific Computing, 21, 323-339. https://doi.org/10.1137/S1064827597325785

Brugnano, L., & Trigiante, D. (2000). Block implicit methods for ODEs. In Recent Trends in Numerical Analysis. Nova Science, New York, 81-105.

Conte, S. D., & de Boor, C. (1981). Elementary numerical analysis, an algorithmic approach, third ed. McGraw-Hill, Tokyo, Japan.

Rosser, J. D. (1967). A Runge-Kutta for all seasons. SIAM, 9, 417-452. https://doi.org/10.1137/1009069

Motsa, S. S. (2022). Hybrid block methods for IVPs using Mathematica. In Proceeding of the 14th Annual Workshop on Computational Mathematical and Modelling (pp. 4-8). University of KwaZulu-Natal, Pietermaritzburg Campus, Durban, South Africa.

Fatunla, S. O. (1989). Numerical methods for initial value problems in ordinary differential equations. Academic Press Inc, London. https://doi.org/10.1016/B978-0-12-249930-2.50012-6

Lambert, J. D. (1991). Numerical methods for ordinary differential equations. Wiley, New York.

Henrici, P. (1962). Discrete variable methods in ODEs. John Wiley, New York.

Ogunfeyitimi, S. E., & Ikhile, M. N. O. (2021). Multi-block generalized Adams-type integration methods for differential Algebraic equations. International Journal of Applied Computational Mathematics, 7, 1-29. https://doi.org/10.1007/s40819-021-01135-x

Akinfenwa, O. A., & Jator, S. N. (2015). Extended continuous block backward Differentiation formula for stiff systems. Fasciculi Mathematici, 55, 5-18. https://doi.org/10.1515/fascmath-2015-0010

Ismail, G., & Ibrahim, I. (1998). A new higher-order effective P-C methods for stiff systems. Journal of Mathematics and Computers in Simulation, 47, 541-552. https://doi.org/10.1016/S0378-4754(98)00136-0

Fotta, A. U., & Alabi, T. J. (2015). Block method with one hybrid point for the solution of first-order initial value problems of ordinary differential equations. International Journal of Pure and Applied Mathematics, 103, 511-521. https://doi.org/10.12732/ijpam.v103i3.12

Rufai, M. A., Duromola, M. K., & Ganiyu, A. A. (2016). Derivation of one-sixth hybrid block method for solving general first order ordinary differential equations. IOSR-JM, 12, 20-27. https://doi.org/10.9790/5728-1205022027

Wu, X., & Xia, J. (2001). Two low accuracy methods for stiff systems. Applied Mathematics and Computation, 123, 141-153. https://doi.org/10.1016/S0096-3003(00)00010-2

Ehigie, J., Jator, S., Sofoluwe, A. B., & Okunuga, S. A. (2014). Boundary-value technique for initial value problems with continuous second derivative multistep method of Enright. Computers and Applied Mathematics, 33, 81-93. https://doi.org/10.1007/s40314-013-0044-4

Akinfenwa, O. A., Jator, S. N., & Yoa, N. M. (2011). An eight order backward differentiation formula with continuous coefficients for stiff ordinary differential equations. World Academy of Science, Engineering and Technology, 74, 848-853.

Adid, A., & Conte, D. (2020). Implementation of second derivative general linear methods. Calcolo, 20. https://doi.org/10.1007/s10092-020-00370-w

Ramos, H., & Rufai, M. A. (2023). One-step method with three intermediate points in a variable step-size mode for stiff differential systems. Journal of Mathematical Chemistry, 61, 673-688. https://doi.org/10.1007/s10910-022-01427-7

Suleiman, M. B., Musa, H., Ismail, F., & Senu, N. (2013). A new variable step size block backward differentiation formula for solving stiff initial value problems. International Journal of Computer Mathematics, 90, 2391-2408. https://doi.org/10.1080/00207160.2013.776677

How to Cite
Esuabana, I. M., & Ogunfeyitimi, S. E. (2024). High Order Continuous Extended Linear Multistep Methods for Approximating System of ODEs. Earthline Journal of Mathematical Sciences, 14(3), 501-533. https://doi.org/10.34198/ejms.14324.501533