On the implementation of ESIRK methods for stiff IVPs

The effective order singly-implicit methods (ESIRK) are designed for solving stiff IVPs. These generalizations of SIRK methods are shown to have some computational advantages over the classical SIRK methods by moving the abscissae inside the integration interval [6]. In this paper, we consider some of the important computational aspects associated with these methods. We show that the ESIRK methods can be implemented efficiently by the comparsion with the standard stiff solvers RADAU5 and LSODE.