利用Mathematica计算Newton多项式的分裂差异

Question

我有一个任务是在函数 f ，点数 n ，起点 x0 和delta ħ

所以有一些带有以下签名的方法：

NewtonPoly [n_，x0_，h_，f_] = ??

我想我知道如何构造多项式本身，但我无法迭代地计算分割的差异（任务的一部分）。我这样做：

divDiff = Table[0, {n}, {n}];
Do[(
ix = x0 + i *h;
divDiff[[ix]][[ix]] = f[x0 + i*h];
), {i, 0, n}];
Do[
Do[(
 end = x0 + j*h + i;
start = x0 + j*h;
divDiff[[start]][[
  end]] = (divDiff[[start + 1]][[end]] - 
    divDiff[[start]][[end - 1]])/(end - start);
), {j, 0, n - i}];
, {i, 1, n}];

不幸的是我收到了错误：“部件分配中的divDiff [[ix]]不是符号。” 在生成符号和数字时，我似乎无法绕过头......

如果有人可以提供帮助，我会很高兴。