function BatchDist
% 簡單蒸餾計算
%
%   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

alpha = 2.5;    % 相對揮發度
x0 = 0.60;      % 開始時的苯含量    
var = 2;        % var = F0/F

[x,fval] = fzero(@DistEq,x0,[],alpha,x0,var)

% ------------------------------------------------------------------
function f = DistEq(x,alpha,x0,var)
f = log(var) - 1/(alpha-1)*(log(x0/x)+alpha*log((1-x)/(1-x0)));