For precise modeling of electromagnetic devices, we have to model material hysteresis. A Genetic Algorithm, Differential Evolution with three different strategies, teaching–learning-based optimization and Artificial Bee Colony, were used for testing seven different modified mathematical expressions, and the best combination of mathematical expression and solving method was used for hysteresis modeling. The parameters of the hysteresis model were determined based on the measured major hysteresis loop and first-order reversal curves. The model offers a simple determination of the magnetization procedure in the areas between measured curves, with the only correction of two parameters based on only two known points in the magnetization process. It was tested on two very different magnetic materials, and results show good agreement between the measured and calculated curves. The calculated curves between the measured curves have correct shapes. The main difference between our model and other models is that, in our model, each measured curve, major and reversal, is described with different parameters. The magnetization process between measured curves is described according to the nearest measured curve, and this ensures the best fit for each measured curve. In other models, there is mainly only one curve, a major hysteresis or magnetization curve, used for the determination of the parameters, and all other curves are then dependent on this curve. Results confirm that the evolutionary optimization method offers a reliable procedure for precise determination of the parameters.