I write a three-phase flash calculation for water/pentane/Co2 which has three phases, including vapor, water-rich liquid, and hydrocarbon-rich phases, but this program works up to 12 bar. When the pressure is more than 12 bar, it does not work and converge. It is also very dependent on initial guesses of K values. What is the problem with this program? (I do calculations based on the attached article and my equation of state is VPT EOS)