Simulation of nonlinear reaction-diffusion equations