clc
clear

exec(pwd()+'/mct_main.sci',-1)


xd = 0.95
xf = 0.5
xw = 0.1
q  = .5
RR = 3.5

// generate ideal VLE data
alpha = 1.7 //relative volatility
xeq   = linspace(0,1,200)
yeq   = alpha*xeq./(1+(alpha-1).*xeq)
VLE   = [xeq;yeq]


info_level = 1

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