The steady two-dimensional magnetohydrodynamic stagnation flow towards a nonlinear stretching surface is studied. The no-slip condition on the solid boundary is replaced with a partial slip condition. A scaling group transformation is used to get the invariants. Using the invariants, a third-order ordinary differential equation corresponding to the momentum is obtained. An analytical solution is obtained in a series form using a homotopy analysis method. Reliability and efficiency of series solutions are shown by the good agreement with numerical results presented in the literature. The effects of the slip parameter, the magnetic field parameter, the velocity ratio parameter, the suction velocity parameter, and the power law exponent on the flow are investigated. The results show that the velocity and shear stress profiles are greatly influenced by these parameters.