A novel numerical approach for solving delay differential equations arising in population dynamics

In this paper, the initial-value problem for a class of first order delay differential equations, which emerges as a model for population dynamics, is considered. To solve this problem numerically, using the finite difference method including interpolating quadrature rules with the basis functions, we construct a fitted difference scheme on a uniform mesh. Although this scheme has the same rate of convergence, it has more efficiency and accuracy compared to the classical Euler scheme. The different models, Nicolson's blowfly and Mackey-Glass models, in population dynamics are solved by using the proposed method and the classical Euler method. The numerical results obtained from here show that the proposed method is reliable, efficient, and accurate.