Numerical solution of threshold problems in epidemics and population dynamics

• Z. Bartoszewski ^[1] ; Z. Jackiewicz ^[2] ; Y. Kuang ^[2]
  1. [1] Gdańsk University of Technology

     Gdańsk University of Technology

     Gdańsk, Polonia

     
  2. [2] Arizona State University

     Arizona State University

     Estados Unidos

     
• Localización: Journal of computational and applied mathematics, ISSN 0377-0427, Vol. 279, Nº 1 ((1 May 2015)), 2015, págs. 40-56
• Idioma: inglés
• DOI: 10.1016/j.cam.2014.10.020
• Texto completo no disponible (Saber más ...)
• Resumen

  • 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.