The flow and heat transfer characteristics of incompressible viscous flow over a nonlinearly stretching sheet with the presence of viscous dissipation is investigated numerically. The similarity transformation reduces the time‐independent boundary layer equations for momentum and thermal energy into a set of coupled ordinary differential equations. The obtained equations, including nonlinear equation for the velocity field f and differential equation by variable coefficient for the temperature field θ, are solved numerically by using the fourth order of Runge‐Kutta integration scheme accompanied by shooting technique with Newton‐Raphson iteration method. The effect of various values of Prandtl number, Eckert number and nonlinear stretching parameter are studied. The results presented graphically show some behaviors such as decrease in dimensionless temperature θ due to increase in Pr number, and curve relocations are observed when heat dissipation is considered.