The steady two dimensional magnetohydrodynamic (MHD) boundary layer flow and heat transfer over a stretching/shrinking permeable wedge is numerically investigated. The partial differential equations governing the flow and heat transfer are transformed into a system of ordinary differential equations using a similarity transformation. These equations are then solved numerically using the boundary value problem solver, bvp4c in Matlab software. It is found that dual solutions exist for a certain range of the shrinking strength. A stability analysis is performed to identify which solution is stable and physically reliable.