骰子期望与方差

def solve():
    import numpy as np
    import matplotlib.pyplot as plt
    a, b = 1, 6
    x = np.arange(a, b+1)
    p = np.ones(6)/6
    E = np.sum(x*p)
    V = np.sum((x-E)**2*p)
    print(f'E={E:.4f}（公式={(a+b)/2:.4f}），Var={V:.4f}（公式={((b-a+1)**2-1)/12:.4f}）')
    plt.bar(x, p, color='steelblue', edgecolor='black')
    plt.xlabel('点数'); plt.ylabel('概率'); plt.title('骰子离散均匀分布')
    plt.xticks(x); plt.ylim(0, 0.2); plt.tight_layout(); plt.show()