The principal concern of the present work is to investigate the effects of thermally developed Brownian motion and thermophoresis diffusion in non-Newtonian nanofluid through an inclined stretching surface with effect of chemical reaction and thermal radiation. By using compatible similarity transformations the governing equations of momentum, energy and concentration profiles are converted to a set of non-linear differential equations. A Mathematica package BVPh 2.0 is applied to examine the given system of equations. Influences of developing parameters such as magnetic parameter, radiation, Lewis number, Prandtl number, mixed convective parameter, Brownian motion and thermophoresis parameter on velocity, temperature and nanoparticles concentration are displayed through graphical experiment. A comparative study between OHAM and previously published results is made. It is noticed that OHAM can overcome the earlier restriction, assumptions and limitation of classical perturbation schemes. The implementation of OHAM is reasonably simple through Mathematica package need to define the governing equations, boundary conditions; suitable auxiliary linear operator and initial guess etc. The key advantage of the suggested technique is that it can be used direct way in highly nonlinear differential equations without using discretization, linearization and round-off errors. (C) 2017 Elsevier Ltd. All rights reserved.