1 function R=test_diff_div_dist(T,X,Y)
\r
3 % test_diff_div_dist : fonction test pour le calcul du polynôme interpolateur (vectoriel)
\r
6 % *********************************************************
\r
8 % R=test_diff_div_dist(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).
\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 (deux à deux distinctes)
\r
16 % * Y : contient les valeurs f(x_i), pour 0 <=i<=n
\r
18 % variables de sortie
\r
19 % * R : contient le tableau rectangulaire des images de T par p_n,
\r
20 % polynôme d'interpolation de f sur le support x_i, pour 0 <=i<=n
\r
23 % ************ Fonctions auxiliaires utilisées ************
\r
25 % diff_div_dist, eval_horner
\r
27 % *********************************************************
\r
30 D=diff_div_dist(X,Y);
\r
31 R=eval_horner(T,X(1:end-1),D);
\r