牛顿插值法代码解析：求解差商系数

这段程序是用来求解一个给定的数据点集x1和y1的牛顿插值多项式系数的程序。具体的求解过程是：

1. 首先，定义一个长度为n+1的向量c，用来存储牛顿插值多项式的系数；

2. 然后，从j=2开始循环，依次计算二阶、三阶、四阶……n阶的差商，直到求得n阶差商为止；

3. 在每一次循环中，对于每一个i，计算它与i-j+1之间的差商，并将结果存储在c(i)中；

4. 最终，c(i)中存储的就是牛顿插值多项式的系数。

具体来说，差商的求解过程如下：

1. 对于二阶差商，有：f[x(i),x(i-1)]= [f(x(i))-f(x(i-1))]/[x(i)-x(i-1)]，其中f(x)是函数在x点的取值；

2. 对于三阶差商，有：f[x(i),x(i-1),x(i-2)]= [f[x(i),x(i-1)]-f[x(i-1),x(i-2)]]/[x(i)-x(i-2)]；

3. 对于四阶差商，有：f[x(i),x(i-1),x(i-2),x(i-3)]= [f[x(i),x(i-1),x(i-2)]-f[x(i-1),x(i-2),x(i-3)]]/[x(i)-x(i-3)]；

4. 以此类推，直到求得n阶差商为止。