% repressilator_plot.m - run the repressilator model
% RMM, 24 Jan 07

sol = ode45(@repressilator, [0 20000], [1 0 0 200 0 0]);
plot(sol.x/60, sol.y(4,:), sol.x/60, sol.y(5,:), sol.x/60, sol.y(6,:));

title('Elowitz and Leibler Repressilator model (courtesy MIT SBC)')
xlabel('time (min)');
ylabel('Protein concentration (# molecules)');
legend('cI', 'lacI', 'tetR');
