Stability Analysis of Diagonally Implicit Two Derivative Runge-Kutta methods for Solving Delay Differential Equations

The stability properties of fourth and fifth-order Diagonally Implicit Two Derivative Runge-Kutta method (DITDRK) combined with Lagrange interpolation when applied to the linear Delay Differential Equations (DDEs) are investigated. This type of stability is known as P-stability and Q-stability. Their stability regions for (lambda, mu is an element of R) and (mu is an element of C, lambda = 0) are determined. The superiority of the DITDRK methods over other same order existing Diagonally Implicit RungeKutta (DIRK) methods when solving DDEs problems are clearly demonstrated by plotting the efficiency curves of the log of both maximum errors versus function evaluations and the CPU time taken to do the integration.