The problem of steady laminar mixed convection boundary layer flow along vertical thin needle with variable surface heat, mass and motile microorganism flux in the presence of gyrotactic microorganism is considered in this study. The dimensionless leading equations of continuity, momentum, concentraton and motile microorganism conservation are reduced to ordinary differential equations with the help of similarity transformations. The transformed governing equations are then numerically solved by using MATLAB BVP4C function. The research is reached to excellent argument by comparison in few cases between the results obtained from MATLAB and Maple algorithm with the help of dsolve command. Numerical calculations are carried out for various values of the dimensionless parameters of the problem which includes mixed convection parameter λ, power law index m, buoyancy parameters N1, N2 Lewis parameter Le, bioconvection lewis parameter Lb, Bioconvection peclet number Pe and also the parameter a representing the needle size. It is also shown from the results that the surface (wall) temperature, surface fluid concentration, surface motile microorganism concentration and the corresponding velocity, temperature, concentration and motile microorganism profiles are significantly induced by these parameters. The results are pictured and discussed in detail.