Abstract
We describe a numerical algorithm for the simulation of the spread of rabies in a spatially distributed fox population. The model considers both territorial and wandering rabid foxes and includes a latent period for the infection. The resulting systems are mixtures of partial differential and integral equations. They are discretized in the space variable by central differences of second order and by the composite trapezoidal rule. In a second step, they are discretized in time by explicit continuous Runge–Kutta methods of fourth order for ordinary and delay differential systems. The results of the numerical calculations are compared for latent periods of fixed and exponentially distributed length and for various proportions of territorial and wandering rabid foxes. The speeds of spread observed in the simulations are compared to spreading speeds obtained by analytic methods and to observed speeds of epizootic frontlines in the European rabies outbreak 1940 to 1980.
Original language | English (US) |
---|---|
Journal | Mathematics and Computers in Simulation |
DOIs | |
State | Accepted/In press - Jan 1 2018 |
Keywords
- Continuous Runge–Kutta method
- Diffusing versus territorial rabid foxes
- Latent period
- Method of lines
- Spreading speed
ASJC Scopus subject areas
- Theoretical Computer Science
- General Computer Science
- Numerical Analysis
- Modeling and Simulation
- Applied Mathematics