0001 function objectiveAbbr=checkObjective(model)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015 objRxnInd=find(model.c~=0);
0016 objectiveAbbr=model.rxns{objRxnInd};
0017 if isempty(objRxnInd)
0018 warning('There is no objective!')
0019 else
0020 fprintf('%s\t%s\t%s\t%s\t\t\t%s\n','Coefficient','Metabolite','#','Reaction','#')
0021 for n=1:length(objRxnInd)
0022 objMetInd=find(model.S(:,objRxnInd(n)));
0023 for m=1:length(objMetInd)
0024 Sij=model.S(objMetInd(m),objRxnInd(n));
0025 if length(model.mets{objMetInd(m)})<4
0026 fprintf('%6.4g\t\t%s\t\t\t%i\t%s\t%i\n',Sij,model.mets{objMetInd(m)},objMetInd(m),model.rxns{objRxnInd(n)},objRxnInd(n))
0027 else
0028 if length(model.mets{objMetInd(m)})<8
0029 fprintf('%6.4g\t\t%s\t\t%i\t%s\t%i\n',Sij,model.mets{objMetInd(m)},objMetInd(m),model.rxns{objRxnInd(n)},objRxnInd(n))
0030 else
0031 if length(model.mets{objMetInd(m)})<12
0032 fprintf('%6.4g\t\t%s\t%i\t%s\t%i\n',Sij,model.mets{objMetInd(m)},objMetInd(m),model.rxns{objRxnInd(n)},objRxnInd(n))
0033 end
0034 end
0035 end
0036 end
0037 end
0038 end
0039