High order second derivative diagonally Implicit Multistage Integration Methods for ODEs

Construction of second derivative diagonally implicit multistage integration methods (SDIMSIMs) as a subclass of second derivative general linear methods with Runge–Kutta stability property requires to generate the corresponding conditions depending of the parameters of the methods. These conditions which are a system of polynomial equations can not be produced by symbolic manipulation packages for the methods of order p ≥ 5. In this paper, we describe an approach to construct SDIMSIMs with Runge–Kutta stability property by using some variant of the Fourier series method which has been already used for the construction of high order general linear methods. Examples of explicit and implicit SDIMSIMs of order five and six are given which respectively are appropriate for both non-stiff and stiff differential systems in a sequential computing environment. Finally, the efficiency of the constructed methods is verified by providing some numerical experiments.