1

我想计算这个积分，但 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), (E, 1.12, oo))

我曾尝试使用 evalf() 但它没有用。
如果我像这样删除-1参数；sp.integrate(E**2/(exp(E/(k*Ts)))然后它起作用了：（我不知道为什么？）

Output: 
0.174655209356127

我怎么解决这个问题？提前致谢！

4

0 回答 0