Document Type


Publication Date



We study numerical solution of systems of delay-differential equations in which the delay function, which depends on the unknown solution, is defined implicitly by the threshold condition. We study discrete variable numerical methods for these problems and present error analysis. The global error is composed of the error of solving the differential systems, the error from the threshold conditions and the errors in delay arguments. Our theoretical analysis is confirmed by numerical experiments on threshold problems from the theory of epidemics and from population dynamics.

Copyright Statement

This is an author-produced, peer-reviewed version of this article. © 2009, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License ( The final, definitive version of this document can be found online at Journal of Computational and Applied Mathematics, doi: 10.1016/

Included in

Mathematics Commons