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

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


xd = 0.99
xf = 0.5
xw = 0.01
q  = 0.5
R  = 3.5

info_level = 1

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