clear all;
clc;

sbioloadproject hw5

noise = getconfigset(m1);
set(noise, 'SolverType', 'ssa', 'StopTime', 100)
simDataObj = sbiosimulate(m1);

%Set RNAP Destruction and Creation Rate
RNAP_Des = sbioselect(m1,'Name','RNAP Destruction');
RNAP_Des.Value = 20;

RNAP_Crea = sbioselect(m1,'Name','RNAP Creation');
RNAP_Crea.Value = 20;

yfp = [];
cfp = [];

for i=1:500,
    simDataObj = sbiosimulate(m1);
    [times, data, names] = selectbyname(simDataObj, {'CFP', 'YFP'});
    yfp = [yfp data(end,1)];
    cfp = [cfp data(end,2)]; 
end

plot(yfp,cfp,'r*')