Team:Paris/Modeling/plot

% Main function for the simulations

clear, close all;

%% 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=[1;0;0;0;0;0;0;0;0];%Initial state

V=simul(t,V0);%Numerical simulation of the system

% Visualization of the output

Z=[V(:,1),V(:,5),V(:,6),V(:,9)];

figure;

plot(t,V');

legend('FlhDC','FliA','FliL','FlgA','FlhB','AHL interieur','AHL ext','mTetR','TetR');

%legend('FlhDC','FlhB','AHL int','TetR');

% xlabel('Time');

% ylabel('Concentrations')