1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
|
clear;
clc;
//To find Approx Value
function[A]=approx(V,n)
A=round(V*10^n)/10^n;//V-Value n-To what place
funcprot(0)
endfunction
//Example 10.6
//Caption : Program To Find the Composition of Vapor and Liquid Phases
//Given Values
z=[0.1,0.2,0.7];
K=[10,1.76,0.52];
V_01=[0.1:0.1:1]
for(i=1:10)
for(j=1:3)
y_01(i,:)=approx((z.*K)./(1+(V_01(i)*(K-1))),3);
end
Sum(i)=sum(y_01(i,:));
end
for(i=1:10)
if(Sum(i+1)<1)
a1=V_01(i);
i=-1;
break;
end
end
V_001=[a1:0.01:a1+0.1]
for(i=1:10)
for(j=1:3)
y_001(i,:)=approx((z.*K)./(1+(V_001(i)*(K-1))),3);
end
Sum(i)=sum(y_001(i,:));
end
for(i=1:10)
if(Sum(i+1)<1)
a01=V_001(i);
i=-1;
break;
end
end
V_0001=[a01:0.001:a01+0.01];
for(i=1:10)
for(j=1:3)
y_0001(i,:)=approx((z.*K)./(1+(V_0001(i)*(K-1))),3);
end
Sum(i)=sum(y_0001(i,:));
end
for(i=1:10)
if(Sum(i+1)<1)
a001=V_0001(i);
i=-1;
break;
end
end
V=a001;
y_02=approx((z.*K)./(1+(a1*(K-1))),3);
y_03=approx((z.*K)./(1+((a1+0.1)*(K-1))),3);
y_0273=approx((z.*K)./(1+(V*(K-1))),3);
x_0273=approx(y_0273./K,3);
Ans=[z',K',y_02',y_03',y_0273',x_0273';1,0,sum(y_02),sum(y_03),sum(y_0273),approx(sum(x_0273),2)];
disp(Ans,' z K y,V=0.2 y,V=0.3 y,V=0.273 x,V=0.273')
disp('NOTE : Last Row represents the summation')
disp('Hence for V = 0.273 E(yi) = 0 and E(xi) = 0')
disp(V,'Fraction of Vapor (V) = ')
//End
|
