      modelKOfixed = changeRxnBounds(modelKO,BOF,FBAsol0.f-0e-0,'l'); %fix the growth rate      modelKOfixed = changeObjective(modelKOfixed,target); %set target as the objective      FBAsol0 = optimizeCbModel(modelKOfixed,'min',0,true); %find minimum target rate at this growth rate      