Abstract
Diagonally implicit multistage integration methods (DIMSIMs) are widely utilized in finding the solution to any problems in the subject of ordinary differential equations. These methods are selected from the general linear methods, which is considerable potential for efficient implementations. The extrapolation is derived from the stability of the explicit Runge-Kutta methods. In this paper, the combination of DIMSIMs with Richardson extrapolation of different orders shows that numerical solutions give higher accuracy when the extrapolation is applied with the base method.