x = Symbol('x')
M_x = integrate(x * 0.5 * cos(x), (x, -np.pi/2, np.pi/2)) # Мат. очікування (теоретич.)
D_x = integrate((x - M_x)**2 * 0.5 * cos( x), (x, -np.pi/2, np.pi/2)) # Дисперсія (теоретич.)
print(
f"Мат. очікування (теор): {M_x}\n" +
f"Дисперсія (теор): {D_x}"
)
#
X_sq = 0
for k in range(9):
P_k = integrate(0.5 * cos(x), (x, bins[k], bins[k + 1]))
X_sq += (256 * P_k - count[k])**2 / (256 * P_k)
print("Імовірність того, що експериментальні дані не суперечать гіпотезі:", 1-chi2.cdf(float(X_sq), 7))