写一个三次样条插值的代码

以下是一个使用 Python 的三次样条插值代码：

import numpy as np
from scipy.interpolate import CubicSpline

# 创建样本数据
x = np.array([0, 1, 2, 3, 4, 5])
y = np.array([0, 1, 4, 9, 16, 25])

# 利用样本数据创建三次样条插值函数
cs = CubicSpline(x, y)

# 生成新数据并进行插值计算
new_x = np.linspace(0, 5, 11)
new_y = cs(new_x)

# 输出插值结果
print(new_y)

在这个例子中，我们使用 `numpy` 构建了样本数据 `x` 和 `y`，表示了一个二次函数 $y=x^2$。然后，我们使用 `scipy.interpolate` 模块的 `CubicSpline` 创建了一个三次样条插值函数 `cs`，用来对样本数据进行插值。最后，我们用新的 `x` 值生成新的 `y` 值并输出。