For Preisach model, Everett function from the transient curves is needed to simulate the hysteresis phenomena. However it becomes very difficult to get the function if the it would be made only from experiments. In this paper, a simple and stable procedure using least square method and logarithm function to determine the Everett function which follows the Gauss distribution for interaction field axis is proposed. The characteristics of the parameters used in this procedure are also presented. The proposed method is applied to implement hysteresis loops. The simulation for hysteresis loop is compared with experiments and good agreements could be shown.