Second derivative methods with RK stability

Abstract General linear methods are extended to the case in which second derivatives, as well as first derivatives, can be calculated. Methods are constructed of third and fourth order which are A-stable, possess the Runge–Kutta stability property and have a diagonally implicit structure for efficient implementation.