The Preisach model needs a distribution function or Everett function to simulate the hysteresis phenomena. To obtain these functions, many experimental data obtained from the first order transition curves are usually required. In this paper, a simple procedure to determine the Preisach density function using the Gaussian distribution function and genetic algorithm is proposed. The Preisach density function for the interaction field axis is known to have Gaussian distribution. To determine the density and distribution, genetic algorithm is adopted to decide the Gaussian parameters. With this method, just basic data like the initial magnetization curve or saturation curves are enough to get the agreeable density function. The results are compared with experimental data and we got good agreements comparing the simulation results with the experiment ones.