% This program calls a function which implements the basic model of the repressilator,

% Close all windows and clear all variables
clear all;
close all;

% Set the simulation to run over 3 hours (3*60*60 seconds)
tspan = [0 12*60*60];

% Set the initial conditions of the 5 species in molecules per cell
% [mRNA_cI, mRNA_lacI, mRNA_tetR, protein_cI, protein_lacI, protein_tetR] 
IC = [0, 0, 0, 0, 100, 0];

% Run the simulation!!
[t result] = ode23s(@repressilator, tspan, IC, []);

% Plot the concentrations of the 2 protein species
figure(2);
plot(t/60, result(:,4), 'b-', t/60, result(:,5), 'r-', t/60, result(:,6), 'g-');
title('Oscillation of the 3 repressor proteins');
legend('cI', 'lacI', 'tetR');
xlabel('Time (min)');
ylabel('Proteins per cell');