Here is a mathcad code I wrote on calculating the temperature field of a 3D multilayer structure.
First the analytical solution of the 3D structure subjected to impulse excitation is solved in Laplace Domain and the then the temperature field subjected to arbitrary time dependent power input is obtained using convolution theory.