Solving delay differential equations using intervalwise partitioning by Runge-Kutta method

Embedded singly diagonally implicit Runge-Kutta method is used to solve sti ff systems of delay differential equations. The delay argument is approxima ted using Newton divided difference interpolation. Initially the whole syst em is considered as nonstiff and solved using simple iteration, in the even t that stiffness is indicated, the whole system is considered as stiff and solved using Newton iteration. Numerical results based on one technique to detect stiffness is tabulated and compared with the numerical results when the system is considered as stiff from the beginning. (C) 2001 Elsevier Sci ence Inc. All rights reserved.