Abstract
This study presents the numerical solutions of boundary layer flow and heat transfer over a stretching sheet with viscous dissipation and internal heat generation. Thermal boundary condition on the surface, namely prescribed heat flux (PHF) is used. The governing nonlinear partial differential equations are transformed into nonlinear ordinary differential equations by applying the similarity transformations before reduced to the system of first order ordinary differential equations. Then the system of first order ordinary differential equations is solved numerically using an implicit finite difference scheme, known as the Keller-box method. The numerical solutions are generated using MATLAB. Temperature profiles and the temperature gradient for some values of the Prandtl number, Eckert number and heat/ source sink parameter are presented in figures and discussed in details.