Re: [问题] integrate (-inf,inf) 实际的积分范围

楼主: holycity (迪儿~)   2020-01-17 15:18:08
照着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
作者: papple23g (逆道者)   2020-01-17 17:15:00
平方根项建议使用sqrt函数,如fun_2=1/sqrt(t**2+1)

Links booklink

Contact Us: admin [ a t ] ucptt.com