  The use of partition function zeros in the study of phase transition is
growing in the last decade mainly due to improved numerical methods as well as
novel formulations and analysis. In this paper the impact of different
parameters choice for the energy probability distribution (EPD) zeros recently
introduced by Costa et al is explored in search for optimal values. Our results
indicate that the EPD method is very robust against parameter variations and
only small deviations on estimated critical temperatures are found even for
large variation of parameters, allowing to obtain accurate results with low
computational cost. A proposal to circumvent potential convergence issues of
the original algorithm is presented and validated for the case where it occurs.
