function R=test_diff_div(T,X,Y)	

% test_diff_div : fonction test pour le calcul du polynôme interpolateur support quelconque
%
%
% *********************************************************
%
%	R=test_diff_div(T,X,Y)	
%     calcul de la valeur du polynôme d'interpolation défini par le nuage de point X,Y
%     sur le tableau T (utilisation du calcul des différences divisées et 
%     de l'algorithme d'évaluation de Horner) pour support quelconque
%
%	variables d'entrées : 
%     * T : contient le tableau rectangulaire des réels où est évalué p_n
%     * X : contient les valeurs x_i, pour 0 <=i<=n (dans l'ordre croissant)
%     * Y : contient les valeurs f^(r)(x_i), pour 0 <=i<=n, de telle sorte que :
%             si x_k est multiple avec x_{k-1}~=x_k=x_{k+1}=...=x_{k+r}~=x_{k+r+1} alors
%                y_k=f(x_k), y_{k+1}=f'(x_k), ..., y_{k+r}=f^(r)(x_k).
%
%	variables de sortie
%     * R : contient le tableau rectangulaire des images de T par p_n,
%           polynôme d'interpolation de f sur le support x_i, pour 0 <=i<=n
%
% 
% ************ Fonctions auxiliaires utilisées ************
%
%     diff_div, eval_horner
%
% *********************************************************


D=diff_div(X,Y);
R=eval_horner(T,X(1:end-1),D);

