The behavior of a water-based nanofluid containing motile gyrotactic micro-organisms passing an isothermal nonlinear stretching sheet in the presence of a non-uniform magnetic field is studied numerically. The governing partial differential equations including continuity, momentums, energy, concentration of the nanoparticles, and density of motile micro-organisms are converted into a system of the ordinary differential equations via a set of similarity transformations. New set of equations are discretized using the finite difference method and have been linearized by employing the Newton’s linearization technique. The tri-diagonal system of algebraic equations from discretization is solved using the well-known Thomas algorithm. The numerical results for profiles of velocity, temperature, nanoparticles concentration and density of motile micro-organisms as well as the local skin friction coefficient Cfx, the local Nusselt number Nux, the local Sherwood number Shx and the local density number of the motile microorganism Nnx are expressed graphically and described in detail. This investigation shows the density number of the motile micro-organisms enhances with rise of M, Gr/Re2, Pe and Ω but it decreases with augment of Rb and n. Also, Sherwood number augments with an increase of M and Gr/Re2, while decreases with n, Rb, Nb and Nr. To show the validity of the current results, a comparison between the present results and the existing literature has been carried out.