Я пользовался масштабированной модифицированной функицей Бесселя CBESI (Amos, D. E. A Portable Package for Bessel Functions of a Complex Argument and Nonnegative Order / D. E. Amos // ACM Transactions on Mathematical Software. — 1986, September. — Vol. 12, № 3. — P. 265–273.) из NetLib (http://netlib.org/amos/). Компилировал в dll с помощью G95 из MinGW. Отличие от Матлаба в очень широком диапазоне аргументов < 3.4E-16 (судя по документации Матлаба, в нём используется как раз эта библиотека). Судя по статьям, эта библиотека вообще очень популярна. Кстати, ещё и работает быстро, по крайней мере у меня эта функция узким местом не была.
--------------------
|