clear all 


% This is the main function for problem 1 HW 1 


% Time span to run simulation 
tspan = [0, 120] 

% Initial mRNA and protein concentration 
r0 = [0; 0]
 
figure(2)

% Call to the ODE function closed loop (with negative feedback) 
[t,r] = ode45(@lacIClosed, tspan, r0, []) 

plot(t, r(:,1), '-r', t, r(:,2), '-k')
legend('lacI mRNA autoreg', 'lacI protein autoreg')
xlabel('time in min') 
ylabel('concentration in number per cell')
title('LacI mRNA and Concentration with negative feedback')


