function LXXT = LagrangeTchebychev(a,b,N,f,XX); 

% LXXT = LagrangeTchebychev(a,b,N,f,XX); 
%
%Etant donné un segment [a,b], un entier strictement positif N et 
%une fonction de [a,b] dans R ou C déclarée pr�alablement 'inline', cette 
%routine calcule le polynôme de Lagrange de degré N interpolant 
%les valeurs de f aux images par t --> (a+b)/2 + t(b-a)/2 des N+1 zéros
%dans [-1,1] du polynôme de Tchebychev T_(N+1), puis évalue ce polynôme
%aux entrées de la liste XX de points de [a,b] (TP5, exercice 1, Question 3) 

K=[1:N+1]; 
X=(a+b)/2 + (b-a)/2 * cos((2*K-1)*pi/(2*(N+1))); 
Y=f(X); 
LXXT=LagrangeNewton(X,Y,XX); 

