Team:Paris/Modeling/first parameters/premiermodele

% Main function for the simulations

%% File to use to obtain :

% 1 - The ouput of the system

% 2 - A visualization of that output

% 3 - The value of a criteria describing whether there are oscillations

tmax=200;

t=[0:0.1:tmax];

Prop=0.3;

V0=[10;0;0;0;0];

V=simul(t,V0,40,20);

% This paragraph was useful to describe the FIFO process but has to be

% modified in order to take into account the appearing oscillations.

%

% Max=V(end,:); Find the final value, in case of 'no oscillations signals'

% Max2=max(V); Find the maximum

% Seuil=Prop*Max; Define a threshold

%

% tempsSeuil=tSeuil(V,Seuil,t); Find the time corresponding to the first

% point above the threshold.

[c,i]=critere(V(:,3)); % Compute the criteria estimating the quality of the oscillations

% Visualization of the output

figure;

plot(t,V');

legend('FlhDC','FliA','GFP','RFP','CFP');

xlabel('Time');

ylabel('Concentrations')

Z1=V(:,3);