%runnning distillation model 'multicomA' for any mixture of NC components
clear all
[t,x]=ode15s('multicomA_lv',[0 20000],0.5*ones(1,4*41));
lengthx=size(x); 
% Vector 'Xss' returns the steady-state molar compositions of components 1 to NC-1 at all NC stages and molar liquid hodups at all NC stages
Xss=x(lengthx(1),:)';
t
