Wolfram mathematica 用Mathematica计算牛顿多项式的分差

我的任务是创建牛顿多项式，给定函数f、点数n、起点x0和δh

因此，一些具有以下签名的方法：

NewtonPoly[n_x0_x0，h_x0，f_x0]=

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

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]]不是符号。” 当符号生成时，当数字出现时，我似乎无法回过头来思考

---

如果有人能提供帮助，我将非常高兴。

您有很多语法和概念错误。你可以开始读了，谢谢，我会调查的。