This paper deals with the numerical investigation of unsteady flow in a porous medium towards a stretching sheet. Twowater based nanofluids namely and are taken into consideration. With the auxiliary conditions the governing partial differentialequations are converted to ordinary differential equations with the pertinent corresponding conditions. MATLAB function “bvp4c” is applied to solve the resulting governing equations together with the boundary conditions. The effect of various parameters namely: solid volume fractions, Prandtl number, porosity parameter, thermal radiation and the unsteadiness parameter have been discussed. The numerical results obtained for different parameters are presented through plots. Comparison with known results is presented and it is found excellent.