clc, clf

xd = 0.85
xf = 0.19
xw = 0.01
q  = 1
RR = 4
info_level = 2

//equilibrium data
xeq = [0.0 0.016 0.0315 0.06 0.0855 0.1465 0.206 0.236 0.3495 0.4675 0.4875 0.58 0.6525 0.7 0.7175 0.789 0.842 0.8749 0.8967 0.9485 0.9727 1.0]
yeq = [0.0 0.147 0.2505 0.3765 0.43 0.5005 0.5415 0.56 0.5945 0.641 0.6425 0.689 0.725 0.7495 0.768 0.8111 0.8488 0.8768 0.8973 0.944 0.9692 1.0]
VLE = [xeq;yeq]

//plot(xeq,yeq), return

mct_main(xd,xf,xw,q,RR,VLE,info_level)
