∂c(r,t)/∂t=L(c(r,t) )+D∇^2 c(r,t)
how to discretize the time derivative?thanks