from math import* from sympy import* X=[1,-3,8] Y=[6,-8,-4] #X est le vecteur qui contient le x_i et Y celui qui contient les y_i ; # P sera le polynôme d’interpolation de Lagrange def Flagrange(X,Y) : n=len(X) # n+1 est le nombre de points. n sera donc le degré du polynôme P x=symbols("x") P=0. #x est l’indéterminée du polynôme P for i in range(0,n) : # On affecte à L le flottant 1. , ce qui aura pour conséquence de renvoyer # des flottants pour la suite du programme L=1. for j in range(0,n) : if i==j : continue L=L*(x-X[j])/(X[i]-X[j]) # on a effectué les différences divisées P=P+L*Y[i] ; return simplify(P)