<div class="snippet-clipboard-content notranslate position-relative overflow-auto" data-snippet-clipboard-copy-content="&gt;&gt;&gt; import scipy.special
&gt;&gt;&gt; scipy.special.hankel1(0, 10**8)
(3.2060295340412074e-05+7.3063911655217044e-05j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**9)
(2.4687471886269185e-05-5.2104226538976152e-06j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**10)
(nan+nan*j)"><pre class="notranslate"><code class="notranslate">&gt;&gt;&gt; import scipy.special
&gt;&gt;&gt; scipy.special.hankel1(0, 10**8)
(3.2060295340412074e-05+7.3063911655217044e-05j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**9)
(2.4687471886269185e-05-5.2104226538976152e-06j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**10)
(nan+nan*j)
</code></pre></div>
<p dir="auto">I'm not sure but i think this should be done by either Mathematica or MATLAB but i got nan when use scipy.</p>
<div class="snippet-clipboard-content notranslate position-relative overflow-auto" data-snippet-clipboard-copy-content="&gt;&gt;&gt; import scipy.special
&gt;&gt;&gt; scipy.special.hankel1(0, 10**8)
(3.2060295340412074e-05+7.3063911655217044e-05j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**9)
(2.4687471886269185e-05-5.2104226538976152e-06j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**10)
(nan+nan*j)"><pre class="notranslate"><code class="notranslate">&gt;&gt;&gt; import scipy.special
&gt;&gt;&gt; scipy.special.hankel1(0, 10**8)
(3.2060295340412074e-05+7.3063911655217044e-05j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**9)
(2.4687471886269185e-05-5.2104226538976152e-06j)
&gt;&gt;&gt; scipy.special.hankel1(0, 10**10)
(nan+nan*j)
</code></pre></div>
<p dir="auto">I'm not sure but i think this should be done by either Mathematica or MATLAB but i got nan when use scipy.</p>