Volume 2, Issue 1, ICCK Journal of Applied Mathematics
Volume 2, Issue 1, 2026
Submit Manuscript Edit a Special Issue
Article QR Code
Article QR Code
Scan the QR code for reading
Popular articles
ICCK Journal of Applied Mathematics, Volume 2, Issue 1, 2026: 87-110

Open Access | Research Article | 04 March 2026
Time Dependent Algorithm using Modified Three-Point Superclass of Block Backward Differentiation Formula for Solving Stiff Ordinary Differential Equations
1 Department of Mathematics and Statistics, Al-Qalam University Katsina, Katsina State, Nigeria
2 Department of Mathematics, Umaru Musa Yar’adua University Katsina, Katsina State, Nigeria
* Corresponding Author: Buhari Alhassan, [email protected]
ARK: ark:/57805/jam.2026.213905
Received: 15 January 2026, Accepted: 03 February 2026, Published: 04 March 2026  
Abstract
In this paper, we present a modified three-point superclass of Block Backward Differentiation Formula (BBDF) for the efficient numerical solution of stiff systems of ordinary differential equations (ODEs). The principal enhancement of this work is a structural modification of the classical BBDF that forms a new parameterized superclass of methods, leading to improved stability and reduced error constants compared with the standard three-point BBDFs. The proposed scheme is formulated as a fully implicit block method capable of simultaneously producing three solution approximations within each integration step. A detailed theoretical analysis is conducted to establish the order of accuracy, consistency, zero-stability, and convergence of the method. Stability analysis based on Dahlquist test equation confirms that the scheme is A-stable and suitable for a wide range of stiff ODEs. The nonlinear systems arising from the implicit formulation are efficiently solved using Newton iteration. Numerical experiments on benchmark stiff ODE problems, including systems with rapidly decaying transient components, demonstrate that the modified three-point superclass BBDF achieves higher accuracy and lower computational cost when compared with existing block implicit methods, such as the standard BBDF and diagonally implicit three-point block BDF schemes. Overall, the proposed method provides a robust and computationally efficient alternative for the numerical integration of stiff ODEs, with potential applications in chemical kinetics, control theory, and biological modeling.

Graphical Abstract
Time Dependent Algorithm using Modified Three-Point Superclass of Block Backward Differentiation Formula for Solving Stiff Ordinary Differential Equations

Keywords
block backward differentiation formula
stiff ODEs
convergence
consistency
a-stability

Data Availability Statement
Data will be made available on request.

Funding
This work was supported without any funding.

Conflicts of Interest
The authors declare no conflicts of interest.

AI Use Statement
The authors declare that generative AI (ChatGPT, version 5.2) was used for language assistance in the preparation of this manuscript. All content was reviewed and approved by the authors.

Ethical Approval and Consent to Participate
Not applicable.

