for baller=0.2:0.0001:0.6

% E2      E3  E5  E6  E7  V2  V3  V4  V5  V6  V7  H3  H4  
A=[1       0   0   0  -1   0   0   0   0   0   0   0   0;  
 -0.55     1   0   0   0   0   0   0   0   0   0   0   0;
   0      -1   1   0   0   0   0   0   0   0   0   0   0;
   0       0 0.03 -1   0   0   0   0   0   0   0   0   0;
   0       0 0.97  0  -1   0   0   0   0   0   0   0   0;
   0       0   0   0   0   1   0   0   0   0  -1   0   0;
 0.631558*baller  0   0   0   0   0  -1   0   0   0   0   0   0;
   0       0   0   0   0   0 0.15 -1   0   0   0   0   0;
   0       0   0   0   0   0 0.85  0  -1   0   0   0   0;
   0       0   0   0   0   0   0   0 0.92 -1   0   0   0;
   0       0   0   0   0   0   0   0 0.08  0  -1   0   0;
 0.368442*baller  0   0   0   0   0   0   0   0   0   0  -1   0;
   0       0   0   0   0   0   0   0   0   0   0  -1   1;];

b = [2500 0   0   0   0   0   0   0   0   0   0   0   0];

j = eliminate(A, b);

E1 = 0; E4 = 0; V1 = 0; H1 = 0; H2 = 0; H5 = 0; H6 = 0; H7 = 0;
E2 = j(1);
E3 = j(2);
E5 = j(3);
E6 = j(4);
E7 = j(5);
V2 = j(6);
V3 = j(7);
V4 = j(8);
V5 = j(9);
V6 = j(10);
V7 = j(11);
H3 = j(12);
H4 = j(13);
Rec = E7 + V7 + H7;

format long g;

% clc
% disp('The flow rate of E1 is:');
% disp(E1);
% disp('The flow rate of E2 is:');
% disp(E2);
% disp('The flow rate of E3 is:');
% disp(E3);
% disp('The flow rate of E4 is:');
% disp(E4);
% disp('The flow rate of E5 is:');
% disp(E5);
% disp('The flow rate of E6 is:');
% disp(E6);
% disp('The flow rate of E7 is:');
% disp(E7);
% disp('The flow rate of V1 is:');
% disp(V1);
% disp('The flow rate of V2 is:');
% disp(V2);
% disp('The flow rate of V3 is:');
% disp(V3);
% disp('The flow rate of V4 is:');
% disp(V4);
% disp('The flow rate of V5 is:');
% disp(V5);
% disp('The flow rate of V6 is:');
% disp(V6);
% disp('The flow rate of V7 is:');
% disp(V7);
% disp('The flow rate of H1 is:');
% disp(H1);
% disp('The flow rate of H2 is:');
% disp(H2);
% disp('The flow rate of H3 is:');
% disp(H3);
% disp('The flow rate of H4 is:');
% disp(H4);
% disp('The flow rate of H5 is:');
% disp(H5);
% disp('The flow rate of H6 is:');
% disp(H6);
% disp('The flow rate of H7 is:');
% disp(H7);

%disp('The flow rate of the Recycle stream is:');
%disp(Rec);
%disp('at');
%disp(baller*100);

hold on

plot(baller*100,Rec);

end



