%running distillation model 'colamod_wilson' for a ternary mixture

% Arbitrary initial conditions 
[t,x]=ode15s('cola_lv_wilson',[0 20000],0.5*ones(1,3*41));
lengthx=size(x); 

% Vector 'Xinit' returns molar compositions of component A and B in each stage and liquid molar holdups in each stage  
Xinit=x(lengthx(1),:)';
