The numerical solution of the nonlinear system of stiff differential equations by the modified matrix-exponential method

Document Type : Research Paper

Authors

Department of Mathematics, Science of Mathematics Faculty, University of Mazandaran, Babolsar, Iran

Abstract

In this paper, the modified matrix exponential (MME) method under zero-order hold (ZOH) assumption, is applied to solve systems of stiff ordinary differential equations. Some examples are given to illustrate the accuracy and effectiveness of the method. We compare our results with results obtained by matrix exponential (ME) method and by the Matlab ode23 solver.

Keywords

[1] A. Abdi, S.A. Hosseini and H. Podhaisky, Adaptive linear barycentric rational finite differences method for stiff ODEs, J. Comput. Appl. Math. 357 (2019), 204–214.
[2] F. Allffi-Pentini, V.D. Fonzo and V. Paris, A novel algorithm for the numerical integration of systems of ordinary differential equations arising in chemical problem, J. Math. Chem. 33 (2003), 1–15.
[3] M. Calvo, S. Gonzalez-Pinto and J.I. Montijano, Runge–Kutta methods for the numerical solution of stiff semilinear system, BIT Numer. Math. 40 (2000), 611–639.
[4] M.S.H. Chowdhury, I. Hashim and Md. Alal Hosen, Solving linear and non-linear stiff system of ordinary differential equations by multistage Adomian decomposition method, Proc. The Third Intl. Conf. Adv. Appl. Sci. Environ. Technol. 125 (2015), 79–82.
[5] M. Fattahi and M. Matinfar, The numerical solution of the second kind of Abel equations by the modified matrix-exponential method, Int. J. Nonlinear Anal. Appl. In Press, (2023), 1–7.
[6] G.Y. Kulikov, R. Weiner and G.M. Amiraliyev, Doubly quasi-consistent fixed-step size numerical integration of stiff ordinary differential equations with implicit two-step peer methods, J. Comput. Appl. Math. 340 (2018), 256–275.
[7] U. Lepik, Haar wavelet method for solving stiff differential equations, Math. Model. Anal. 14 (2009), 79–82.
[8] S. Li, Canonical Euler splitting method for nonlinear composite stiff evolution equations, Appl. Math. Comput. 289 (2016), 220–236.
[9] Y. Park and K.T. Chong, The numerical solution of the point kinetics equation using the matrix exponential method, Ann. Nuclear Energy 55 (2013), 42–48.
Volume 15, Issue 4
April 2024
Pages 339-347
  • Receive Date: 24 February 2023
  • Revise Date: 31 May 2023
  • Accept Date: 11 June 2023