1. A method of operating a computer system to create a model of a ferroelectric element for use in a computerized simulation of an electronic circuit including the ferroelectric element, comprising the steps of: measuring polarization behavior of at least one ferroelectric element; from the measured polarization behavior, deriving a probability distribution function of positive and negative coercive voltages at which domains of the ferroelectric element are polarized into polarization states, the probability distribution function represented by a numerical grid of a plurality of gridcells, each gridcell associated with positive and negative coercive voltage values and having a likelihood value assigned thereto; calculating an expected number of domains for each of the plurality of gridcells based on the likelihood value for the gridcell and a number of domains in the modeled ferroelectric element; randomly selecting a number of domains to be assigned to each gridcell for the modeled ferroelectric element according to its expected number of domains; and superposing polarization characteristics associated with each of the selected domains for the modeled ferroelectric element.