function g=ELV2Wilson(w)
global A
global B
global C
global k
global P
global y
global Ren
% La variable w contiene las incógnitas w=[T x_1]
T=w(1);    % T en grados kelvin para consistencia dimensional
xx1=w(2);
x=[xx1 1-xx1];
% Escribimos primero las ecuaciones auxiliares
% Para presión de vapor usamos el modelo de Antoine
Pvap=10.^(A-B./(T+C));% mmHg
% Para el coeficiente de actividad usamos el modelo de Wilson
Lambda12=k(4)./k(3).*exp(-k(1)./(Ren.*T));
Lambda21=k(3)./k(4).*exp(-k(2)./(Ren.*T));
gamma(1)=exp(-log(x(1)+Lambda12.*x(2))+x(2).*(Lambda12./(x(1)+Lambda12.*x(2))-Lambda21./(Lambda21.*x(1)+x(2))));
gamma(2)=exp(-log(x(2)+Lambda21.*x(1))-x(1).*(Lambda12./(x(1)+Lambda12.*x(2))-Lambda21./(Lambda21.*x(1)+x(2))));
% Ahora escribimos las ecuaciones a resolver
g(1)=x(1).*gamma(1).*Pvap(1)-y(1).*P;
g(2)=x(2).*gamma(2).*Pvap(2)-y(2).*P;

Error using ELV2Wilson (line 10)
Not enough input arguments.

Published with MATLAB® R2012b