In this study we analyzed the momentum and heat transfer characteristics of MHD boundary layer flow over an exponentially stretching surface in porous medium in the presence of radiation, non-uniform heat source/sink, external pressure and suction/injection. Dual solutions are presented for both suction and injection cases. The heat transfer analysis is carried out for both prescribed surface temperature (PST) and prescribed heat flux (PHF) cases. The governing equations of the flow are transformed into system of nonlinear ordinary differential equations by using similarity transformation and solved numerically using bvp4c Matlab package. The impact of various non-dimensional governing parameters on velocity, temperature profiles for both PST and PHF cases, friction factor and rate of heat transfer is discussed and presented with the help of graphs and tables. Results indicate that dual solutions exist only for certain range of suction or injection parameters. It is also observed that the exponential parameter have tendency to increase the heat transfer rate for both PST and PHF cases.