In this study, first-order slip effect with viscous dissipation and thermal radiation in micropolar fluid on a linear shrinking sheet is considered. Mathematical formulations of the governing equations of the problem have been derived by employing the fundamental laws of conservations which then converted into highly non-linear coupled partial differential equations (PDEs) of boundary layers. Linear transformations are employed to change PDEs into non-dimensional ordinary differential equations (ODEs). The solutions of the resultant ODEs have been obtained by using of numerical method which is presented in the form of shootlib package in MAPLE 2018. The results reveal that there is more than one solution depending upon the values of suction and material parameters. The ranges of dual solutions are S ≥ S c i , i = 0 , 1 , 2 and no solution is S < S c i where S c i is the critical values of S . Critical values have been obtained in the presence of dual solutions and the stability analysis is carried out to identify more stable solutions. Variations of numerous parameters have been also examined by giving tables and graphs. The numerical values have been obtained for the skin friction and local Nusselt number and presented graphically. Further, it is observed that the temperature and thickness of the thermal boundary layer increase when thermal radiation parameter is increased in both solutions. In addition, it is also noticed that the fluid velocity increases in the case of strong magnetic field effect in the second solution.