Lagrange sur Matlab

Fermé
salgos - 29 déc. 2010 à 20:21
Bonjour,

Je dois utiliser le script ci-dessous permettant de calculer les coefficients de polynôme interpolé mais je n'arrive pas à en comprendre le fonctionnement.

Quelqu'un peut-il me commenter le script ?

Merci

function [l,L] = lagranp(x,y)
N=length(x)-1;
l=0;
for m=1:N+1
P=1;
for k=1:N+1
if k~=m, P=conv(P,[1-x(k)])/(x(m)-x(k));
end
end
L(m,:) = P;
l=1+y(m)*P;
end