clear
load modelo2
mutante=model;

mutante = changeRxnBounds(mutante,'EX_glc(e)',-8,'l');
mutante = changeRxnBounds(mutante,'ATPM',8.39,'b');
%mutante = changeRxnBounds(mutante,'Biomass_Ecoli_core_w_GAM',0.1,'l');

mutante = addReaction(mutante,'sintPHB',{'accoa[c]','nadph[c]','h[c]','3HB','coa[c]','nadp[c]'},[-2 -1 -1 1 2 1],false,0,70);
mutante = addExchangeRxn(mutante,{'3HB'},-20,70);

mutante = changeObjective(mutante,'sintPHB');
solucao = optimizeCbModel(mutante,'max');
polimero = solucao.f
cresc = findRxnIDs(mutante,{'Biomass_Ecoli_core_w_GAM'});
taxa = solucao.x(cresc)