function yi=N(x,y,xi) 　　n=length(x); 　　m=length(y);　 　　Y=zeros(n);　 %计算均差表Y 　　Y(:,1)=y'; 　　for k=1:n-1 for i=1:n-k if abs(x(i+k)-x(i))<eps error('the DATA is error!'); return; end Y(i,k+1)=(Y(i+1,k)-Y(i,k))/(x(i+k)-x(i)); end 　　end 　　yi=0; %计算牛