Modelling of hysteresis in soft magnetic materials is necessary for accurate transient studies of transformers and electrical machines. The phenomena such as magnetization inrush currents and ferroresonance are directly resultants of the hysteresis properties. A proper model for numerical transient analyses should be easy regarding implementation and parameter estimation. In [1], authors introduced a differential hysteresis model derived from the Preisach classic theory. The model has accuracy comparable with Preisach models. Furthermore, in numerical simulations it can be inverted easily. In this paper, a practical method for estimation of the model functions by analytical expressions is suggested. The required input data is decreased to only the major hysteresis loop that is usually available that is easy to measure. Also, details of the numerical implementation are described. The presented method is verified by experimental measurements on several grain oriented and non-oriented electrical steels.