HIGHER ORDER GAUSS-LEGENDRE QUADRATURE RUNGE-KUTTA TYPE METHOD FOR SOLVING STIFF AND OSCILLATORY PROBLEMS

Authors

  • Zainab Adamu Nigerian Defence Academy, Kaduna
  • Yakubu Muntaka Lawal Nigerian Defence Academy, Kaduna
  • Aliyu Jaafar Bello Colleges of Education (NCCE), Abuja

DOI:

https://doi.org/10.33003/fjs-2026-1004-4861

Keywords:

Implicit, More efficient, A-stable, Collocation methods, Perturbed Gaussian points, Order and Error constant

Abstract

In this Paper a four (4) stage super convergent implicit Runge-Kutta type methods of order eight (8) have been constructed. [h1] [m2] [m3] This is an improvement of existing lower order four and six Gauss Legendre method for solving Ordinary Differential Equations (ODEs). Legendre polynomials of fourth degree basic function was used to generate the special Gaussian points for the construction of an implicit four points Gauss Legendre Runge-Kutta type methods. Collocation and matrix inversion approach is used to obtain continuous formulas for the method. The continuous formula is evaluated at the special Gaussian points of fourth-degree Legendre polynomials to yield a block discrete scheme which is converted to Runge-Kutta evaluation functions for the integration of Ordinary Differential Equations of ODEs, especially highly stiff and oscillatory problems of first order ODEs which are found in science and engineering models. Experimental problems used, show that the developed methods is A-stable, consistent, more efficient with less error estimate than the lower stage Gauss Legendre Runge-Kutta methods of order 4 and six respectively. 


References

Abdul-Hassan, N. Y., Kadum, Z. J., & Ali, A. H. (2024). An efficient third-order scheme based on Runge–Kutta and Taylor series expansion for solving initial value problems. Algorithms, 17(3), 123.

Adeosun, A. T., Mohd Kasim, A. R., Akolade, M. T., & Tijani, Y. O. (2026). Gauss–Legendre weighted residual method on the stability of Eyring–Powell hybrid nanofluid flow past a shrinking surface with slip conditions. International Journal of Numerical Methods for Heat & Fluid Flow, 36(2), 836-857.

Agam, S.A., and Yahaya, Y.A. (2014). A new three stage implicit Runge-Kutta type method with error estimation for first order ordinary differential equations. Sci. Technol2015,1(1), 19-24.

Agam, S.A (2014), Singly and multiply implicit Runge- kutta Type methods for solutions of first, second and third order Ordinary Differential Equations, Ph.D. thesis (unpublished).

Babaei, M. Optimized Gauss–Legendre–Hermite 2-point (O-GLH-2P) method for nonlinear time-history analysis of structures. Meccanica 59, 305–332 (2024). https://doi.org/10.1007/s11012-023-01752-4

Brugnano, L., Iavernaro, F., & Susca, T. (2010). Numerical comparisons between Gauss-Legendre methods and Hamiltonian BVMs defined over Gauss points. arXiv preprint arXiv:1002.2727.

Butcher, J.C. (1964). Implicit Runge-Kutta processes. Math Comp.18:50-64.

Butcher, J.C. (1966). Convergence of Numerical Solution of Ordinary Differential Equation. Mathematics of Computation, 1-10.

Butcher, J.C and Jackiewicz (1997) Implementation of diagonally implicit multi stage Integration methods of ODEs. SIAM J. Numer. Anal 34:2119-2141.

Butcher, J.C. (2008). Numerical methods for ordinary Differential Equations,

second Edition. John Wiley & Sons Ltd.

Chapra, S.C., & Canale, R.P. (2006), (5th edition). Numerical Methods for

Engineers. New York, NY, USA: McGraw-Hill.

Guo, B., & Wang, Z. (2018). Legendre-Gauss collocation methods for ordinary

differential equations. Springer science-Journal of scientific computing.

Hairer, E. and Wanner, G. (1996). Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problem. New York: Springer-Verlag Berlin Heidelberg GmbH.

Iserles, A. (1996). A First Course in the Numerical Analysis of Differential

Equations. Cambridge: Cambridge University Press. Retrieved from

https://en,wikipedia.org/wiki/Gauss-Legendre_method. (08/05/2023)

Kuntzmann, J. (1961). Neure Entwickhengsweisen der methoden Von Runge und kutta. Z. Angew Math Mech. 41: T28-T31.

Lie, I., and Norsett, S.P. (1989). Super Convergence for multistep collocation

methods. Journal of mathematics computer, 52:65-79.

Moshtaghi, N., & Saadatmandi, A. (2021). Polynomial–Sinc collocation method combined with the Legendre–Gauss quadrature rule for numerical solution of distributed order fractional differential equations. Revista de la Real Academia de Ciencias Exactas, Físicas y Naturales. Serie A. Matemáticas, 115(2), 47.

Onumanyi, P., Awoyemi, D.O., & Sirisena, U.W. (1994). New Linear Multi-Step with Continuous coefficient for first order initial value problems. Journal of Mathematical Society of Nigeria, 13:37-51.

Sasser, E. J. (2005). History of ordinary differential equations. The first hundred years. Cincinnati: Dover publications, Inc.

Teschl, G. (2012. Ordinary Differential Equations and Dynamical Systems. American Mathematical Society. ISBN 978-0-8218-8328-0.

Comparison of Numerical Solution of Problem 1

Downloads

Published

25-02-2026

How to Cite

Adamu, Z., Lawal, Y. M., & Bello, A. J. (2026). HIGHER ORDER GAUSS-LEGENDRE QUADRATURE RUNGE-KUTTA TYPE METHOD FOR SOLVING STIFF AND OSCILLATORY PROBLEMS. FUDMA JOURNAL OF SCIENCES, 10(4), 248-253. https://doi.org/10.33003/fjs-2026-1004-4861

Most read articles by the same author(s)