  Reaction-diffusion equations are often used in epidemiological models. In
this paper we generalize the algorithm of Meerschaert and Tadjeran for
fractional advection-dispersion flow equations to a coupled system of
fractional reaction-diffusion like equation which arise from vector bourne
disease modeling.
