In this work we expose, first, a method of successive approximations for the solution of the system of equations of the thermal boundary layer of the laminar flow in non stationary regime in three dimensions. Then, we find analytical solutions of the differential equations resulting from this method of successive approximations for any Prandtl numbers. Finally, these analytical solutions are also computed numerically and presented graphically for certain values of the Prandtl number.