Použijeme zadané Antoineovy rovnice

> rov1:=log[10](p1)=6.40627-1702.80/(t+203.);

> rov2:=log[10](p2)=6.22069-1732.22/(t+200.3);

jejímž řešením dostáváme tlaky nasycených par obou složek (převedené na Pa)

> p1:=1000*solve(rov1,p1);

> p2:=1000*solve(rov2,p2);

Dosadíme do vztahu pro celkový tlak vyplývající Raoultova zákona

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

a dosadíme složení kapalné fáze

> x2:=1-x1;

> x1:=0.9;

Závislost tlak nad kapalinou na teplotě lze znázornit následujícím grafem

> plot(p,t=150..200);

Pro zadaný tlak pak vypočteme hledanou teplotu varu

> pv:=100000;

> t:=solve(p=pv,t);

a složení koexistující parní fáze

> y1:=p1*x1/p;

> y2:=1-y1;

>