Abstract
We have recently developed a generic approach for solving neutral delay differential equations based on the use of a continuous Runge–Kutta formula with defect control and investigated its convergence properties. In this paper, we describe a method, DDVERK, which implements this approach and justify the strategies and heuristics that have been adopted. In particular we show how the assumptions related to error control, stepsize control, and discontinuity detection (required for convergence) can be efficiently realized for a particular sixth-order numerical method. Summaries of extensive testing are also reported.
Similar content being viewed by others
References
S.C. Corwin and S. Thompson, DKLAG6: Solution of systems of functional differential equations with state dependent delays, Technical Report 96–002, Computer Science Department, Radford University, Radford (1996).
J.R. Dormand and P.J. Prince, Runge-Kutta triples, Comput. Math. Appl. 12A(9) (1986) 1007–1017.
W.H. Enright, Analysis of error control strategies for continuous Runge-Kutta methods, SIAM J. Numer. Anal. 26 (1989) 588–599.
W.H. Enright, A new error-control for initial value solvers, Appl. Math. Comput. 31 (1989) 288–301.
W.H. Enright, The relative efficiency of alternative defect control schemes for high-order continuous Runge-Kutta formulas, SIAM J. Numer. Anal. 30(5) (1993) 1419–1445.
W.H. Enright and H. Hayashi, Convergence analysis of the solution of retarded and neutral delay differential equations by continuous numerical methods, SIAM J. Numer. Anal. 35(2) (1998).
W.H. Enright and H. Hayashi, The evaluation of numerical software for delay differential equations, in: The Quality of Numerical Software: Assessment and Enhancement, ed. R.F. Boisvert (Chapman and Hall, London, 1997) pp. 179–192.
W.H. Enright, K.R. Jackson, S.P. Nørsett and P.G. Thomsen, Interpolants for Runge-Kutta formulas, ACM Trans. Math. Software 12 (1986) 193–218.
W.H. Enright, K.R. Jackson, S.P. Nørsett and P.G. Thomsen, Effective solution of discontinuous IVPs using a Runge-Kutta formula pair with interpolants, Appl. Math. Comput. 27 (1988) 313–335.
I. Gladwell, L.F. Shampine, L.S. Baca and R.W. Brankin, Practical aspects of interpolation in use Runge-Kutta codes, SIAM J. Sci. Statist. Comput. 8 (1987) 322–341.
H. Hayashi, Numerical solution of retarded and neutral delay differential equations using continuous Runge-Kutta methods, Ph.D. thesis, Department of Computer Science, University of Toronto, Toronto, Canada (1996).
M.K. Horn, Fourth-and fifth-order, scaled Runge-Kutta algorithms for treating dense output, SIAM J. Numer. Anal. 20(3) (1983) 558–568.
T.E. Hull, W.H. Enright and K.R. Jackson, User's guide for DVERK - a subroutine for solving non-stiff ODE's, Technical Report 100, Department of Computer Science, University of Toronto (1976).
C.A.H. Paul, A user-guide to ARCHI: An explicit Runge-Kutta code for solving delay and neutral differential equations and parameter estimation problems, Technical Report 283 (extended), Department of Mathematics, University of Manchester, Manchester (1997).
L.F. Shampine, Interpolation for Runge-Kutta methods, SIAM J. Numer. Anal. 22 (1985) 1014–1027.
L.F. Shampine, Some practical Runge-Kutta formulas, Math. Comp. 46 (1986) 135–150.
J.H. Verner, Differentiable interpolants for high-order Runge-Kutta methods, SIAM J. Numer. Anal. 30(5) (1993) 1446–1466.
T.E. Wheldon, J. Kirk and H.M. Finlay, Cyclical granulopoiesis in chronic granulocytic leukemia: A simulation study, Blood 43 (1974) 379–387.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Enright, W., Hayashi, H. A delay differential equation solver based on a continuous Runge–Kutta method with defect control. Numerical Algorithms 16, 349–364 (1997). https://doi.org/10.1023/A:1019107718128
Issue Date:
DOI: https://doi.org/10.1023/A:1019107718128