This work intends to contribute to the growing literature on the direct numerical solution of special and general second-order initial value problems (IVPs) of ordinary differential equations (ODEs) by reformulating an A-stable Runge-Kutta (RK) method as a type of implicit Runge–Kutta–Nystrӧm methods (RKN). Unlike traditional methods that reduce second-order IVPs of ODEs to equivalent systems of first-order IVPs, the approach in this work preserves the original problem structure while also benefiting from the sixth-order accuracy and A-stability of the GLRK method. Using an extended Butcher array, the resulting RKN coefficients are obtained explicitly, ensuring they meet the consistency and order conditions. The linear stability analysis reveals a broad stability region, which is necessary for handling stiff and oscillatory systems. Application of the method to numerical examples shows that it offers better accuracy than most existing methods, while maintaining a similar computational cost. This reformulation strategy paves the way for deriving high-order A-stable RKN methods based on existing RK schemes.
| Published in | Engineering Mathematics (Volume 10, Issue 1) |
| DOI | 10.11648/j.engmath.20261001.12 |
| Page(s) | 14-23 |
| Creative Commons |
This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited. |
| Copyright |
Copyright © The Author(s), 2026. Published by Science Publishing Group |
A-stable, Stiff and Oscillatory Differential Equation, Runge–Kutta–Nystrӧm Methods, Numerical Analysis
| [1] | Ramos, H., Singh, G. A tenth order A-stable two-step hybrid block method for solving initial value problems of ODEs. Applied Mathematics and Computation, 310,(2017), 75-88. |
| [2] | Lambert J. D., Computational Methods in Ordinary Differential Equations, John Wiley, New York, (1973). |
| [3] | Dahlquist G., Convergence and stability in the numerical integration of ordinary differential equations, Math. Scand. 4 (1956) 33-53. |
| [4] | Milne W. E., Numerical Solution of Differential Equations, John Wiley & Sons, New York, (1953). |
| [5] | Henrici P., Discrete variable methods in ordinary differential equations, New York; John Wiley & Sons. (1962). |
| [6] | Ramos, H. Development of a new Runge-Kutta method and its economical implementation. Computational and Mathematical Methods, 1 (2019). (2), e1016. |
| [7] | Zhai, S., Wang, H., Li, J. Exponentially fitted Runge–Kutta–Nystrӧm methods for oscillatory second-order differential equations. Journal of Computational and Applied Mathematics, 388, (2021) 113307. |
| [8] | Rabiei, F., Van Daele, M. Trigonometrically fitted Runge–Kutta–Nystrӧm methods with minimal phase-lag. Applied Numerical Mathematics, 79 (2014), 1–14. |
| [9] | Senu, N., Ibrahim, Z. B., Ismail, F. Zero-dissipative Runge–Kutta–Nystrӧm methods for oscillatory problems. Applied Mathematics and Computation, 216 (2010) (8), 2392-2400. |
| [10] | Olabode , B. T., Momoh, L. A-stable optimal hybrid linear multistep formulas in Runge-Kutta form and their application to dynamical systems. Journal of Mathematical Analysis and Modeling, 5 (2024)(2), 64-80. |
| [11] | Olabode Olabode B. T., Momoh A. L. Continuous Hybrid Multistep Methods with Legendre Basis Function for Direct Treatment of second order stiff Odes. America Journal of Computational and Applied Mathematics, 6 (2016) (2), 38-49, |
| [12] | Ariffin, N., Senu, N., Ismail, F. Diagonally implicit Runge–Kutta–Nystrӧm methods for second-order ordinary differential equations. Mathematical Problems in Engineering, (2013)., Article ID 213157. |
| [13] | Agarwal, P., Ibrahim, I. H. A new type OF hybrid multistep multiderivative formula for solving stiff IVPs. Advances in Difference Equations, (2019)(1). |
| [14] | Abdulsalam, H., Senu, N., Ibrahim, Z. B. Efficient Runge–Kutta–Nystrӧm methods for second-order initial value problems. Applied Numerical Mathematics, 167 (2021), 42–58. |
| [15] | Chen, Z., Wang, Y. Two-derivative Runge–Kutta–Nystrӧm methods for second-order differential equations. Numerical Methods for Partial Differential Equations, 31 (2015)(4), 1234-1256. |
| [16] | Kovalnogov, V., Podhaisky, H. Efficient adaptive Runge–Kutta–Nystrӧm methods with reuse of rejected steps. Numerical Algorithms, 93 (2023)(3), 987-1012. |
| [17] | Zhang, Y., Liu, X. Modified Runge–Kutta–Nystrӧm methods for nonlinear second-order delay boundary value problems. Journal of Computational Mathematics, 41(2023)(2), 245–263. |
| [18] | Rabiei, F., Van Daele, M. Two-derivative Runge–Kutta–Nystrӧm methods for oscillatory problems. Journal of Computational and Applied Mathematics, 236 (2012) (6), 1374–1386. |
| [19] | Imoni S. O., Otunta F. O., Ramamohan T. R. Embedded implicit Runge–Kutta Nystrӧm method for solving second-order differential equations, International Journal of Computer Mathematics, 8311 (2006), 777-784 |
| [20] | Abdelrahim, R., Omar, Z. Direct solution of second-order ordinary differential equation using a single-step hybrid block method of order five. Mathematical and Computational Applications, 21 (2016)(2), 12. Retrieved from |
| [21] | Alkasassbeh, M., Omar, Z. Implicit one-step block hybrid third-derivative method for the direct solution of initial value problems of second-order ordinary differential equations. Journal of Applied Mathematics, (2017), 1-8. Retrieved from |
| [22] | Jator, S. N. Block third derivative method based on trigonometric polynomials for periodic initial-value problems. Afrika Matematika, 27 (2015) (3-4), 365-377. Retrieved from |
APA Style
Lukuman, M. A. (2026). Insight into Reformulation of an A-stable Runge-Kutta (RK) Method as an Implicit Runge-Kutta-Nystrom (RKN) Method for Directly Solving Second-order ODEs. Engineering Mathematics, 10(1), 14-23. https://doi.org/10.11648/j.engmath.20261001.12
ACS Style
Lukuman, M. A. Insight into Reformulation of an A-stable Runge-Kutta (RK) Method as an Implicit Runge-Kutta-Nystrom (RKN) Method for Directly Solving Second-order ODEs. Eng. Math. 2026, 10(1), 14-23. doi: 10.11648/j.engmath.20261001.12
Share This Article
Twitter
Linked In
Facebook
Copy
Download@article{10.11648/j.engmath.20261001.12,
author = {Momoh Adelegan Lukuman},
title = {Insight into Reformulation of an A-stable Runge-Kutta (RK) Method as an Implicit Runge-Kutta-Nystrom (RKN) Method for Directly Solving Second-order ODEs
},
journal = {Engineering Mathematics},
volume = {10},
number = {1},
pages = {14-23},
doi = {10.11648/j.engmath.20261001.12},
url = {https://doi.org/10.11648/j.engmath.20261001.12},
eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.engmath.20261001.12},
abstract = {This work intends to contribute to the growing literature on the direct numerical solution of special and general second-order initial value problems (IVPs) of ordinary differential equations (ODEs) by reformulating an A-stable Runge-Kutta (RK) method as a type of implicit Runge–Kutta–Nystrӧm methods (RKN). Unlike traditional methods that reduce second-order IVPs of ODEs to equivalent systems of first-order IVPs, the approach in this work preserves the original problem structure while also benefiting from the sixth-order accuracy and A-stability of the GLRK method. Using an extended Butcher array, the resulting RKN coefficients are obtained explicitly, ensuring they meet the consistency and order conditions. The linear stability analysis reveals a broad stability region, which is necessary for handling stiff and oscillatory systems. Application of the method to numerical examples shows that it offers better accuracy than most existing methods, while maintaining a similar computational cost. This reformulation strategy paves the way for deriving high-order A-stable RKN methods based on existing RK schemes.
},
year = {2026}
}
TY - JOUR T1 - Insight into Reformulation of an A-stable Runge-Kutta (RK) Method as an Implicit Runge-Kutta-Nystrom (RKN) Method for Directly Solving Second-order ODEs AU - Momoh Adelegan Lukuman Y1 - 2026/04/08 PY - 2026 N1 - https://doi.org/10.11648/j.engmath.20261001.12 DO - 10.11648/j.engmath.20261001.12 T2 - Engineering Mathematics JF - Engineering Mathematics JO - Engineering Mathematics SP - 14 EP - 23 PB - Science Publishing Group SN - 2640-088X UR - https://doi.org/10.11648/j.engmath.20261001.12 AB - This work intends to contribute to the growing literature on the direct numerical solution of special and general second-order initial value problems (IVPs) of ordinary differential equations (ODEs) by reformulating an A-stable Runge-Kutta (RK) method as a type of implicit Runge–Kutta–Nystrӧm methods (RKN). Unlike traditional methods that reduce second-order IVPs of ODEs to equivalent systems of first-order IVPs, the approach in this work preserves the original problem structure while also benefiting from the sixth-order accuracy and A-stability of the GLRK method. Using an extended Butcher array, the resulting RKN coefficients are obtained explicitly, ensuring they meet the consistency and order conditions. The linear stability analysis reveals a broad stability region, which is necessary for handling stiff and oscillatory systems. Application of the method to numerical examples shows that it offers better accuracy than most existing methods, while maintaining a similar computational cost. This reformulation strategy paves the way for deriving high-order A-stable RKN methods based on existing RK schemes. VL - 10 IS - 1 ER -
Department of Mathematical Sciences/ School of Physical Sciences, Federal University of Technology, Akure, Nigeria
