照著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 --

※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 1.164.168.128 (臺灣) ※ 文章網址: https://webptt.com/m.aspx?n=bbs/Python/M.1579245490.A.3FC.html ※ 編輯: holycity (1.164.168.128 臺灣), 01/17/2020 15:30:46 ※ 編輯: holycity (1.164.168.128 臺灣), 01/17/2020 15:31:02