Abstract
This article analyses a mode of two-dimensional fox-rabies dynamics with advection terms. The fox population is divided into two species: susceptible S and infective I. Finite-difference techniques are employed to compute the numerical solutions of this model. Iterative schemes with second-order accuracy are derived and their convergence is investigated by using the maximum principle analysis. An illustrative numerical example with graphical reis given to highlight the efficiency of the derived numerical schemes.