入门级matlab程序，拉格朗日插值法

程序对照书本上的公式编写

function y0=Lagrange(x,y,x0)
n=length(x);
l=ones(1,n);
y0=0;
for j=1:n
    for k=1:n
        if j~=k
           l(j)=l(j)*(x0-x(k))/(x(j)-x(k));
        end
    end
end
for i=1:n
    y0=y0+l(i)*y(i);
end

输入的x，y为数组，用于拟合，y0为插值函数在x0点的值。