% demo_affiche_courbe : ce script  permet d'afficher des courbes
% correspondant à des approximations d'équations différentielles.
%
% paramètres d'entrée : 
% * nombre de courbes ;
% * la borne supérieure de l'intervalle d'étude T ; 
% * pour chacune de ces courbes, X un vecteur contenant les valeurs, déja calculées,
%   aux instants 0,h,2h,...,Nh=T.
%

%nombre de points pour l'affichage défini par ngra 
ngra=1000;

T=input('entrez la borne supérieur de l''intervalle : ');
nc=input('entrez le nombre de courbes d''erreurs à tracer : ');
YY=cell(nc,1);
xx=cell(nc,1);
yy=cell(nc,1);
for k=1:nc
   YY{k}=input(['entrez le vecteur des valeurs de la courbe numero ',num2str(k),': ']);
end
for k=1:nc
   N=length(YY{k})-1; 
   if (N<ngra)
      h=T/N;
      xx{k}=0:h:T;      
      yy{k}=YY{k};
   else
      h=T/N;
      p=fix(N/ngra);
      auxiy=YY{k};
      yy{k}=[auxiy(1:p:N+1) auxiy(N+1)];
      x=0:h:T;
      xx{k}=[x(1:p:N+1) T];
   end
end
chlegend='legend(';
for k=1:nc-1
   chlegend=[chlegend,'''courbe ',num2str(k),''','];   
end
chlegend=[chlegend,'''courbe ',num2str(nc),''',0)'];   
clf;
hold on;
for k=1:nc
   plot(xx{k,1},yy{k,1},couleur7(k));
end
eval(chlegend);
title('valeurs approchées');
hold off;




