Bessel function の J₀(z) と J₁(z) は、次の様な級数で計算できます。引数のzは、実数です。

同様に Neumann function である N₀(z) と N₁(z) も次の様になります。

計算にあたって注意したいことは、zの値です。倍精度で計算しても、z=15が限度です。

そこで、z が大きくなったら J_n(z) と N_n(z) は、どうなるか気になるところです。これについて、書物に解がありましたので、参考までに紹介しておきます。下がその解です。