Sympy 无法计算积分

问题描述

我想计算这个积分,但 sympy 无法解决它:

import sympy as sp
from sympy import oo,pi,exp

q = 1.602176565*(10**(-19))     # Joule - Elektonvolt 
h = 4.135667516*(10**(-15))     # eV.s - Planck Constant
c = 299792458                   # m/s - Light Speed
k = 8.617343*(10**(-5))         # eV.K^-1 - Boltzmann Constant
Ts = 6000 #kelvin - heat of sun
sin2Qsun = 1/46000 # 

E = sp.Symbol('E')
integral_result = sp.integrate(E**2/(exp(E/(k*Ts))-1),(E,1.12,oo))
print(integral_result)
Output: 
Integral(E**2/(exp(1.93408416801637*E) - 1),oo))

我曾尝试使用 evalf() 但它没有用。
如果我像这样删除参数中的 -1sp.integrate(E**2/(exp(E/(k*Ts))) 然后它有效:(我不知道为什么?)

Output: 
0.174655209356127

我该如何解决这个问题?提前致谢!

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)