clear
load modelo
aula=model

aula=changeRxnBounds(aula,'Biomass_Ecoli_core_w_GAM',0.45,'l');
%aula=changeRxnBounds(aula,'EX_glc(e)',-10,'b');
%aula=changeObjective(aula,'Biomass_Ecoli_core_w_GAM');
aula=changeObjective(aula,'EX_glc(e)');

solucao=optimizeCbModel(aula,'max');
solucao.x;

listaGlc=findRxnsFromMets(aula,{'glc_D[e]'});
IDsGlc=findRxnIDs(aula,listaGlc);
for m=1:length(IDsGlc)
    fluxGlc(m)=solucao.x(IDsGlc(m));
end