Re: [問題] integrate (-inf,inf) 實際的積分範圍

看板Python作者 (迪兒~)時間5年前 (2020/01/17 15:18), 5年前編輯推噓1(100)
留言1則, 1人參與, 5年前最新討論串3/3 (看更多)
照著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://www.ptt.cc/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

01/17 17:15, 5年前 , 1F
平方根項建議使用sqrt函數,如fun_2=1/sqrt(t**2+1)
01/17 17:15, 1F
文章代碼(AID): #1U8L-oFy (Python)
文章代碼(AID): #1U8L-oFy (Python)