NRTL
ln(γi) = (summation from j=1 to 3 of [xj × τji × Gji] / summation from k=1 to 3 of [xk × Gki]) + summation from j=1 to 3 of [xj × Gij / summation from k=1 to 3 of [xk × Gkj]] × {τji - summation from m=1 to 3 of [xm × τmj × Gmj] / summation from k=1 to 3 of [xk × Gkj]}
where Gkj = exp(-αij × τij)
and αij and τij are binary NTRL parameters at a given temperature.
τ12 = binary parameter c1
τ21 = binary parameter c2
α12 = binary parameter k
α21 = binary parameter k
τii = 0
αii is not considered.