initCobraToolbox;

% calculate metabolic flux in standard conditions (maximal growth)
model = readCbModel('C:\Users\Public\gb-2009-10-6-r69-s4');
solution=optimizeCbModel(model);
printFluxVector(model, solution.x, true, false);

% calculate metabolic flux during maximum urease activity
model = changeObjective(model, {'rxn00101'}, 1);
solution=optimizeCbModel(model);
printFluxVector(model, solution.x, true, false);