Yl. Ratis et Pf. Decordoba, A CODE TO CALCULATE (HIGH-ORDER) BESSEL-FUNCTIONS BASED ON THE CONTINUED FRACTIONS METHOD, Computer physics communications, 76(3), 1993, pp. 381-388
We have developed a fast code to calculate Bessel functions of integer
and fractional order based on the continued fractions method. This al
gorithm is specially useful in the case of Bessel functions of high or
der because it does not require any recalculation using normalization
relations.