1. Solve: T1sat(P), T2sat(P), T0 = x1T1sat + x2T2sat
2. Table: Temp, A, γ1, γ2, α11, α21, x1γ1α11, x2γ2α21, P1sat, P2sat, y1, y2
3. Iteration 0: Temp = T0 → A0 = f(T0) → γi0 = f(x,T0) → αii = 1, αik = exp(Ai - (Bi/(Ci+T0)) - Ak + (Bk/(Ck+T0))) → x1γ1α11, x2γ2α21 → P1sat = P/(x1γ1α11 + x2γ2α21)
4. Iteration 1: Temp = T(P1sat) → A = f(T) → γi = f(x,T) → αii = 1, αik = exp(Ai - (Bi/(Ci+T)) - Ak + (Bk/(Ck+T))) → x1γ1α11, x2γ2α21 → P1sat = P/(x1γ1α11 + x2γ2α21) → P2sat = f(T)
5. Repeat iterations until convergence
6. yi = 𝛾𝑖xiPisat/P