In this paper we study the partial slip effects on the flow and heat transfer of an incompressible non-Newtonian fluid over a nonlinear stretching sheet. The velocity slip boundary condition based on the Sisko constitutive fluid model is introduced. Suitable dimensionless variables are used to convert the governing partial differential equations into ordinary differential equations. Numerical solutions of these equations are obtained by the Runge–Kutta Fehlberg method. Additionally, the exact analytical solutions are presented in some special cases. The computational results for the velocity, temperature, skin-friction coefficient, and Nusselt number are presented in graphical and tabular forms. To validate the numerical results obtained, a comparison is made with the exact analytical solutions. The analysis of the results obtained shows that enhancement in the velocity slip parameter reduces the velocity as well as the momentum boundary layer thickness. However, quite the opposite is true with the temperature and corresponding thermal boundary layer thickness.