利用python求牛顿插值多项式

好的，下面是利用 Python 求解牛顿插值多项式的代码：

import numpy as np

def newton_interpolation(x, y):
    n = len(x)
    f = np.zeros([n, n])
    f[:, 0] = y

    for j in range(1, n):
        for i in range(j, n):
            f[i][j] = (f[i][j-1] - f[i-1][j-1]) / (x[i] - x[i-j])

    def p(z):
        res = f[0][0]
        for i in range(1, n):
            tmp = f[i][i]
            for j in range(i):
                tmp *= (z - x[j])
            res += tmp
        return res

    return p

其中，`x` 和 `y` 分别是已知的数据点的横纵坐标，返回一个插值函数 `p`，可以通过调用 `p(z)` 计算未知点的纵坐标。