In this communication, an analytical solution for the thermal transfer of Newtonian fluid flow with quadratic order thermal and velocity slips is presented for the first time. The flow of a Newtonian fluid over a stretching sheet which is embedded in a porous medium is considered. Karniadakis and Beskok’s quadratic order slip boundary conditions are taking into account. A closed form of analytical solution of momentum equation is used to derive the analytical solution of heat transfer equation in terms of confluent hyper-geometric function with quadratic order thermal slip boundary condition. Accuracy of present results is assured with the numerical solution obtained by Iterative Power Series method with shooting technique. The impacts of porous medium parameter, tangential momentum accommodation coefficient, energy accommodation coefficient on velocity and temperature profiles, skin friction coefficient and reduced Nusselt number are discussed. The Nusselt number increases with the higher estimations of tangential momentum and energy accommodation coefficients.