In this article, we study the effect of thermal radiation and heat generation/absorption on hydromagnetic threedimensional free convention flow of nanofluid over a vertical stretching surface. Using a set of self similar transformations the model equations are solved numerically by applying a fifth-order Runge-Kutta-Fehlberg scheme with a shooting technique using appropriate boundary conditions. The impact of non-dimensional parameters of the x- (φ’ and M), y- (N), z- (θ), velocity components, pressure (G), temperature (φ) and mass (H) are presented and discussed graphically. The results compare favorably with previously published work. Results indicate that an increase in Brownian motion and thermophoresis parameters increases the temperature and concentration profiles. © 2017 by American Scientific Publishers All rights reserved.