The source is here.
Solves the 3D time-dependent heat equation
by direct convolution with the Green's function,