%%% Aula - estado estacionario - 05/06/2020

%%% MFA

%%% Ordem metabolitos: G6P, G3P, PIR, AcCoA, CO2, O2, R5P

N = [1 -1 0 0 0 0 0 0 -1 0 1;
     0 1 -1 0 0 0 0 0 0 1 -1;
     0 1 1 -1 0 0 0 0 0 0 0;
     0 0 0 1 -2 -1 0 0 0 0 0;
     0 0 0 1 0 2 -1 0 1 0 0;
     0 -.5 -.5 -.5 0 -2 0 1 0 -1 0;
     0 0 0 0 0 0 0 0 1 -1 0];
     
%%% exemplo 1
%id_m = [1; 5];          % indices dos fluxos medidos em N (colunas)
%id_d = [2; 3; 4; 6];    % indices dos fluxos a serem determinados em N (colunas)
%
%v_m = [1.97; 1.6];
%
%Nm = N(1:4,id_m);      % Modiquei para funcionar com a matriz completa 
%Nd = N(1:4,id_d);      % considerei so os 4 primeiros metabolitos como fizemos em aula com a matriz menor
%
%v_d = -Nd^-1*Nm*v_m;

%%% exemplo 2
%id_m = [1; 5; 7];
%id_d = [2; 3; 4; 6];
%
%v_m = [1.97; 1.6; 4.26];
%
%Nm = N(1:5,id_m);     % Modiquei para funcionar com a matriz completa 
%Nd = N(1:5,id_d);     % considerei so os 5 primeiros metabolitos como fizemos em aula com a matriz menor
%
%v_d = -Nd\Nm*v_m;
%
%v_m_est = [v_d(1);
%           (v_d(3) - v_d(4))/2;
%            v_d(3) + 2*v_d(4)];
%
%bal_c = (v_m_est(2)*4 + v_m_est(3)*1)/(v_m_est(1)*6)

%%% exercicio 1
%id_m = [1; 5; 7; 8];
%id_d = [2; 3; 4; 6];
%
%v_m = [1.97; 1.6; 4.26; 3.57];
%
%Nm = N(1:6,id_m);       % Modiquei para funcionar com a matriz completa 
%Nd = N(1:6,id_d);       % considerei so os 6 primeiros metabolitos como fizemos em aula com a matriz menor
%
%v_d = -Nd\Nm*v_m;
%
%v_m_est = [v_d(1);
%           (v_d(3) - v_d(4))/2;
%            v_d(3) + 2*v_d(4)
%            .5*v_d(1) + .5*v_d(2) + .5*v_d(3) + 2*v_d(4)];
%
%bal_c = (v_m_est(2)*4 + v_m_est(3)*1)/(v_m_est(1)*6)

%%% exemplo 3
%id_m = [1; 5; 7; 8];
%id_d = [2; 3; 4; 6; 9; 10; 11];
%
%v_m = [1.97; 1.6; 4.26; 3.57];
%
%Nm = N(:,id_m);
%Nd = N(:,id_d);
%
%v_d = -Nd^-1*Nm*v_m;