newton插值法python

Newton插值法是一种多项式插值方法，可以用于在给定数据点的情况下，构造一个经过这些点的多项式函数。下面是Python实现Newton插值法的代码：

def newton_interpolation(x, y):
    n = len(x)
    a = y.copy()
    for j in range(1, n):
        for i in range(n-1, j-1, -1):
            a[i] = (a[i] - a[i-1]) / (x[i] - x[i-j])
    def f(t):
        result = a[-1]
        for i in range(n-2, -1, -1):
            result = result * (t - x[i]) + a[i]
        return result
    return f

其中，x和y分别是数据点的横坐标和纵坐标，返回值是一个函数，可以用于计算插值多项式在任意点的函数值。