An alternative higher-order Eggshell method is proposed for determining forces acting on a ferromagnetic body in magnetic field. The strongly nonlinear problem is solved by the Newton method effectively implemented in our codes Agros2D and Hermes. The computation exhibits substantially higher parameters (accuracy, time of solution) than the algorithms used even in available professional codes. The methodology is illustrated with a typical example.