
62 y=[y y9];
63 x10=1.0;
64 x=[x x10];
65 temp0=inline(‘1.0.*(10.^(4.72583-(1660.652./(T-
1.461))))+(1-1.0).*(10.^(4.08245-(1346.382/(T-
53.508))))-1.013’);
66 temp=fzero(temp0,initval);
67 y10=x10.*(10.^(4.72583-(1660.652./(temp-1.461))))/
(x10.*(10.^(4.72583-(1660.652./(temp-1.461))))+(1-
x10).*(10.^(4.08245-(1346.382/(temp-53.508)))));
68 y=[y y10];
69 fprintf(‘x \t\t\t y \t\t\t’)
70 for i=1:11
71 fprintf(‘\n%f \t %f \t %f \n’,x(i),y(i))
72 end
73 fprintf(‘\n’);
74 diagonal =(0:1);
75 plot(x,y,diagonal,diagonal)
76 grid on
77 axis square
78 axis equal
79 xlabel(‘Benzene in Liquid (mole fraction)’);
80 ylabel(‘Benzene in vapor (mole fraction)’);
81 title(‘Benzene ...