% Simulación de un separador con reciclo

format short g

% Corriente de entrada
f4=[1409.5; 1276.5; 119.1; 36.0; 2.1] % kmol/h
T4=681.3 % C
P4=37.7 %bar

F4=sum(f4)
x4=f4/F4

C4={F4;x4;P4;T4}

% Especificaciones de equipos

x1v101=[0.1;0.1;0.9;0.9;0.1]
x1d6=0.8

% Inicialización

C16i=C4;
v16i=cav(C16i);
VDEP={C4,x1v101,x1d6};

% Newton 
for i=1:6
   res(:,i)=citer(VDEP,v16i);
   [J,FAC]=numjac('citer',VDEP,v16i,citer(VDEP,v16i),1e-10,[],[]);
   deltac=-J\res(:,i);
   v16=v16i+deltac;
   v16i=v16;
end
C16=vac(v16);
CE=mezclador(C4,C16);
[C5,CC]=separador(CE,x1v101);
[C6,C16]=divisor(CC,x1d6);

res
tbalance(C4,CE,C5,CC,C6,C16)

% Sustitución directa
for i=1:6
   CE=mezclador(C4,C16i);
   [C5,CC]=separador(CE,x1v101);
   [C6,C16]=divisor(CC,x1d6);
   res2(:,i)=compara(C16i,C16);
   C16i=C16;
end


res2
tbalance(C4,CE,C5,CC,C6,C16)




