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

// this example rev.3 incorrectly solves for R min
xeq = [0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0.85,0.90,0.95,1]
yeq = [0,0.3372,0.4521,0.5056,0.5359,0.5589,0.5794,0.5987,0.6177,0.6371,0.6558,0.6765,0.6986,0.7250,0.7550,0.7840,0.8167,0.8591,0.8959,0.9474,1]
VLE = [xeq;yeq]



xd = 0.82
xf = 0.5
xw = 0.05
q  = 0.1
R  = 2

info_level = 1

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