% script demo_trace_poly_gauss : tracé des polynômes de gauss P_0, P_1,...,P_n 
%
%
%    pour l'une des quatre familes :  
%             Gauss-Legendre,
%             Gauss-Tchebytchev,
%             Gauss-Hermite,
%             Gauss-Laguerre.
%
%
%
% ************ Fonctions auxiliaires utilisées ************
%  
%   calcul_poly_gauss 
%
% *********************************************************


n=input('entrez l''entier naturel n : ');
N=input('entrez le nombre de points pour le graphique : ');
disp('choix de la famille : ');
disp('c=1 : Gauss-Legendre');
disp('c=2 : Gauss-Tchebytchev');
disp('c=3 : Gauss-Hermite');
disp('c=4 : Gauss-Laguerre');
c=input('entrez c : ');
if (c==3) | (c==4) 
   A=input('entrez la borne supérieure de l''intervalle d''etude : ');
else
   A=1;
end
choix=input('voulez vous voir les légendes ? entrez o ou n : ','s');
if (c==4)
   h=A/(N);
   x=0:h:A;
else
   h=2*A/(N);
   x=-A:h:A;   
end
P=calcul_poly_gauss(c,n);
clf;
hold on;
if strcmp(choix,'o')
   touslesnoms=[];
   toutesleslegendes=[];
   Q=5;
   Qtestt=0;
   for k=0:n
      nom =num2str(k);
      long=length(nom);
      Qtestt=Qtestt+1;
      touslesnoms=[touslesnoms,' ',nom];
      legende=[nom,blanks(Q-long+1)];
      toutesleslegendes=[toutesleslegendes;legende];
   end
end
for k=0:n
   cou=couleur7(k);
   plot(x,polyval(P(k+1,1:k+1),x),cou);      
end
if strcmp(choix,'o')
   legend(strvcat(toutesleslegendes),0);
end   
hold off;
