The novel heat conduction equation which takes into
account both the heat flux relaxation (the finite speed of
thermal wave propagation) and the relaxation of internal heat
source capacity is solved numerically by the MacCormack
predictor - corrector method. The temperature field in the
body due to a step change in heat flux at the surface is
analysed for various ratios of relaxation times. It is shown
that the relaxation of internal heat source capacity can have
an appreciable influence on the temperature distribution in the
body.