A block method of the Adams type is employed to obtain a numerical solution for the boundary layer problems arising from agnetohydrodynamics (MHD) flow. The governing equations are modeled as a set of third and second order differential equations subjected to boundary condition. The block method is implemented to provide efficient numerical approximation of the boundary layer equations. Numerical experiments are carried out on three test problems, and the solutions obtained using the method are compared with existing solutions in the literature.