The study examined a three-dimensional unsteady Magnetohydrodynamic non-Newtonian nanofluid flow with magnetic induction, Lorentz force, viscous dissipation and thermophoresis between two parallel horizontal plates. In this study, fluid’s dynamic viscosity and thermal conductivity parameters have been assumed to vary depending on temperature changes. The density has been assumed to be incompressible and also the study assumes that the gravitational effects are negligible. The governing equations: continuity, Navier-Stokes, Energy, Magnetic Induction and Concentration equations for the non-Newtonian nanofluid flow have been developed and non-dimensionalized. Dimensionless parameters arising from the dimensionless equations have also been determined. Finite difference numerical approximation method has been used to approximate the systems of the governing equations in difference form. Profiles for the flow variables have been presented and discussed. Results show that increasing thermophoresis parameter increases the specie concentration while increasing Schmidt number and chemical reaction parameter reduces concentration profiles. Magnetic induction profiles rise with an increase in Reynolds number but declines with an increase in magnetic Prandtl number. Temperature and velocity profiles increase with an increase in Reynolds number. The study of electrically conducting fluids with the consideration of Lorentz force, thermophoresis, viscous dissipation, chemical reaction, variable dynamic viscosity, variable thermal conductivity and magnetic induction is very useful in designing heat and mass transfer appliances. It is also significant in cooling and overheating control systems.