summaryrefslogtreecommitdiff
path: root/632/CH7/EX7.14/example7_14.sce
blob: e1d9b82caade14ecdac0dd86ee45d37378e1c52a (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
//clc()
//the three phase temperature is first find out, which comes to be 342K, the corresponding Ps1 = 71.18, Ps2 = 30.12
T = [342 343 348 353 363 373];
Ps2 = [30.12 31.06 37.99 47.32 70.11 101.3];
Ps1 = [71.18 72.91 85.31 100.5 135.42 179.14];
P = 101.3;//kPa
for i = 1:4
    y1(i) = 1 - (Ps1(i))/P;
end
for i = 1:6
    y2(i) = 1 - (Ps2(i))/P;
end
plot2d(y2,T);
plot2d(1-y1,T,rect = [0,320,1,380]);
xtitle('Temperature - compositon diagram','x, y (mole fraction of benzene)','Temperature')