Two-dimensional MHD mixed convection boundary layer flow of heat and mass transfer stagnation-point flow of a non-Newtonian power-law nanofluid towards a stretching surface in the presence of thermal radiation and heat source/sink is investigated numerically. The non-Newtonian nanofluid model incorporates the effects of Brownian motion and thermophoresis. The basic transport equations are made dimensionless first and the complete nonlinear differential equations with associated boundary conditions are solved numerically by finite element method (FEM). The numerical calculations for velocity, temperature, and nanoparticles volume fraction profiles for different values of the physical parameters to display the interesting aspects of the solutions are presented graphically and discussed. The skin friction coefficient, the local Nusslet number and the Sherwood number are exhibited and examined. Our results are compatible with the existing results for a special case.