整理文档很辛苦，赏杯茶钱您下走！
还剩 ... 页未读，继续阅读 >>
免费阅读已结束，点击下载阅读编辑剩下 ... 页
阅读已结束，您可以下载文档离线阅读编辑
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为差分表
