>

Vytvoření seznamů experimentálních dat (P=tlak (MPa), A=adsorb. množství)

> P:=[0.1,0.2,0.28,0.41,0.98,1.39,1.93,2.75,3.01,3.51]; A:=[0.261,0.385,0.501,0.669,1.104,1.322,1.694,2.203,2.243,2.472];

Přepočtení experimentálních dat do tvaru vhodného pro linearizaci

> lnP:=map(x->ln(x),P);lnA:=map(x->ln(x),A);

>

Lineární korelace metodou nejmenších čtverců (závislost ln a na ln p )

> with(stats):

> fit[leastsquare[[x,y],y=a*x+b, {a,b}]]([lnP,lnA]);

> assign(%);

Přiřazení hodnoty směrnice symbolu K a posunutí symbolu Q

> K:=lcoeff(y); Q:=tcoeff(y);

Výsledek (parametry Freundlichovy izotermy)

> k:=exp(Q); n:=1/K;

Grafické zobrazení získaných dat

> p:=seq([lnP[i],lnA[i]],i=1..nops(P));

> plot([y,[p]], x = -2..2, color=[red,blue], style=[line,point], labels=["ln p", "ln a"] );

Grafické znázornění Freundlichovy izotermy

> aF := k*pF^(1/n);

> plot(aF, pF=0..3);

>

>

>