Solving Stiff Systems of ODEs by Explicit Methods with Conformed Stability Domains

Abstract

The Cauchy problem for a stiff system of ODEs is considered. The explicit m-stage first order methods of the Runge-Kutta type are designed with stability domains of intermediate numerical schemes conformed with the stability domain of the basic scheme. Inequalities for accuracy and stability control are obtained. A numerical algorithm based on the first-order method and the five-stage fourth order Merson method is developed. The algorithm is aimed at solving large-scale systems of ODEs of moderate stiffness with low accuracy. It has been included in the library of solvers of the ISMA simulation environment. Numerical results showing growth of the efficiency are given.