Matlab实现拉格朗日插值算法

function L=Lagrange(x,y,a)
%实现拉格朗日插值
%   [x,y] 输入的样本点
%   a 待插值点
clc;
n=length(x);
s=0;
for i=1:n
    x1=x;
    x1(i)=[];   %作用于内循环（替代if),防止相同样本点相减
    fz=1;fm=1;
    for j=1:n-1
        fz=fz*(a-x1(j));
        fm=fm*(x(i)-x1(j));
    end
    q=fz/fm;
    s=s+q*y(i);
end
L=s;

end