A numerical study has been carried out to analyze the constant heat flux, internal heat generation, variable viscosity and thermal radiation effects on the flow and heat transfer of a Newtonian fluid over an exponentially stretching porous sheet. Using a similarity transformation, the governing partial differential equations are transformed into coupled, non-linear ordinary differential equations with variable coefficients. Numerical solutions to these equations subject to appropriate boundary conditions are obtained by using an efficient Chebyshev spectral method. The effects of various physical parameters such as viscosity parameter, the suction parameter, the radiation parameter, internal heat generation or absorption parameter and the Prandtl number on velocity and temperature are discussed by using graphical approach. Moreover, numerical results indicate that in the presence of constant heat flux, the skin-friction coefficient as well as Nusselt number is strongly affected by the viscosity parameter, suction parameter, radiation parameter and the internal heat generation parameter.