The aim of this work is to obtain the numerical solutions for the boundary layer flow of heat transfer of incompressible
viscous fluid. The governing partial differential equations are converted into ordinary differential equation by using a similarity
transformation. The nonlinear equation governing the flow problem is modeled and then solved numerically by means of a
successive linearization method (SLM). The numerical results are derived in tables for comparisons. The important result of this
comparison is to show the high precision of the SLM in solving system of nonlinear differential equations. Graphical outcomes of
various parameters such as Prandtl number (Pr) and Eckert number (Ec) on the flow, field are discussed and analyzed. Besides
this the present results have been tested and compared with the available published results in a limiting manner and an excellent
agreement is found.