In this article, the effects of swimming gyrotactic microorganisms for magnetohydrodynamics nanofluid using Darcy law are investigated. The numerical results of nonlinear coupled mathematical model are obtained by means of Successive Local Linearization Method. This technique is based on a simple notion of the decoupling systems of equations utilizing the linearization of the unknown functions sequentially according to the order of classifying the system of governing equations. The linearized equations, that developed a sequence of linear differential equations along with variable coefficients, were solved by employing the Chebyshev spectral collocation method. The convergence speed of the SLLM technique can be willingly upgraded by successive applying over relaxation method. The comparison of current study with available published literature has been made for the validation of obtained results. It is found that the reported numerical method is in perfect accord with the said similar methods. The results are displayed through tables and graphs.