integrate.quad (lambda x : x * exp (-x ** 2/2)/sqrt (2 * pi), 0.0, 100000)가 왜 0입니까?

Question

이 함수는 0에서 inf 사이에서 양수입니다. 1000보다 큰 값을주는 것이 0이 아닌 이유는 합리적이지 않습니다.

import scipy.integrate as integrate 
from math import * 

integrate.quad(lambda x: x*exp(-x**2/2)/sqrt(2*pi), 0.0, 1000) 
Out[52]: 
(0.3989422804014328, 1.6471510195390376e-11) 
integrate.quad(lambda x: x*exp(-x**2/2)/sqrt(2*pi), 0.0, 100000) 
Out[54]: 
(0.0, 0.0) 

Answer 1

약 x보다 큼. 39, exp(-x**2/2)은 0.0을 제공합니다. 적분의 상한이 100000 일 때, 함수가 0이 아닌 적분 구간의 분율은 너무 작아서 알고리즘이 결코 그것을 보지 못합니다. 지금까지 quad까지는 기능이 동일합니다.