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.