Numerical solution of a nonlinear reaction diffusion equation