function [C,L,R]=lagrang(X,Y,x)
%C 最终插值函数系数，次数由大到小
%L 每一个y(i)对应x系数，次数由高到低，构成的矩阵
%X，Y插值的原始数据，x插值的数，r=f(x)
n=length(X);
L=zeros(n,n);
R=0;
for k=1:n
    V=1;
    for j=1:n
        if k~=j
            V=conv(V,poly(X(j)))/(X(k)-X(j));
        end
    end % 得到系数
    L(k,:)=V;    
end
C=Y*L;
for i=1:n
    R=R+C(i)*x^(n-i);
end

1