The present paper investigates the steady laminar development flow of a Newtonian incompressible electrically conducting fluid in a plane channel subjected to a uniform transverse external magnetic field. In order to solve the problem, the analytical integral momentum method and numerical finite volume method (FVM) are conducted. In integral momentum solution, the Pohlhausen’s fourth-degree velocity profile is assumed. Applying the boundary conditions leads to the appearance of a new dimensionless parameter named auxiliary parameter (AP). It is shown that the value of AP depends on the magnitude of magnetic intensity and the pressure gradient. For each specified value of AP, distinct velocity distribution and as a result, a correlation for computing the magnetic development length is obtained. Besides, the FVM is applied to solve the same problem and the correlation for estimating the magnetic development length is obtained. It can be concluded that for a particular engineering application of magnetohydrodynamics channel flows, the associated value of AP can be determined by the numerical investigation. The results of numerical solution are compared with the analytical results and it can be concluded that the results for AP = − 6 are in agreement with the numerical results. Further, the effect of Hartmann and Reynolds number on the magnetic development length and Lorentz force are illustrated. It was shown that as the Hartmann number increases, the value of development length and Lorentz force are reduced.