In this paper the viscoelastic flow and heat transfer over a non-linearly stretching sheet with the power law velocity of the form uw= cxn is investigated for the first time. A prescribed power-law surface temperature distribution of the form Tw= T ¥ + Axn is considered. Mathematical model is constructed through the constitutive equations of second grade fluid. The arising non-linear boundary value problem has been treated analytically by a powerful optimal homotopy analysis method (OHAM). The solutions are found in excellent agreement with the obtained numerical solutions in the case of Newtonian fluid. The results show that velocity and skin friction coefficient have direct relationship with the power-law index n. Further the thermal boundary layer becomes thinner when larger values of n are taken into account.