[MATLAB] Lagrange sur Matlab

[salgos]

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