In the present study, squeezing flow of micropolar nanofluid between parallel infinite disks in the presence of magnetic field perpendicular to plane of the disks is taken into account. The constitutive equations that govern the flow configuration are converted into nonlinear ordinary differential with the help of suitable similarity transforms. HAM package BVPh2.0 has been employed to solve the nonlinear system of ordinary differential equations. Effects of different emerging parameters like micropolar parameter K, squeezed Reynolds number R, Hartmann number M, Brownian motion parameter Nb, thermophoresis parameter Nt, Lewis number Le for dimensionless velocities, temperature distribution and concentration profile are also discussed graphically. In the presence of strong and weak interaction (i.e. n = 0 and n = 0.5), numerical values of skin friction coefficient, wall stress coefficient, local Nusselt number and local Sherwood number are presented in tabulated form. To check the validity and reliability of the developed algorithm BVPh2.0 a numerical investigation is also a part of this study.