This paper investigates the steady, two dimensional, and magnetohydrodynamic flow of copper and alumina/water hybrid nanofluid on a permeable exponentially shrinking surface in the presence of Joule heating, velocity slip, and thermal slip parameters. Adopting the model of Tiwari and Das, the mathematical formulation of governing partial differential equations was constructed, which was then transformed into the equivalent system of non-linear ordinary differential equations by employing exponential similarity transformation variables. The resultant system was solved numerically using the BVP4C solver in the MATLAB software. For validation purposes, the obtained numerical results were compared graphically with those in previous studies, and found to be in good agreement, as the critical points are the same up to three decimal points. Based on the numerical results, it was revealed that dual solutions exist within specific ranges of the suction and magnetic parameters. Stability analysis was performed on both solutions in order to determine which solution(s) is/are stable. The analysis indicated that only the first solution is stable. Furthermore, it was also found that the temperature increases in both solutions when the magnetic parameter and Eckert number are increased, while it reduces as the thermal slip parameter rises. Furthermore, the coefficient of skin friction and the heat transfer rate increase for the first solution when the magnetic and the suction parameters are increased. Meanwhile, no change is noticed in the boundary layer separation for the various values of the Eckert number in the heat transfer rate.