python - 如何解决溢出错误: math range error

Question

我收到溢出错误

import numpy as np

pi = np.pi

from scipy.integrate import quad
from math import exp

hbar = 1.055e-34
boltz = 1.381e-23
c = 2.998e8

def z(x):
    return (x**3)/(exp(x)-1)

B=quad(z,0,np.inf)       
A= ((boltz**4)*B)/(4*(pi**2)*(c**2)*(hbar**3))

print (A) 

它在第 11 行给了我一个溢出错误，即 return (x**3)/(exp(x)-1)

Answer 1

你达到了机器精度，Python 就崩溃了。

>>> def z(x):
...     return (x**3)/(exp(x)-1)
...
>>> z(709)
4.336616682334302e-300
>>> z(710)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "<stdin>", line 2, in z
OverflowError: math range error

只要积分到 ~700 就可以了。