% 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);
"
