The paper presents a fixed point iterative method for the magnetic field computation in hysteretic media. The method is applied to a 2D field computation by the finite element method, hysteresis being modeled by a simple analytical scalar model. Convergence of the computation is verified. A test problem, unsolved by classical methods, is treated using both the scalar and vector magnetic potentials. The obtained results show the algorithm efficiency, which is then applied to magnetic recording modeling.