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