function BoilingPoint
% 計算混合液體的沸點及平衡蒸汽組成
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/02/26 $

clear all
clc

x = [0.5      0.3      0.2    ];    % 液相組成
a = [7.84135  8.08840  8.11404];
b = [1750     1985     2129   ];
P = 760;    % 大氣壓，mmHg
T0 = 353;   % K

T = fzero(@Equation,T0,[],x,a,b,P);
y = x.*P0(T,a,b)./P;
fprintf('\n\n  Results：')
fprintf('\n\t混合物沸點為：%.2f %s\n',T,'K')
fprintf('\t平衡蒸汽組成為：')
disp(y)

% ------------------------------------------------------------------
function f = Equation(T,x,a,b,P)
y = x.*P0(T,a,b)./P;
f = sum(y) - 1;

% ------------------------------------------------------------------
function f = P0(T,a,b)
f = exp(log(10)*(a-b./T));