“…Yet, as the number of phase is not known in advance, it is tricky to select the phase number relevant set of equations that is to be solved and numerical problem difficulties may occur due to non-smooth problems or pseudo solutions Morbidelli, 1981, 1982 ;Bullard and Biegler, 1993). To overcome this difficulty, equationsolving method usually involve two sequential or imbricated loops : first the determination of the number and type of coexisting phases relying on a stability test like the tangent plane criterion (Baker et al, 1981, Michelsen, 1982a, 1982b, Cairn and Fürzer, 1990a, 1990b, Gupta et al, 1991, Widagdo et al, 1992, Sun and Seider, 1995Jalali and Seader, 2000) and second, the calculation of phases composition and ratio solving the correct set of governing equations. For more than two coexisting phases like liquid-liquid-vapor equilibria, these two loops may be repeated several times until the correct phase configuration is found.…”