问题描述
我想计算这个积分,但 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() 但它没有用。
如果我像这样删除参数中的 -1
; sp.integrate(E**2/(exp(E/(k*Ts)))
然后它有效:(我不知道为什么?)
Output:
0.174655209356127
我该如何解决这个问题?提前致谢!
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)