  This contribution is devoted to cover some technical aspects related to the
use of the recently proposed energy probability distribution zeros in the study
of phase transitions. This method is based on the partial knowledge of the
partition function zeros and has been shown to be extremely efficient to
precisely locate phase transition temperatures. It is based on an iterative
method in such a way that the transition temperature can be approached at will.
The iterative method will be detailed and some convergence issues that has been
observed in its application to the 2D Ising model and to an artificial spin ice
model will be shown, together with ways to circumvent them.
