Numerical Solution of Diffusion Equations