We consider a central difference scheme for the numerical solution of a system of coupled reaction-diffusion equations. We show that the scheme is almost second-order convergent, uniformly in the perturbation parameter, thus improving previous results [2]. We present the results of numerical experiments to confirm our theoretical results.