Team:Alberta/Modeling/DualDeletionTest

From 2009.igem.org

% Remember to run start.m;
b=0.01;
deltamodel = model;
solution1 = optimizeCbModel(deltamodel, 'max', 'false', 'false');
for a=1:length(model.rxns);
deltamodel = model;
deltamodel = changeRxnBounds(model, model.rxns(a), 0, 'b');
deltasolution = optimizeCbModel(deltamodel);
   if deltasolution.f<b;
       v(a)=model.grRules(a);
   end
end
v = v'; % Copy v and model.rxns into a txt file to compare resulting data
Retrieved from "http://2009.igem.org/Team:Alberta/Modeling/DualDeletionTest"