The Preisach model neds a density function or Everett function for the hysterisis operator to simulate the hysteresis phenomena. To obtain the function, many experimental data for the first order transition curves are required. However, it needs so much efforts to measure the curves, especially for the hard magnetic materials. By the way, it is well known that the density function has the Gaussian distribution for the interaction axis on the Preisach plane. In this paper, we propose a simple technique to determine the distribution function or Everett function analytically. The initial magnetization curve is used for the distribution of the Everett function for the coercivity axis. A major, minor loop and the initial curve are used to get the Everett function for the interaction axis using the Gaussian distribution function and acceptable results were obtained.