</code> calcula o array <code>besselarray</code> tal que
<code>besselarray [i] = bessel_k [i + v - int(v)] (z)</code> para <code>i</code> de zero a <code>int(v)</code>.
</p>
<p><code>bessel_k</code> é definida como
</p><pre class="example">           %pi csc(%pi v) (bessel_i(-v, z) - bessel_i(v, z))
           -------------------------------------------------
                                  2
</pre>
<p>quando <em>v</em> não for inteiro.  Se <em>v</em> for um inteiro <em>n</em>,
então o limite  com <em>