python + numpy + np.polyfit()(最小二乘多项式拟合曲线)
import numpy as np
x = np.arange(1994, 2004, 1)
y = np.array([67.052, 68.008, 69.803, 72.024, 73.400, 72.063, 74.669, 74.487, 74.065, 76.777])
# 10代表拟合10次多项式,可以自由更改
z1 = np.polyfit(x, y, 10) # 系数的集合
p1 = np.poly1d(z1) # 多项式