The main aim of the current study is to determine the effects of the temperature dependent viscosity and thermal conductivity on magnetohydrodynamics (MHD) flow of a non-Newtonian fluid over a nonlinear stretching sheet. The viscosity of the fluid depends on stratifications. Moreover, Powell–Eyring fluid is electrically conducted subject to a non-uniform applied magnetic field. Assume a small magnetic reynolds number and boundary layer approximation are applied in the mathematical formulation. Zero nano-particles mass flux condition to the sheet is considered. The governing model is transformed into the system of nonlinear Ordinary Differential Equation (ODE) system by using suitable transformations so-called similarity transformation. In order to calculate the solution of the problem, we use the higher order convergence method, so-called shooting method followed by Runge-Kutta Fehlberg (RK45) method. The impacts of different physical parameters on velocity, temperature and concentration profiles are analyzed and discussed. The parameters of engineering interest, i.e., skin fraction, Nusselt and Sherwood numbers are studied numerically as well. We concluded that the velocity profile decreases by increasing the values of S t , H and M. Also, we have analyzed the variation of temperature and concentration profiles for different physical parameters.