second derivative general linear methods (sglms) as an extension of general linear methods (glms) have been introduced to improve the stability and accuracy properties of glms. the coefficients of sglms are given by six matrices, instead of four matrices for glms, which are obtained by solving nonlinear systems of order and usually runge--kutta stability conditions. in this p...