Simulation of reaction-diffusion equations