Abstract
A second-order, unconditionally-stable, finite-difference scheme is developed for the numerical solution of the SI model of fox-rabies dynamics. The local stability of the scheme, by direct inspection of the eigenvalues dependent on the time step size and on two parameters, is shown to be unconditionally stable.