This paper is concerned with local error estimation in the numerical integration of stiff systems of ordinary differential equations by means of Runge-Kutta methods. With implicit Runge-Kutta methods it is often difficult to embed a local error estimate with the appropriate order and stability properties. In this paper local error estimation based on the information from the last two integration steps (that are supposed to have the same steplength) is proposed. It is shown that this technique, applied to Radau IIA methods, lets us get estimators with proper order and stability properties. Numerical examples showing that the proposed estimate improves the efficiency of the integration codes are presented.

No abstract

No abstract

No abstract

No abstract