This paper subdivided the interior solid into triangular shape of equal size to calculate the temperature distribution around the square heat source of it, and compared calculated values with measured ones. The result obtained are as follows. 1) It was found that we can calculate the temperature distribution around the square heat source of interior solid by the variational method of the finite element method as the calculated values were almost accord with the measured ones. 2) The temperature distributed were higher when the distances between heat source were farther and lower when those nearer. 3) Vertical surface temperature distribution is remarkably efficient by thermal conductivities.