這個問題推薦使用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