function [ dRdt ] = lacIClosed( t, R )


am = 24  % transcription rate in bp/min/cell, Bionumbers ID=100060
dm = 0.2  % degradation rate in /min/cell 	BioNumbers ID= 100060 
p = 2    % translation rate in lacI mRNA/min/cell 	BioNumbers ID=100060
d = 0.18  % degradation rate of protein in /min/cell, 
         % estimated to be much larger than mRNA degradation (not from bionumbers)

K = 1

amp = am*(K + ((p*am)/(dm*d))/K)

% Dyamics for LacI mRNA
dRdt(1,:) = amp*K/(K+R(2)) - dm*R(1)
 
% Dyamics for LacI protein 
dRdt(2,:) = p*R(1) - d*R(2)

end