MATLAB程序设计：拉格朗日插值法

function lagerangeinterpolation(X,Y,x)
%X为输入的自变量矩阵，Y为输入的因变量矩阵，x为插值点
n=length(X);
l=[];
for k=1:n
    for i=1:n+1
        if i~=k
            l(i)=(x-X(i))/(X(k)-X(i));
        end
    end
    P(k)=l(k)*Y(k);
end

end