Přšdpokládáme, že relativní těkavost (poměr p1/p2) je konstantní

> p1:=2*p2;

Tlak a složení parní fáze výpočteme z Raoultova zákona

> p:=p1*x1+p2*x2;

> y1:=factor(p1*x1/p);

> x2:=1-x1;

Zadáme složení směsi u paty kolony

> x:=0.05;

a postupně vždy vypočteme složení parní fáze, které použijeme jako složení

kapaliny na dalším patře kolony

> i:=0: X[i]:=x:
while x<0.9 do
x:=subs(x1=x,y1);
i:=i+1:
X[i]:=x:
print(i,x);
od:

Průběh destilace podél kolony pak můžeme také graficky znázornit

> n:=i; i:='i';

> s:=[seq([[X[i],X[i]],[X[i],X[i+1]],[X[i+1],X[i+1]]],i=0..n-1)];

> with(plots):

> P1:=plot({x1,y1},x1=0..1,color=[red,black],scaling=CONSTRAINED,axes=BOXED,thickness=2):

> P2:=plot(s,color=green,scaling=CONSTRAINED,axes=BOXED,thickness=2):

> display([P1,P2],labels=['x1','y1']);

>