function fi=Lagran_(x,y,xi) %根據給出的數據向量x和y，構造拉格朗日插值多項式，%并計算在xi處的值。
fi=zeros(size(xi));  % 全零陣
np1=length(y);
for i=1:np1
    z=ones(size(xi)); %全1陣
    for j=1:np1
        if i~=j
            z=z.*(xi-x(j))/(x(i)-x(j)); %插值基函數
        end
    end
    fi=fi+z*y(i);
end
return