牛顿（Newton）插值法的Matlab实现

算法

关于Newton插值多项式的具体内容和算法见《数值计算方法》—丁丽娟，P130-134

该算法核心是差商表的建立

上图来自《数值计算方法》—丁丽娟，P133-134

程序

Matlab代码如下：

//  输入量
xi: 离散样点的横坐标值
yi: 离散样点的纵坐标值
x: 插值多项式中自变量符号

// 输出量：
y: Newton插值多项式

// An highlighted block
function p= Newton_fun(x,xi,yi)
n=length(xi);
f=zeros(n,n);

% 对差商表第一列赋值
for k