The problem of an unsteady magnetohydrodynamic stagnation point flow of an incompressible viscous fluid over a shrinking sheet is discussed in the presence of thermal radiation and boundary slip, which has not been documented to date in the literature. The governing boundary-layer equations are transformed to high order nonlinear and ordinary differential equations by similarity transformations and then solved numerically. The effects of magnetic parameter, unsteadiness parameter, radiation parameter, velocity, and thermal slip parameters on velocity and temperature are analyzed and discussed. It is found that dual solutions of both velocity and temperature fields exist for negative values of the velocity ratio parameter. The results indicate that dual solution domains of velocity and temperature expand as the unsteadiness parameter increases. C⃝ 2015 Wiley Periodicals, Inc. Heat Trans Asian Res, 00(00): 1-16, 2015; Published online in Wiley Online Library (wileyonlinelibrary.com/journal/htj).