% CONSTITUTIVE EXPRESSION SIMPLE MODEL RESOLUTION:

% Parameters definition
Cr = 400; % molec.min-1
p = 2.38; % min-1
dm = 0.247; % min-1
dp = 0.0063; % min-1
mu = 0.017; % min-1
Kmax = 1.6e8; % cells

param = [Cr p dm dp mu Kmax];

%Initial conditions
n0 = [1 1 1];

% ODE Options:
options = odeset('RelTol',1e-6,'Stats','off','Abstol',1e-6);
[t,n]=ode45(@mc_simple,[0 300],n0,options,param);

% Graphs:
fig = figure('Color','w');

subplot(3,1,1)
plot(t,n(:,1),'k','LineWidth',1)
xlabel('Time (min)','FontName','Lato')
ylabel('mRNA (molec)','FontName','Lato')
title('mRNA vs Time','FontName','Lato','FontSize',15)
grid on

subplot(3,1,2)
plot(t,n(:,2),'r','LineWidth',1)
xlabel('Time (min)','FontName','Lato')
ylabel('PoI (molec)','FontName','Lato')
title('Protein of Interest vs Time','FontName','Lato','FontSize',15)
grid on

subplot(3,1,3)
plot(t,n(:,3),'LineWidth',1)
xlabel('Time (min)','FontName','Lato')
ylabel('Number of cells (molec)','FontName','Lato')
title('Cells vs Time','FontName','Lato','FontSize',15)
grid on