Second derivative Runge-Kutta collocation methods based on Lobatto nodes for stiff systems

Second derivative Runge-Kutta collocation methods for the numerical solution of sti system of rst order initial value problems in ordinary dierential equations are derived and studied. The inclusion of the second derivative terms enabled us to derive a set of methods which are convergent with large regions of absolute stability. Although the implementation of the methods remains iterative in a precisely dened way, the advantage gained makes them suitable for solving sti system of equations with large Lipschitz constants. The derived methods are illustrated by the applications to some test problems of sti system found in the literature and the numerical results obtained conrm the potential of the second derivative methods.