Polynômes de Legendre

%-------------------------------------------------------
% utilisation polynôme de Legendre
% (c) JPaul Molina
%                                µ²
% sol de (1-t²)y"-2ty'[v(v+1)- ----- ]y = 0
%                              1-t²
%                      v 
%            1        d        v
%  P (t) = -------  -----[t²-1]
%   v      (2^v)v!      v
%                     dt
%--------------------------------------------------------------

echo on;
%---------------------------------------------------------
%  Polynôme de Legendre
% voir contenu fichier pour explications
%---------------------------------------------------------
echo off;
clear all;
t=-5:0.01:5; 		% domaine
nu=2;     		% nu c'est le 'v' précédent
y=legendre(nu,t);
s=char( 'r','b','g' ); 	%  définit un vecteur de couleurs
clf;
for i=1:nu+1      	%  µ varie de 0 à nu
   plot(t,real(y(i,:)),s(i)); 	% trace en couleur différente 
   str=sprintf('mu=%d',i-1);
   z=y(i,:);		% astuce pour mettre le bon texte avec la bonne courbe
   text(t(1,2), z(1,2),str);
   hold on;
end

ylabel('P(t)');
xlabel('t');
title('Polynômes de Legendre');
grid;