%牛顿插值.txt

来自「方便初学者！插值与拟合的matlab实现！ 可方便同学做课设！」· 文本 代码 · 共 16 行

function [C,D]=newpoly(x,y) %牛顿插值

n=length(x);

D=zeros(n);

D(:,1)=y';

for j=2:n  %计算差商矩阵

    for k=j:n

        D(k,j)=(D(k,j-1)-D(k-1,j-1))/(x(k)-x(k-j+1)); 

    end