Hi,
I am trying to write a code for VLS equilibrium flash unit, the system is a binary system of CH4 and CO2 at high pressures over and low temperatures, a solid phase of CO2 should be present, however, while the code has complied, the steady-state solution fails every time. May you please help me in fixing this issue. The code is attached to this question