10.3 将17至29岁的运动员每两岁一组分成7组，每组两人测量其旋转定向的能力

import numpy as np
import statsmodels.formula.api as smf
import pylab as plt

x = np.arange(17, 30, 2); a = np.loadtxt('data10_3.txt')

plt.rc('text', usetex=True); plt.rc('font', size=16)

plt.style.use( 'default')
plt.plot(x, a[0], '*', label='\(y_1\)')
plt.plot(x, a[1], 'o', label='\(y_2\)')
x = np.hstack([x, x]); d = {'y': a.flatten(), 'x': x}
re = smf.ols('y~x+I(x**2)', d).fit()
print(re.summary()); print('残差的方差：', re.mse_resid)
plt.legend(); plt.show()

print("学号后两位：04")