From 2009.igem.org
model = changeRxnBounds(model,{'EX_glc(e)','EX_gly(e)','EX_ala_D(e)','EX_ala_L(e)','EX_arg_L(e)','EX_asn_L(e)','EX_asp_L(e)','EX_cys_D(e)','EX_cys_L(e)','EX_gln_L(e)','EX_glu_L(e)','EX_his_L(e)','EX_ile_L(e)','EX_leu_L(e)','EX_lys_L(e)','EX_met_D(e)','EX_met_L(e)','EX_phe_L(e)','EX_pro_L(e)','EX_ser_D(e)','EX_ser_L(e)','EX_thr_L(e)','EX_trp_L(e)','EX_tyr_L(e)','EX_val_L(e)','EX_adn(e)','EX_gsn(e)','EX_cytd(e)','EX_thymd(e)','EX_uri(e)','EX_thm(e)','EX_pnto_R(e)','EX_xyl_D(e)','EX_skm(e)','EX_26dap_M(e)'},[-99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 -99],'l');
solution=optimizeCbModel(model);
fValue = solution.f
"
