from scipy.stats import t
from numpy import sqrt
mean_m, mean_f = 45, 34
sd_m, sd_f = 9, 10
N = 100
se = sqrt((sd_m ** 2)/N + (sd_f ** 2)/N)
t_value = (mean_m - mean_f)/se
p = t.sf(t_value, N+N-2)
print(f'p={p}')
if p >= 0.05:
print('Мы НЕ можем отклонить нулевую гипотезу')
else:
print('Мы можем отклонить нулевую гипотезу')