A method is proposed for analytical approximation of the Preisach distribution function. The Preisach functions of several materials are determined from measured Everett integrals and the features of parameter determination by mean square error minimization over the major or a symmetrical minor loop investigated. The parameters of a model that relaxes the congruency property of the classical Preisach model are also identified.