blob: 478edc6efcaae9d197472780b3b50b8728332594 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
//clc()
T = [371.4 378 383 388 393 398.6]
Pas = [101.3 125.3 140 160 179.9 205.3]
Pbs = [44.4 55.6 64.5 74.8 86.6 101.3]
Ptotal = 101.3;//kPa
for i = 1:6
x(i) = (Ptotal - Pbs(i))/(Pas(i) - Pbs(i));
end
for i = 1:6
y(i) = x(i) * Pas(i) / Ptotal;
end
w = x;
plot(x,w);
plot(x,y,'-o');
xtitle('Equilibrium curve','x, mole fraction in liquid','y, mole fraction in vapour');
|
