This paper describes the finite element procedure including the magnetic hysteresis phenomena. The magnetization-dependent Preisach model is employed to simulate the magnetic hysteresis and applied to each elements. Magnetization is calculated by the Fibonacci search method for the applied field in the implementation of the magnetization-dependent model. This can calculate the magnetization very accurately with small iteration numbers. The magnetic field intensity and the magnetization corresponding to the magnetic flux density obtained by the finite element analysis(FEA) are computed at the same time under the condition that these balues must satisfy the constitutive equation. In order to reduce the total calculation cost, pseudo-permeability is used for the input for the FEA. It is found that the presented method is very useful in combining the hysteresis model with the finite element method