用python写一个关于拉格朗日线性插值算法的代码

以下是一个简单的Python代码实现拉格朗日线性插值算法的示例：

def lagrange_interpolation(x, y, x_new):
    """
    x: 已知的x值
    y: 已知的y值
    x_new: 需要估算的x值
    """
    n = len(x)
    result = 0.0
    for i in range(n):
        temp = y[i]
        for j in range(n):
            if i != j:
                temp *= (x_new - x[j]) / (x[i] - x[j])
        result += temp
    return result

使用示例：

# 已知数据
x = [1, 2, 3, 4]
y = [1, 4, 9, 16]

# 估算新的数据
x_new = 2.5
y_new = lagrange_interpolation(x, y, x_new)
print("x_new = %f, y_new = %f" % (x_new, y_new))

输出结果为：

x_new = 2.500000, y_new = 6.250000

这个示例中，我们输入了4个已知点的x和y值，然后输入一个需要估算的新的x值（2.5），最终输出了估算出的y值（6.25）。注意，这个算法只适用于线性插值，对于高阶插值可能会出现插值误差。