diffusion equation using finite element

Hi,

I would like to write a simple subroutine to solve diffusion equation d(ro)/dt = k∇2(ro). using finite element method.

Any suggestions ?

Thanks,

Ayoub