Hi

I was wondering whether anyone knows of the existence of c\c++ source code for Bessel functions of integer order and complex argument (which I'll denote Jnu(z) for 1st kind, Ynu(z) for second kind)? I found some code here

http://perso.orange.fr/jean-pierre.moreau/c_bessel.html

I ran some tests on cbessj.cpp, and I found that the results were identical to those produced by Matlab when nu and z were both small. However, it seems to fail when nu and z are even moderately large. In the work I'm doing I need these functions to be computed accurately for nu=25 and even higher. Is anyone aware of such source code?

Thank you in advance for any comments and suggestions

catMan