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]
 

% Call to the ODE function open loop (no negative feedback) 
[t,r] = ode45(@lacIOpen, tspan, r0, []) 

figure(1)
plot(t, r(:,1), t, r(:,2))
legend('lacI mRNA unreg ', 'lacI protein unreg')
xlabel('time in min') 
ylabel('concentration in number per cell')
title('LacI mRNA and Concentration, no negative feedback')
hold on; 



