求牛顿插值法差商matlab

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

function [c] = newton_diff(x, y)
% x: 插值节点的横坐标
% y: 插值节点的纵坐标
% c: 差商表

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` 是差商表。函数返回的 `c` 可以用于计算牛顿插值多项式。