Bonjour ,

J'ai besoin d'un coup de main

Je dois programmer une fonction qui prend en entrée deux vecteurs (un vecteur pour les abcisses des points considérés et un vecteur pour les cordonnées ) et donne en sortie le polynome interpolateur en ces points calculé par la méthode d Lagrange .

donc je doit appliquer ça a la fonction f(x)=exp(x) sur l'intervalle [-1,1] avec points d'interpolation équidistants puis tracer la fonction considérée et le polynome sur le meme graphique:et utiliser linspace pour les définir

voila mon essai
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
function poly=lagrange(X,Y,XX)
z=zeros(size(XX));
n=length(Y);
for i=1:n
    z=ones(size(XX));
    for j=1:n
        if i~=j
            z=z.*(XX-X(j)/X(i)-X(j));
        end
    end
    poly=poly+Z(i).*Y(i);
end