function Z=poly_eval_lagr(T,X,Y)

% poly_lagrange : calcul du polynôme de Lagrange en T (vectoriel)
%
% *********************************************************
%
%  Z=poly_eval_lagr(T,X,Y)
%     calcul de la valeur du polynôme p_n d'interpolation défini par le nuage de point X,Y
%     en T  (à partir des polynômes de Lagrange l_i).
%
%	variables d'entrées : 
%     * T : contient le tableau ligne des réels où on évalue p_n
%     * X : contient les centres x_i, pour 0<=i<=n (deux à deux distincts)
%     * Y : contient les valeurs  f(x_i), pour 0<=i<=n 
%
%	variables de sortie
%     * Z : contient l'image du tableau T par p_n (tableau ligne)
%
% 
% 
% ************ Fonctions auxiliaires utilisées ************
%
%     poly_lagrange
%
% *********************************************************
%





% Contrôles d'entrée
%
% nombre d'arguments
if nargin~=3
   error('nombre d''arguments de la fonction incorrect');
end
na=length(X);
nb=length(Y);
if (na~=nb)
   error('X et Y n''ont pas la même taille');
end
n=na-1;

% Corps d'algorithme
Auxi=Y(1)*poly_lagrange(0,T,X);
for i=2:n+1
	Auxi=Auxi+Y(i)*poly_lagrange(i-1,T,X);
end
Z=Auxi;
