%1-1 拉格朗日插值
%x,y:數組形勢的數據表
%xi:待計算y值的橫坐標數組
%yi:用Lagrange插值算出的y值數組
function fi=Lagrange(x,f,xi)
fi=zeros(size(xi));
np1=length(f);
for i=1:np1
    z=ones(size(xi));
    for j=1:np1
        if i~=j,
           z=z.*(xi-x(j))/(x(i)-x(j));
       end
    end
    fi=fi+z*f(i);
 end