Polynome interplateur de Lagrange

. On considère 2N+1 points équirépartis autour de 0. Par ces 2N+1 points, il ne passe qu'un seul et unique polynome de degré 2N.

Prenons par exemple le polynome suivant:

> F:=sum(a[i]*X^i,i=0..2);

Nous cherchons à faire passer ce polynome par la famille de points suivante . Il suffit de poser les trois équations, de les résoudre, et d'exprimer la valeur du polynome interpolateur du second degré.

> eq_set := {'f(q*i) = subs(X=q*i,F)'$i=-1..1};

> solve(eq_set,{'a[i]'$i=0..2});

> F2:=subs(",F);