Team:DTU-Denmark/plot mutant2

From 2013.igem.org

Timer series plots of Mutant2

clear all
% create model
model=SBmodel('mutant_2.txt')
SBPDmakeMEXmodel(model,'modelmex');

%calculate steady state
SBsteadystate(model)

conc = [5.872e-2 5.872e-1 5.872 29.358];
enz = [2.115e-4 1.0575e-03 2.1115e-3 2.115e-1];
count = 0;
for i=1:4
    count = 0;
    for j=1:4
        count = count+1;
        % make parametervector
        parametervector = makeparamvecSBPD('modelmex',{'NO20','E'},[conc(i); enz(j)]);
        output=modelmex([0:0.001:10],[],parametervector);
        % save results in separate vectors
        % check in which order they are (double click on output in workspace)
        NO(:,count) = output.statevalues(:,1);
        N2O(:,count) = output.statevalues(:,2);
    end
    t = output.time';
    % plot concentrations over time:
    figure
    plot(t,N2O(:,1),t,N2O(:,2),'--',t,N2O(:,3),':',t,N2O(:,4),'-.','LineWidth',2)
    if i == 1
        title('1 mg/l nitrite (58.72 uM)')
    elseif i == 2
        title('10 mg/l nitrite (587.2 uM)')
    elseif i == 3
        title('100 mg/l nitrite (5.872 mM)')
    elseif i == 4
        title('500 mg/l nitrite (29.358 mM)')
    end
    
    xlabel('Time in [min]','FontSize',16)
    ylabel('N2O concentration in [mM]','FontSize',16)
    hleg1 = legend('100 enzymes per cell','500 enzymes per cell','1000 enzymes per cell','100 000 enzymes per cell');
    set(hleg1,'Location','SouthEast','FontSize',12)
end