References
  1. Bavarian, M., Soroush, M., Kevrekidis, I. G., & Benziger, J. B. (2010). Mathematical modeling, steady-state and dynamic behavior, and control of fuel cells: A review. Industrial & engineering chemistry research, 49(17), 7922-7950.
    [CrossRef]   [Google Scholar]
  2. Dameris, M., Grewe, V., Ponater, M., Deckert, R., Eyring, V., Mager, F., ... & Giorgetta, M. A. (2005). Long-term changes and variability in a transient simulation with a chemistry-climate model employing realistic forcing. Atmospheric Chemistry and Physics, 5(8), 2121-2145.
    [CrossRef]   [Google Scholar]
  3. Zhao, Z., Tan, D., & Li, K. (2019). Transient behaviors of multiscale megawatt power electronics systems—Part I: Characteristics and analysis. IEEE Journal of Emerging and Selected Topics in Power Electronics, 7(1), 7-17.
    [CrossRef]   [Google Scholar]
  4. Subedi, S., Rauniyar, M., Ishaq, S., Hansen, T. M., Tonkoski, R., Shirazi, M., ... & Cicilio, P. (2021). Review of methods to accelerate electromagnetic transient simulation of power systems. IEEE access, 9, 89714-89731.
    [CrossRef]   [Google Scholar]
  5. Lambert, J. D. (1973). Computational Methods in Ordinary Differential Equations. Chichester, New York: Wiley.
    [Google Scholar]
  6. Ibrahim, Z. B., Ijam, H. M., Aksah, S. J., & Abd Rasid, N. (2024). Enhancing accuracy and efficiency in stiff ODE integration using variable step diagonal BBDF approaches. Malaysian Journal of Fundamental and Applied Sciences, 20(5), 1083-1100.
    [CrossRef]   [Google Scholar]
  7. Gear, C. W. (1971). Numerical Initial Value Problems in Ordinary Differential Equations. London: Prentice-Hall.
    [Google Scholar]
  8. Wanner, G., & Hairer, E. (1996). Solving ordinary differential equations II (Vol. 375, p. 98). New York: Springer Berlin Heidelberg.
    [Google Scholar]
  9. Alexander, R. (1977). Diagonally implicit Runge–Kutta methods for stiff ODE’s. SIAM Journal on Numerical analysis, 14(6), 1006-1021.
    [CrossRef]   [Google Scholar]
  10. Mahnken, R., & Westermann, H. (2024). Construction of A-stable explicit last-stage diagonal implicit Runge–Kutta (ELDIRK) methods. Computational Mechanics, 74(3), 485-507.
    [CrossRef]   [Google Scholar]
  11. Liu, X., Xia, Y., Luo, H., & Xuan, L. (2016). A comparative study of Rosenbrock-type and implicit Runge-Kutta time integration for discontinuous Galerkin method for unsteady 3D compressible Navier-Stokes equations. Communications in Computational Physics, 20(4), 1016-1044.
    [CrossRef]   [Google Scholar]
  12. Rothe, S., Hamkar, A. W., Quint, K. J., & Hartmann, S. (2012). Comparison of diagonal-implicit, linear-implicit and half-explicit Runge–Kutta methods in non-linear finite element analyses. Archive of Applied Mechanics, 82(8), 1057-1074.
    [CrossRef]   [Google Scholar]
  13. Ferracina, L., & Spijker, M. N. (2008). Strong stability of singly-diagonally-implicit Runge–Kutta methods. Applied Numerical Mathematics, 58(11), 1675-1686.
    [CrossRef]   [Google Scholar]
  14. Hazizah, M. I., Saufianin, J. A., Ahmad, F. N. R., Norshakila, A. R., Athraa, A., Nurul, H. M. A., & Fatimah, H. (2024). Numerical Solutions of Stiff Chemical Reaction Problems using Hybrid Block Backward Differentiation Formula. Journal of Advanced Research, 25(1), 100-115.
    [Google Scholar]
  15. Azhari, N. A. I., Abd Rasid, N., Ijam, H. M., Arsad, R., Kamal, M. H. A., Salih, N. M., & Aksah, S. J. (2025). Designing an implicit Block Backward Differentiation Formula (BBDF) for stiff ordinary differential equations. Data Analytics and Applied Mathematics (DAAM), 16-30.
    [CrossRef]   [Google Scholar]
  16. Ibrahim, Z. B., Othman, K. I., & Suleiman, M. (2007). Implicit r-point block backward differentiation formula for solving first-order stiff ODEs. Applied Mathematics and Computation, 186(1), 558-565.
    [CrossRef]   [Google Scholar]
  17. Yakubu, S. D., & Sibanda, P. (2025). A novel one-step optimized hybrid block method for solving general second-order ordinary differential equations. Recent Advances in Natural Sciences, 135-135.
    [CrossRef]   [Google Scholar]
  18. Gupta, G. K., Sacks-Davis, R., & Tescher, P. E. (1985). A review of recent developments in solving ODEs. ACM Computing Surveys (CSUR), 17(1), 5-47.
    [CrossRef]   [Google Scholar]
  19. Byrne, G. D., & Hindmarsh, A. C. (1987). Stiff ODE solvers: A review of current and coming attractions. Journal of Computational physics, 70(1), 1-62.
    [CrossRef]   [Google Scholar]
  20. Singla, R., Singh, G., Ramos, H., & Kanwar, V. (2023). Development of a Higher-Order $\mathcal{A$-Stable Block Approach with Symmetric Hybrid Points and an Adaptive Step-Size Strategy for Integrating Differential Systems Efficiently. Symmetry, 15(9), 1635.
    [CrossRef]   [Google Scholar]
  21. Söderlind, G. (2002). Automatic control and adaptive time-stepping. Numerical Algorithms, 31(1), 281-310.
    [CrossRef]   [Google Scholar]
  22. Tao, M., Owhadi, H., & Marsden, J. E. (2010). Nonintrusive and structure preserving multiscale integration of stiff ODEs, SDEs, and Hamiltonian systems with hidden slow dynamics via flow averaging. Multiscale Modeling & Simulation, 8(4), 1269-1324.
    [Google Scholar]
  23. Tao, M., Owhadi, H., & Marsden, J. E. (2010). Nonintrusive and structure preserving multiscale integration of stiff ODEs, SDEs, and Hamiltonian systems with hidden slow dynamics via flow averaging. Multiscale Modeling & Simulation, 8(4), 1269-1324.
    [CrossRef]   [Google Scholar]
  24. Ascher, U. M., & Petzold, L. R. (1993). Stability of computational methods for constrained dynamics systems. SIAM Journal on Scientific Computing, 14(1), 95-120.
    [CrossRef]   [Google Scholar]
  25. Arnold, M., Burgermeister, B., & Eichberger, A. (2007). Linearly implicit time integration methods in real-time applications: DAEs and stiff ODEs. Multibody System Dynamics, 17(2), 99-117.
    [CrossRef]   [Google Scholar]
  26. Musa, H., Suleiman, M. B., Ismail, F., Senu, N., & Ibrahim, Z. B. (2013). An accurate block solver for stiff initial value problems. International Scholarly Research Notices, 2013(1), 567451.
    [CrossRef]   [Google Scholar]
  27. Chu, Y., & You, F. (2015). Model-based integration of control and operations: Overview, challenges, advances, and opportunities. Computers & Chemical Engineering, 83, 2-20.
    [CrossRef]   [Google Scholar]
  28. Geart, C. W., & Watanabe, D. S. (1974). Stability and convergence of variable order multistep methods. SIAM Journal on Numerical Analysis, 11(5), 1044-1058.
    [CrossRef]   [Google Scholar]
  29. 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(11), 2391-2408.
    [CrossRef]   [Google Scholar]
  30. Fronk, C., & Petzold, L. (2024). Training stiff neural ordinary differential equations with implicit single-step methods. Chaos: An Interdisciplinary Journal of Nonlinear Science, 34(12).
    [CrossRef]   [Google Scholar]
  31. Martín-Vaquero, J., & Vigo-Aguiar, J. (2007). Adapted BDF algorithms: higher-order methods and their stability. Journal of Scientific Computing, 32(2), 287-313.
    [CrossRef]   [Google Scholar]
  32. Vijitha, H. Y. K. (1985). Variable stepsize variable order multistep methods for stiff ordinary differential equations (Doctoral dissertation). Iowa State University.
    [Google Scholar]
  33. Enright, W. H., & Pryce, J. D. (1987). Two FORTRAN packages for assessing initial value methods. ACM Transactions on Mathematical Software (TOMS), 13(1), 1-27.
    [CrossRef]   [Google Scholar]
  34. Alhassan, B., Yusuf, H., & Adamu, A. (2024). Analysis of convergence and stability properties of diagonally implicit 3-point block backward differentiation formula for first order stiff initial value problems. UMYU scientifica, 3(2), 202-210.
    [CrossRef]   [Google Scholar]
  35. Ibrahim, Z. B., & Zawawi, I. S. M. (2019). A Stable Fourth Order Block Backward Differentiation Formulas $\alpha$ for Solving Stiff Initial Value Problems. ASM Science Journal, Special Issue 6, 60-66.
    [Google Scholar]
  36. Aminikhah, H., & Hemmatnezhad, M. (2011). An effective modification of the homotopy perturbation method for stiff systems of ordinary differential equations. Applied Mathematics Letters, 24(9), 1502-1508.
    [CrossRef]   [Google Scholar]
  37. Wu, X. Y. (1998). A sixth-order A-stable explicit one-step method for stiff systems. Computers & Mathematics with Applications, 35(9), 59-64.
    [CrossRef]   [Google Scholar]

Cite This Article
APA Style
Alhassan, B., Saadu, A. S., & Musa, H. (2026). Time Dependent Algorithm using ModifiedThree-PointSuperclassofBlockBackward Differentiation Formula for Solving Stiff Ordinary Differential Equations. ICCK Journal of Applied Mathematics, 2(1), 87–110. https://doi.org/10.62762/JAM.2026.213905
Export Citation
RIS Format
Compatible with EndNote, Zotero, Mendeley, and other reference managers
RIS format data for reference managers
TY  - JOUR
AU  - Alhassan, Buhari
AU  - Saadu, Abubakar Sadiq
AU  - Musa, Hamisu
PY  - 2026
DA  - 2026/03/04
TI  - Time Dependent Algorithm using Modified Three-Point Superclass of Block Backward Differentiation Formula for Solving Stiff Ordinary Differential Equations
JO  - ICCK Journal of Applied Mathematics
T2  - ICCK Journal of Applied Mathematics
JF  - ICCK Journal of Applied Mathematics
VL  - 2
IS  - 1
SP  - 87
EP  - 110
DO  - 10.62762/JAM.2026.213905
UR  - https://www.icck.org/article/abs/JAM.2026.213905
KW  - block backward differentiation formula
KW  - stiff ODEs
KW  - convergence
KW  - consistency
KW  - a-stability
AB  - In this paper, we present a modified three-point superclass of Block Backward Differentiation Formula (BBDF) for the efficient numerical solution of stiff systems of ordinary differential equations (ODEs). The principal enhancement of this work is a structural modification of the classical BBDF that forms a new parameterized superclass of methods, leading to improved stability and reduced error constants compared with the standard three-point BBDFs. The proposed scheme is formulated as a fully implicit block method capable of simultaneously producing three solution approximations within each integration step. A detailed theoretical analysis is conducted to establish the order of accuracy, consistency, zero-stability, and convergence of the method. Stability analysis based on Dahlquist test equation confirms that the scheme is A-stable and suitable for a wide range of stiff ODEs. The nonlinear systems arising from the implicit formulation are efficiently solved using Newton iteration. Numerical experiments on benchmark stiff ODE problems, including systems with rapidly decaying transient components, demonstrate that the modified three-point superclass BBDF achieves higher accuracy and lower computational cost when compared with existing block implicit methods, such as the standard BBDF and diagonally implicit three-point block BDF schemes. Overall, the proposed method provides a robust and computationally efficient alternative for the numerical integration of stiff ODEs, with potential applications in chemical kinetics, control theory, and biological modeling.
SN  - 3068-5656
PB  - Institute of Central Computation and Knowledge
LA  - English
ER  - 
BibTeX Format
Compatible with LaTeX, BibTeX, and other reference managers
BibTeX format data for LaTeX and reference managers
@article{Alhassan2026Time,
  author = {Buhari Alhassan and Abubakar Sadiq Saadu and Hamisu Musa},
  title = {Time Dependent Algorithm using Modified Three-Point Superclass of Block Backward Differentiation Formula for Solving Stiff Ordinary Differential Equations},
  journal = {ICCK Journal of Applied Mathematics},
  year = {2026},
  volume = {2},
  number = {1},
  pages = {87-110},
  doi = {10.62762/JAM.2026.213905},
  url = {https://www.icck.org/article/abs/JAM.2026.213905},
  abstract = {In this paper, we present a modified three-point superclass of Block Backward Differentiation Formula (BBDF) for the efficient numerical solution of stiff systems of ordinary differential equations (ODEs). The principal enhancement of this work is a structural modification of the classical BBDF that forms a new parameterized superclass of methods, leading to improved stability and reduced error constants compared with the standard three-point BBDFs. The proposed scheme is formulated as a fully implicit block method capable of simultaneously producing three solution approximations within each integration step. A detailed theoretical analysis is conducted to establish the order of accuracy, consistency, zero-stability, and convergence of the method. Stability analysis based on Dahlquist test equation confirms that the scheme is A-stable and suitable for a wide range of stiff ODEs. The nonlinear systems arising from the implicit formulation are efficiently solved using Newton iteration. Numerical experiments on benchmark stiff ODE problems, including systems with rapidly decaying transient components, demonstrate that the modified three-point superclass BBDF achieves higher accuracy and lower computational cost when compared with existing block implicit methods, such as the standard BBDF and diagonally implicit three-point block BDF schemes. Overall, the proposed method provides a robust and computationally efficient alternative for the numerical integration of stiff ODEs, with potential applications in chemical kinetics, control theory, and biological modeling.},
  keywords = {block backward differentiation formula, stiff ODEs, convergence, consistency, a-stability},
  issn = {3068-5656},
  publisher = {Institute of Central Computation and Knowledge}
}

Article Metrics
Citations:

Crossref

0

Scopus

0

Web of Science

0
Article Access Statistics:
Views: 29
PDF Downloads: 7

Publisher's Note
ICCK stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and Permissions
CC BY Copyright © 2026 by the Author(s). Published by Institute of Central Computation and Knowledge. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
ICCK Journal of Applied Mathematics

ICCK Journal of Applied Mathematics

ISSN: 3068-5656 (Online)

Email: [email protected]

Portico

Portico

All published articles are preserved here permanently:
https://www.portico.org/publishers/icck/