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 embedded in a porous medium is considered. Karniadakis and Beskok’s quadratic order slip boundary conditions are taken into account. A closed-form analytical solution of the momentum equation is used to derive the analytical solution of the heat transfer equation in terms of confluent hypergeometric functions with quadratic order thermal slip boundary condition. Accuracy of the present results is assured with the numerical solution obtained by the 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 higher estimations of tangential momentum and energy accommodation coefficients.