Implementation of Diagonally Implicit Multistage Integration Methods for Ordinary Differential Equations

Abstract

We investigate the implementation of diagonally implicit multistage integration methods (DIMSIMs). The implementation issues addressed are the local error estimation, changing stepsize using the Nordsieck technique, and the construction of continuous interpolants. Numerical experiments with a method of order three indicate that the error estimates that have been constructed are very reliable in both a fixed and a variable stepsize environment.