%script pour caluler certaines solutions de la ricatti


%on nettoie le graphique
clf


%les parametres
a=2;
b=0.4;
c=1;
d=0.1;

%La duree
T=5

% Le pas
h=0.1

% les differentes conditions initiales en t0=0
CI0=[22 , 8 ;  18 , 7 ; 14, 6 ; 10, 5];

Z0=size(CI0);
z0=Z0(1);

%le calcul et le dessin des solutions issues de t0=0 et x0 dans CI0

for i=1:z0

   M=LVEuler(CI0(i,1),CI0(i,2),a,b,c,d,T,h);
   r=M(:,3);
   s=M(:,2);
   t=M(:,1);
   plot(t,r,t,s)
 
%on fait une pause entre chaque trace
%appuyer sur une touche pour passer a la boucle suivante
   pause
end
pause

% nettoyage
clf;

% tracer des courbes parametre (sardines(t), requins(t))

for i=1:z0

   M=LVEuler(CI0(i,1),CI0(i,2),a,b,c,d,T,h);
   r=M(:,3);
   s=M(:,2);
   t=M(:,1);
   plot(s,r)
   hold on
%on fait une pause entre chaque trace
%appuyer sur une touche pour passer a la boucle suivante
   pause
end
