Normal distribution - cdf and ppf

from scipy import stats as st

# математическое ожидание
mu = 1000
# стандартное отклонение
sigma = 100

# задаём нормальное распределение с математическим ожиданием 1000 
# и стандартным отклонением 100
distr = st.norm(mu, sigma) 

# значение, для которого хотим найти вероятность
x = 1000

# считаем вероятность получить значение, равное x или меньше 
result = distr.cdf(x)  # обратная функция - ppf() (percent point function) позволяет по вероятности найти значение

print(result)  # выведет 0.5