%script pour caluler certaines solutions de la ricatti


%on nettoie le graphique
clf

% les differentes conditions initiales en t0=0
CI0=[0.8 , 0.73 , 0.72734 , 0.727338208159, 0.72733820815842170626197571437, 0.7273382081584217062619757143, 0.727338205 , 0.7273, 0.72 , 0.5 , 0 , -0.5 , -2 , -6];

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

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

for i=1:z0

   M=Euler(0,CI0(i));
   u=M(:,2);
   t=M(:,1);
   plot(t,u)
   hold on

end

% rajoutons quelques solutions

% les differentes conditions initiales en x0=-6
T6=[0 ,1 , 2 , 4 , 6 , 8 , 10 , 12];

Z6=size(T6);
z6=Z6(2);

%le calcul et le dessin des solutions issues de t0 dans T6 et x0=-6 

for i=1:z6

   M=Euler(T6(i),-6);
   u=M(:,2);
   t=M(:,1);
   plot(t,u)
   hold on

end


% deux solution supplementaires

% les differentes conditions initiales en x0=9
CI9=[3.02714021,3.027140213];

Z9=size(CI9);
z9=Z9(2);

%le calcul et le dessin des solutions issues de t0=9 et x0 dans CI9 

for i=1:z9

   M=Euler(9,CI9(i));
   u=M(:,2);
   t=M(:,1);
   plot(t,u)
   hold on

end
