这个问题推荐使用sympy
用法如下
from sympy import *
t=symbols('t')
f=t**2*exp(-t**2)
Ans=integrate(f,(t,-oo,oo))
print(Ans)
#>>> sqrt(pi)/2
print(Ans.evalf())
#>>> 0.886226925452758