Numerical solution of threshold problems in epidemics and population dynamics

Research output: Contribution to journalArticlepeer-review

1 Scopus citations


A new algorithm is proposed for the numerical solution of threshold problems in epidemics and population dynamics. These problems are modeled by the delay-differential equations, where the delay function is unknown and has to be determined from the threshold conditions. The new algorithm is based on embedded pair of continuous Runge-Kutta method of order p=4 and discrete Runge-Kutta method of order q=3 which is used for the estimation of local discretization errors, combined with the bisection method for the resolution of the threshold condition. Error bounds are derived for the algorithm based on continuous one-step methods for the delay-differential equations and arbitrary iteration process for the threshold conditions. Numerical examples are presented which illustrate the effectiveness of this algorithm.

Original languageEnglish (US)
Pages (from-to)40-56
Number of pages17
JournalJournal of Computational and Applied Mathematics
StatePublished - May 1 2015


  • Bisection method
  • Continuous
  • Convergence analysis
  • Delay-differential equations
  • Local error estimation
  • RungeKutta methods
  • Threshold conditions

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Numerical solution of threshold problems in epidemics and population dynamics'. Together they form a unique fingerprint.

Cite this