牛顿插值matlab求差商表的代码

以下是牛顿插值法求解差商表的 Matlab 代码：

function [c] = newton_diff(x,y)

% 计算差商表
n = length(x);
c = zeros(n,n);
c(:,1) = y';

for j=2:n
    for i=j:n
        c(i,j) = (c(i,j-1) - c(i-1,j-1))/(x(i)-x(i-j+1));
    end
end

其中，`x` 和 `y` 分别为已知数据点的横、纵坐标向量，`c` 为差商表。在函数中，首先计算第一列的差商，然后利用递推公式求出其余列的差商。