为什么scipy.special.hankel1（0，10**10）是否返回“nan”？

Question

>>> import scipy.special
>>> scipy.special.hankel1(0, 10**8)
(3.2060295340412074e-05+7.3063911655217044e-05j)
>>> scipy.special.hankel1(0, 10**9)
(2.4687471886269185e-05-5.2104226538976152e-06j)
>>> scipy.special.hankel1(0, 10**10)
(nan+nan*j)

我认为这些例子可以与Mathematica或MATLAB一起使用。但是我不能从scipy那里得到正确的结果。我不知道为什么，可能是我犯了什么错误？在

Answer 1

我不认为这是一个bug，而是一个特性。在

您应该通过执行以下操作来启用错误报告功能

scipy.special.errprint(1)

然后你会注意到：

^{pr2}$