Second Derivative Block Hybrid Methods for the Numerical Integration of Differential Systems

Abstract

The second derivative block hybrid method for the continuous integration of differential systems within the interval of integration was derived. The second derivative block hybrid method maintained the stability properties of the Runge–Kutta methods suitable for solving stiff differential systems. The lack of such stability properties makes the continuous solution not reliable, especially in solving large stiff differential systems. We derive these methods by using one intermediate off-grid point in between the familiar grid points for continuous solution within the interval of integration. The new family had a high accuracy, non-overlapping piecewise continuous solution with very low error constants and converged under the suitable conditions of stability and consistency. The results of computational experiments are presented to demonstrate the efficiency and usefulness of the methods, which also indicate that the block hybrid methods are competitive with some strong stability stiff integrators.