Team:Alberta/Modeling/DualDeletionTest

% 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