  We report on a new method to extract thermodynamic properties from the
density of partition function zeroes on finite lattices. This allows direct
determination of the order and strength of phase transitions numerically.
Furthermore, it enables efficient distinguishing between first- and
second-order transitions, elucidates crossover between them and illuminates the
origins of finite-size scaling. The power of the method is illustrated in
typical applications for both Fisher and Lee-Yang zeroes.
