function yy=nalagr(x,y,xx)
% 用途：Lagrange 插值法數值求解
% 格式：yy=nalagr(x,y,xx), x是節點向量，y是節點上的函數值，xx是插值點(可以是多個)
%       yy返回插值
m=length(x);n=length(y);
if m~=n, error('向量x與y的長度必須一致');end
s=0;
for i=1:n
   t=ones(1,length(xx));
   for j=1:n
      if j~=i, 
         t=t.*(xx-x(j))/(x(i)-x(j));
      end
   end
   s=s+t*y(i);
end
yy=s;