disp('datos iniciales')
A1 = [0.55 0 0 0 0 0.10 -1 0];
A2 = [1 -1 0 0 0 1 0 0];
A3 = [0 0 -1 0 0 0 1 -1];
A4 = [0 1 -1 -1 0 0 0 0];
A5 = [0 0 1 0 0 0 -2/3 0];
A6 = [0 0 0 0 -0.1 -0.1 0 1];
A7 = [0 0 0 1 -1 -1 0 0];
A8 = [0 0 0 0 1 0 0 0];
A = [A1;A2;A3;A4;A5;A6;A7;A8] 
b = [0; 0; 0; 0; 0; 0; 0; 1]
disp('resolucion directa')
x=A\b
%*********Factorizando a LU**********
disp('resolucion con factorizacion LU')
[L,U,P]=lu(A);
% PAx=LUx ---> con Ax=b----> Pb=LUx 
b2=P*b; %Queda LUx=b2
%Ux=y ----> Ly=b2
y=L\b2
x=U\y
%********************************

%*****variacion del reciclo**********
y = []; % vector vacio
for k=0.5:0.45/10:0.95, % Bucle de 10 pasos para variar lo que se recicla
A(5,7) = -k; % variacion del reciclo en el splitter
x=A\b; % resolucion del sistema lineal
R = x(6,1)/x(1,1); % calculo de la relacion de reciclo
y = [y R]; % alamcena resultado en el vector y
end
k=0.5:0.45/10:0.95; % eje x para el grafico
plot(k,y)
title('grafico del reciclo')


