牛顿插值matlab程序

function[y0,X]=newtonX(x0,x,y)n=length(x);X=zeros(n);X(:,1)=y;y0=X(1,1);lx=1;forj=2:nfori=j:nX(i,j)=(X(i,j-1)-X(i-1,j-1))/(x(i)-x(i-j+1));endlx=lx*(x0-x(j-1));y0=y0+X(j,j)*lx;end%牛顿插值x为插值节点，y为插值节点处的值，输出X为差分表