牛顿插值多项式python

以下是 Python 实现的牛顿插值多项式代码：

def newton_interpolation(x, y):
    """
    :param x: x 坐标序列
    :param y: y 坐标序列
    :return: 插值多项式
    """
    n = len(x)
    c = y.copy()
    for k in range(1, n):
        for i in range(n - 1, k - 1, -1):
            c[i] = (c[i] - c[i - 1]) / (x[i] - x[i - k])
    p = c[n - 1]
    for i in range(n - 2, -1, -1):
        p = c[i] + (x - x[i]) * p
    return p

其中，`x` 和 `y` 分别为已知的 x 坐标序列和 y 坐标序列，函数返回插值多项式。