Я беру нормальное распределение со средним значением 8704 и стандартным отклонением 1. Таким образом, интеграция функции плотности по [-infinity, infinity] должна составлять 1. Но я получаю гораздо меньше. Я что-то пропустил? Я написал это в Sagemath.
from mpmath import *
import time
mp.dps = 200
mp.pretty = True
mu0=8704
sigma0=1
F1 = lambda x: 1/(sqrt(2*pi)*sigma0)*e^(-(x-mu0)^2/(2*sigma0^2))
J=quad(F1, [-inf,inf])
print J