Numerical analysis of heat flow