%1.此程序為n個節點及函數值的n-1次Lagrange插值多項式求x處的近似值程序
%2.lagrange insert
%3.x0為n維插值節點向量，y0為n維被插函數值向量
%4.x為插值點（也可以是m維向量）
%5.y為插值點x處的函數值（近似值，也可以是與x同維的向量）
%6.本程序于2001.7.3.編寫? 
function y=lagrangen(x0,y0,x)
n=length(x0);
%m=length(x);
%for i=1:m
   z=x;%(i);
   s=0;
   for k=1:n
      L=1;
      for j=1:n
         if j~=k
            L=L.*(z-x0(j))/(x0(k)-x0(j));
         end
      end
      s=s+L*y0(k);
   end
   y=s;%y(i)=s;
%end