两个正态总体N（μ₁,σ₁²），N（μ₂,σ₂²）的均值差μ₁-μ₂的置信区间(σ₁²=σ₂²且均未知时)

#Python代码实现如下：
ltx1=[221,244,243,288,233,220,210,258,245]
ltx2=[268,213,188,189,217,207]
def compute(ltx):
    le=len(ltx)
    sum=0
    S=0
    for i in range(0,le):
        sum+=ltx[i]
    avg=sum/le
    print('avg:',avg)
    for i in range(0,le):
        S+=(ltx[i]-avg)*(ltx[i]-avg)
    s=S/(le-1)
    print('s^2:',s)
    return s,avg
s1,x1=compute(ltx1)
s2,x2=compute(ltx2)
v1=(s1*s1/9+s2*s2/6)**2
v2=(s1*s1/9)**2/8+(s2*s2/6)**2/5
v=v1/v2
print("自由度v:",v)
#查表t₀.₂₅[7]=2.3646
se1=(x1-x2)+2.3646*v1**0.25
se2=(x1-x2)+2.3646*v1**0.25
print('μ₁-μ₂的置信区间为',(se1,se2))