The shear stress of the steady two-dimensional boundary layer flow of a hydromagnetic flow due to a radially stretching surface is investigated. The boundary layer equations governing the flow are transformed into a singular equation by using suitable equivalent transformations. The equation is then turned to nonlinear equations by using Galerkin finite element method. At last, the numerical solutions for the nonlinear equations are estimated through Newton iterative method. It is obtained the shear stress of this fluid corresponding to the parameter M different values. Moreover, the results are compared with previous conclusions through table. It’s shown that the numerical results and previous solution is consistent. This means that the Galerkin finite element method is a good method to solve the hydromagnetic boundary layer flow.