擅长:python、mysql、java
<p>我认为你正确地计算了牛顿系数。但是,这些不是单项系数,牛顿插值公式是</p>
<pre class="lang-none prettyprint-override"><code>p(x) = f(x_0) + f[x_0,x_1]*(x-x_0) + f[x_0,x_1,x_2]*(x-x_0)*(x-x_1) +
... + f[x_0,x_1,...,x_{n-1}]*(x-x_0)*(x-x_1)*...*(x-x_{n-2})
</code></pre>
<p>你必须像牛顿和霍纳方案一样构造多项式</p>
<pre><code>newton = 0
for ck,xk in zip(newton_coefficients[::-1],x_values[::-1]):
newton = newton*np.poly1d([1,-xk]) + np.poly1d([ck])
</code></pre>