Функции Бесселя целого порядка от комплексного аргумента

[SciSimba]

Не могу найти вычислительноудобного способа представления функций Бесселя от комплексного аргумента.
Если для функции первого рода J (переходящую для чисто мнимого аргумента в I) можно воспользоваться хорошим интегральным представлением, то для функции второго рода (неограниченной в нуле) ничего хорошего мною найдено не было.

Если кому-нибудь такая задача уже встречалась, особенно в вычислительном плане (возможно вы сталкивались с пакетами на Fortran или C/C++), если можно, поделитесь информацией.

Спасибо заранее.

P.S. MathCAD, MatLab, Mathematica и т.п. не предлагать. Ищутся книжки, статьи пакеты исходных кодов (открытых).