%ej2_fsolve

% 1 propano
% 2 n-butano
% 3 n-pentano
% 4 n-hexano

F=100 %kmol/h
z=[0.25,0.25,0.25,0.25] % fracci´on molar
P=5.4 %bar
T=273.15+65 % K
%Datos para la ecuaci´on de Wagner

propfis=[
 % T_C   P_C      A          B        C        D    T_min  T_max
  369.8  42.5  -6.72219  1.33236  -2.13868  -1.38551  145  369.8 % propano
  425.2  38.0  -6.88709  1.15157  -1.99873  -3.13003  170  425.2 %n-butano 
  469.7  33.7  -7.28936  1.53679  -3.08367  -1.02456  195  469.7 %n-pentano
  507.5  30.1  -7.46765  1.44211  -3.28222  -2.50941  220  507.5 %n-hexano
];

for i=1:4
    TC=propfis(i,1);
    PC=propfis(i,2);
    A=propfis(i,3);
    B=propfis(i,4);
    C=propfis(i,5);
    D=propfis(i,6);
    Tr=T/TC;
    Psat(i)=PC*exp((A*(1-Tr) + B*(1-Tr)^1.5 + C*(1-Tr)^3 + D*(1-Tr)^6)/Tr);
    K(i)=Psat(i)/P;
end
 Psat, K
 
 % Funci´on de Rachford-Rice
 
 psi=0
 ((K-1)./((K-1)*psi+1))*z'
 
 psi=1
 ((K-1)./((K-1)*psi+1))*z'
 
 % Soluci´on
 
 psi_0=0.5

 %Los datos de la alimentaci´on (F,z), as´i como las condiciones finales se
 %pasan como par´ametros tal como se indica en la llamada a fsolve.
 
 psi=fsolve(@ej2_ft,psi_0,[],F,z,T,P)
 
 V=psi*F
 L=F-V
 
 x=z./((K-1)*psi+1)
 sum(x)
 y=K.*x
 sum(y)