Computing Bessel functions of the second kind in extreme parameter regimes Original Research Article

A useful method of computing the integral order Bessel functions of the second kind Yn(x+iy) when either, the absolute value of the real part, or the imaginary part of the argument z=x+iy is small, is described. This method is based on computing the Bessel functions for extreme parameter regimes when x∼0 (or y∼0) and is useful because a number existing algorithms and methods fail to give correct results for small x or small y. The approximating equations are derived by expanding the Bessel function in Taylor series, are tested and discussed. The present work is a continuation of the previous one conducted in regard to the Bessel function of the first kind. The results of our formalism are compared to the available existing numerical methods used in Mathematica, IMSL, MATLAB, and the Amos library. Our numerical method is easy to implement, efficient, and produces reliable results. In addition, this method reduces the computation of the Bessel functions of the second complex argument to that of real argument which simplify the computation considerably.