Sin-Cos-Taylor-Like method for solving stiff ordinary diffrential equations

Authors

  • Rokiah @ Rozita Ahmad
  • Nazeeruddin Yaacob

DOI:

https://doi.org/10.11113/mjfas.v1n1.13

Keywords:

Stiff ordinary differential equations, Explicit methods, A-stable, L-stable,

Abstract

This paper discusses the derivation of an explicit Sin-Cos-Taylor-Like method for solving stiff ordinary differential equations, which is a formulation of the combination of a polynomial and the exponential function. This new method requires extra work to evaluate a number of differentiations of the function involved. However, the result shows smaller errors when compared to the results from the explicit classical fourth-order Runge-Kutta (RK4) and the Adam-Bashforth-Moulton (ABM) methods. Implicit methods could work well for stiff problems but have certain drawbacks especially when discussing about the cost. Although extra work is required, this explicit method has its own advantages. Besides providing excellent results, the cost of computation using this explicit method is much cheaper than the implicit methods. We also considered the stability property for this method since the stability property of the classical explicit fourth order Runge-Kutta method is not adequate for the solution of stiff problems. As a result, we find that this explicit method is of order-6, which has been developed, and proved to be both A-stable and L-stable.

References

G. K. Gupta, R. Sacks-Davis, P. E. Tescher, ACM Computing Surveys, 17 (1985) 5-47.

R. C. Aiken, Stiff Computation, Oxford University Press, Oxford, 1985.

Butcher, Journal of Computational and Applied Mathematics, 125 (2000) 1 – 29.

J. D. Lambert, Numerical methods for Ordinary Differential Equations, the Initial Value Problem. Chichester: John Wiley & Sons, 1991.

X. Y. Wu, Computers Math. Applications, 35 (1998) 59 – 64.

R. Ahmad, N. Yaacob, Explicit Taylor-Like Method for Solving Stiff Differential Equations, Technical Report. LT/M Bil.3/2004.

R. Ahmad, N. Yaacob, Explicit Sine-Taylor-Like Method for Solving Stiff Differential Equations, Technical Report. LT/M Bil.8/2004.

S. Wolfram, Mathematica: A System For Doing Mathematics By Computer, 2nd ed. Addison-Wesley Publishing Company, 1991

J. D. Lambert, Computational Methods in Ordinary Differential Equations. John Wiley & Sons Ltd, London, 1973.

Downloads

Published

16-06-2014