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 an equation that arises from vector bourne disease modeling.