牛顿插值法（MATLAB）

1.牛顿插值函数

function [ y ] = Newton_insert( x0,y0 )
n=length(x0);//初始值的维数
y=y0;//y的初始值
for i=2:n
    for j=i:n
            y(j,i)=(y0(j,i-1)-y0(j-1,i-1))/(x0(j,1)-x0(j-i+1,1));//计算均差值
    end
    y0=y;
end
end

2.插值执行文件

x0=[-1,0,1,2,3,4,5]';//x的初始值
n=length(x0);
y0=zeros(n);
y0(1:n,1)=[-20,-12,1,15,4,21,41];//y0的初始值
y=Newton_insert(x0,y0)

3.执行结果
y =
-20.0000 0 0 0 0 0 0
-12.0000 8.0000 0 0 0 0 0
1.0000 13.0000 2.5000 0 0 0 0
15.0000 14.0000 0.5000 -0.6667 0 0 0
4.0000 -11.0000 -12.5000 -4.3333 -0.9167 0 0
21.0000 17.0000 14.0000 8.8333 3.2917 0.8417 0
41.0000 20.0000 1.5000 -4.1667 -3.2500 -1.3083 -0.3583
4. 编程过程存在问题：
①x0的初始值是列矩阵，容易忽略，造成计算时矩阵越界和不匹配