Laminar boundary layer slip flow from a stretching surface in a nanofluidsaturated homogenous, isotropic porous medium is studied numerically. A Newtonian heating boundary condition in the presence of thermal radiation is incorporated and a Darcy model utilized for the porous medium. The model used for the nanofluids include the effects of Brownian motion and thermophoresis. A group theoretical analysis is conducted to generate similarity transformations. The governing transport equations are nondimensionalized and rendered into a set of coupled similarity ordinary differential equations using similarity transformations. The transformed equations are then solved using the Runge-Kutta-Fehlberg fourth-fifth order numerical method with shooting technique. It is shown that the physical quantities of interest depend on a number of parameters. The results are presented in tabular and graphical forms. Comparison of the present numerical solutions with published work shows very good agreement. The study finds applications in high-temperature nanotechnological materials processing. C⃝ 2014 Wiley Periodicals, Inc. Heat Trans Asian Res, 44(8): 681-695, 2015; Published online in Wiley Online Library (wileyonlinelibrary.com/journal/htj).