Create thermalDiffusion.cpp