作者holycity (迪儿~)
看板Python
标题Re: [问题] integrate (-inf,inf) 实际的积分范围
时间Fri Jan 17 15:18:08 2020
照着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/cn.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
1F:推 papple23g: 平方根项建议使用sqrt函数,如fun_2=1/sqrt(t**2+1) 01/17 17:15