Abstract
A one-dimensional model of fox-rabies of two nonlinear partial differential equations of hyperbolic type is studied. Finite difference techniques are applied to compute the numerical solutions of the initial/boundary value problem. The convergence of the resulting schemes, which have a second order accuracy in space and time, is investigated. The method is tested for different values of advection rate; numerical and graphical results showed that the method is consistent with the dynamic behaviour of fox-rabies.