Ordinary differential equations (ODEs) are among the most important mathe-matical tools used in producing models in the physical sciences, biosciences, chemical sciences, engineering and many more fields. This has motivated re-searchers to provide efficient numerical methods for solving such equations. Most of these types of differential models are stiff, and suitable numerical methods have to be used to simulate the solutions. This paper starts with a sur-vey on the basic properties of stiff differential equations. Thereafter, we present the explicit one-step algorithm proposed by Fatunla to solve stiff systems of first-order scalar ODEs. As an illustrative example, we consider the Robertson problem (RP) which is known to be stiff. The results obtained with the explicit Fatunla method (EFM) are compared with those computed by the solver RADAU which is based on implicit Runge-Kutta methods. Our results are in good agreement with the latter ones.
Cite this paper
Nyengeri, H. , Ndenzako, E. and Nizigiyimana, R. (2019). The Explicit Fatunla’s Method for First-Order Stiff Systems of Scalar Ordinary Differential Equations: Application to Robertson Problem. Open Access Library Journal, 6, e5291. doi: http://dx.doi.org/10.4236/oalib.1105291.
Hairer, E. and Wanner, G. (1996) Solving Ordinary Differential Equations II, Stiff Differential-Algebraic Problems. 2nd Edition, Springer-Verlag, Berlin.
Niemeyer, K.E. and Sung, C.J. (2014) GPU-Based Parallel Integration of Large Numbers of Independent ODE Systems in Numerical Computations with GPUs. Chapter 8, 159-188.
Akinfenwa, O., Jator, S. and Yoa, N. (2011) An Eight Order Backward Differential Formula with Continuous Coefficients of Stiff Ordinary Differential Equations. International Journal of Mathematical and Computer Sciences, 7, 171-176.
Curtiss, C.F. and Hirschfelder, J.O. (1952) Integration of Stiff Equations. Proceedings of the National Academy of Sciences of the United States of America, 38, 235-243. https://doi.org/10.1073/pnas.38.3.235
Kin, J. and Cho, S.Y. (1997) Computational Accuracy and Efficiency of the Time-Splitting Method in Solving Atmospheric Transport/Chemistry Equations. Atmospheric Environment, 31, 2215-2224. https://doi.org/10.1016/S1352-2310(97)88636-0
Shampire, L.F. and Gear, C.W. (1979) A User’s View of Solving Stiff Ordinary Differential Equations. Society for Industrial and Applied Mathematics, 21, 1-17.
Ebady, A.M.N., Habib, H.M. and El-Zahar, E.R. (2012) A Fourth Order A-Stable Explicit One-Step Method for Solving Stiff Differential Systems Arising in Chemical Reactions. International Journal of Pure and Applied Mathematics, 81, 803-812.
Gupta, G.K., Sacks-Davis, R. and Tischer, P.E. (1985) A Review of Recent Developments in Solving ODEs. Computing Surveys, 17, 5-47. https://doi.org/10.1145/4078.4079
Bjurel, G., Dahlquist, G.G., Lindberg, B., Linde, S. and Oden, L. (1970) Survey of Stiff Ordinary Differential Equations. Report NA 70.11, KTH Royal Institute of Technology, Stockholm, Sweden.
Robertson, H.H. (1967) The Solution of a Set of Reaction Rate Equations. In: Walsh, J., Ed., Numerical Analysis: An introduction, Academic Press, Cambridge, Massachusetts, 178-182.
Little, W., Fowler, H.W., Coulson, J. and Onions, C.T. (1068) The Shorter Oxford English Dictionary on Historical Principles. 3rd Edition, Clarendon Press, Oxford England.
Hairer, E. and Wanner, G. (1999) Stiff Differential Equations Solved by Randau Methods. Journal of Computational and Applied Mathematics, 111, 93-111. https://doi.org/10.1016/S0377-0427(99)00134-X
Anake, T.A., Bishop, S.A., Adesanya, A.O. and Agarana, M.C. (2014) An -Stable Method for Solving Initial-Value Problems of Ordinary Differential Equations. Advances in Differential Equations and Control Processes, 13, 21-35.
Cash, J.R. (2003) Review Paper: Efficient Numerical Methods for the Solution of Stiff Initial-Value problems and Differential Algebraic Equations. Proceedings of the Royal Society A, 459, 797-815. https://doi.org/10.1098/rspa.2003.1130
Sekar, S. and Nalini, M. (2015) Numerical Solution of Linear and Nonlinear Stiff Problems Using Adomian Decomposition Method. IOSR Journal of Mathematics, 11, 14-20.
Gaffney, P.W. (1984) A Performance Evaluation of Some FORTRAN Subroutines for the Solution of Stiff Oscillatory ODEs. ACM Transactions on Mathematical Softwares, 10, 58-72. https://doi.org/10.1145/356068.356073
S?derlind, G., Jay, L. and Calvo, M. (2015) Stiffness 1952-2012: Sixty Years in Search of a Definition. BIT Numerical Mathematics, 55, 531-558. https://doi.org/10.1007/s10543-014-0503-3
S?derlind, G. (2006) The Logarithmic Norm. His-tory and Modern Theory. BIT Numerical Mathematics, 46, 631-652. https://doi.org/10.1007/s10543-006-0069-9
Madro?ero, J. and Piraux, B. (2009) Explicit Time-Propagation Method to Treat the Dynamics of Driven Complex Systems. Physical Review A, 80, Article ID: 033409.