function yi=lagranline(x,y,xi)
%x為向量，全部的插值節點；
%y為向量，插值節點處的函數值；
%xi為標量，自變量x;
%yi為xi處的函數估計值
n=length(x);m=length(y);
if n~=m
    error('the length of x and y must be equal');
    return;
end
for k=1:n-1
    if x(k)<=xi&xi<=x(k+1)
        yi=(xi-x(k+1))/(x(k)-x(k+1))*y(k)+(xi-x(k))/(x(k+1)-x(k))*y(k+1);
        return ;
    end
end