A robust algorithm is outlined for the solution of the classical non-linear isofugacity problem by using equations of state (EOS). The procedure suggested is free of numerical pitfalls from the triple point to the critical point. The combination of differential stability (spinodal curve) and zero pressure reference concepts yields an automatic initialization of the calculations inside a region of convergence where the EOS can always predict pure compound vapor pressures at a given temperature.