The problem of the flow and heat transfer over an unsteady stretching sheet embedded in a porous medium in the presence of thermal radiation is studied theoretically and numerically. The continuity, momentum, and energy equations, which are coupled nonlinear partial differential equations, are reduced to a set of two nonlinear ordinary differential equations. Special attention is given to study the convergence of the proposed method. The error estimation is also given. The effects of various parameters, such as the Darcy parameter, the radiation parameter, and the Prandtl number, on the flow and temperature profiles, as well as on the local skin-friction coefficient and the local Nusselt number are presented and discussed. The results obtained agree very well with the data obtained by the Runge–Kutta method coupled with the shooting technique.