Numerical solutions for laminar heat transfer of a non-Newtonian fluid in the thermal entrance region of a square duct are presented for three thermal boundary conditions. The power-law model characterises the non-Newtonian behavior. The numerical results show that for each flow behavior index the Nusselt number decreases from a maximum value at the entry plane to a limiting value when both velocity and temperature profiles are fully developed. The results are compared with the available solutions for Newtonian fluid and excellent agreement is found. © 1977.