A New Class of Third Derivative Fourth-Step Exponential Fitting Numerical Integrator for Stiff Differential Equations
Abstract
In this paper, we construct a new class of four-step third derivative exponential fitting integrator of order six for the numerical integration of stiff initial-value problems of the type: y<sup>Phys.Rev. ime = f(x,y), \:\: y(x<sub>0</sub></sup> ) = y<sub>0</sub>. The implicit method has free parameters which allow it to be fitted automatically to exponential functions. For the purpose of effective implementation of the new proposed method, we adopted the techniques of splitting the method into predictor and corrector schemes. The numerical analysis of the stability of the new method was discussed; the results show that the new method is A-stable. Finally, numerical examples are presented, to show the efficiency and accuracy of the new method
References
G. Dahlquist, A special Stability Problem for Linear Multi-step Method, BIT 3, 27-43, (1963).
J.D. Lambert, Numerical Method in Ordinary Differential Equations with Initial Value Problem, John Wiley, (1973).
W. H. Enright, Second Derivative Multi-step Methods for Stiff Ordinary Differential Equations, SIAM J. Numerical Analysis II, 321-331, (1974).
D. Voss, A fifth order exponentially fitted formula, SIAM J.Numer Anal., 25 (3), 670-678, (1988).
C.E. Abhulimen and G.E. Omeike, A sixth-order exponentially fitted scheme for the numerical solution of systems of ordinary differential equations, . Journal of Applied Mathematics and Bioinformatics, 1 (1), 175-186, (2011) .
F.O. Otunta and C.E. Abhulimen, A Sixth Order Multi Derivative Multi-step Method for stiff System of ODE, International Journal of Numerical.Maths. (ITNM) 2 (1), 248-268, (2006).
C.E. Abhulimen, Exponential Fitting Predictor-Corrector formula for stiff systems of Ordinary Differential Equations, .International Journal of Computational and Applied Mathematics 4 (2), 115-126, (2009).
F.O. Otunta and C.E. Abhulimen, (2005) A 4th Order Exponentially Fitted Multiderivative method for Stiff IVPs, Nigerian Association of Mathematical Physics, 9, 295-306, (2005) .
W. Liniger and R.A. Willioughby, Efficient Integration methods for stiff systems of ordinary differential equations, SIAM J. NumerAnnal, 7, 47-65, (1970).
T.L. Brown, Some Characteristics of implicit Multi-step Multi-derivative Integration formulas, SIAM J.Numerical Math. 34 (1), 59-84, (1977) .
J.R. Cash, On the Exponential Fitting of Composite Multi-derivative Linear Mult-istep Methods, SIAM J. Numerical Annal. 18 (5), 808-821, (1981).
W. H. Enright and J.O. Pryce, Two FORTRAN Packages for Assessing IV Methods, Technical Report 16/83, Department of Computer Science, University of Toronto, Canada, (1983).
S.A. Okunuga, Fourth order composite two step method for stiff problems, Int. J. Comput. Math., 2, 39-47, (1997).
This is an Open Access article distributed under the terms of the Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, adaptation, and reproduction in any medium, provided that the original work is properly cited.
