In magnetic systems, distribution of electromagnetic force density causes mechanical deformation, which results in noise and vibration. In this paper, Korteweg-Helmholtzs energy method and equivalent magnetic charge method are employed for comparison of their resulting distributions of force density. The force density from the Korteweg-Helmholtzs method is expresses with two Maxwell stresses on the inside and the outside fo magnetic material respectively. The other is calculated using the magnetic Coulombs law. In the numerical model of an electromagnet, their numerical results are compared. The distributions by the two methods are almost the same. And their total forces are also shown to be the same to the one calculated from the conventional Maxwell stress tensor. But the magnetic charge method is easier and more efficient in numerical calculation.