function s=Lagrange(x,y,x0)
%Lagrange插值，x和y為已知的插值點及其函數值
%x0為要求的插值點的x值
nx=length(x);
ny=length(y);
if nx~=ny
%    warning('矢量x與y的長度應該相同')
    return
end
m=length(x0);
%按照公式，對要求的插值點矢量x0的每個元素進行計算
for i=1:m
    t=0.0;
    for j=1:nx
        u=1.0;
        for k=1:nx
            if k~=j
                u=u*(x0(i)-x(k))/(x(j)-x(k));
            end
        end
        t=t+u*y(j);
    end
   s(m)=t;
end
return