This paper is devoted to introduce a numerical simulation using the differential transformation method (DTM) with a theoretical study for the effect of viscous dissipation on the steady flow with heat transfer of Newtonian fluid towards a permeable stretching surface embedded in a porous medium with a second order slip. The governing nonlinear partial differential equations are converted into a system of nonlinear ordinary differential equations (ODEs) by using similarity variables. The resulting ODEs are successfully solved numerically with the help of DTM. Graphic results are shown for nondimensional velocities and temperatures. The effects of the porous parameter, the suction (injection) parameter, Eckert number, first and second order velocity slip parameters and the Prandtl number on the flow and temperature profiles are given. Moreover, the local skin-friction and Nusselt numbers are presented. Comparison of numerical results is made with the earlier published results under limiting cases.