照著p大的做法做有順利輸出一樣的結果，
也有試著修改待積函數f去得到一些合理的輸出
不過測試待積函數時遇到一點問題，
測試code如下:
from sympy import *
t = symbols('t')
fun_1 = (t**2)**(-0.5)
fun_2 = (t**2 + 1)**(-0.5)
int_1 = integrate(fun_1,(t,0,100000000))
int_2 = integrate(fun_2,(t,0,100000000))
print(int_1)
print(int_2)
兩組積分式子只差在t**2後面是否有加一個常數
fun_1沒加常數的情況下可以正常輸出沒問題
fun_2加了常數的版本輸出會跑出
_CoeffExpValueError: expr not of form a*x**b: 0.564189583547756
猜測是否可能是格式不符合sympy的要求，不過也檢查不出個所以然來...
想請教一下p大或知道sympy的高手我是否哪邊疏忽了QQ
※ 引述《papple23g (逆道者)》之銘言：
: 這個問題推薦使用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