Модифицированную функцию Бесселя, чем лучше посчитать? Опции

[Pathfinder]

Ищу библиотеку, которая позволяет вычислить модифицированную функцию Бесселя с экспоненциальным масштабным множителем. Пока что использую GSL из GNUWin32, но она не считает функцию с порядком выше 32-го.
Насколько я понял, в Матлабе используется вот эта библиотека: http://www.netlib.org/amos, но она на фортране ... боюсь, проблематично её будет с MSVC-шным проектом использовать.

[Pathfinder]

mdmitry
Насколько я понял, в Scientific Linux используется GSL, и как раз её я использую сейчас. Правда приходится пользоваться версией 1.8, поскольку свежие не понял, как под винду собрать...
За наводку про F2C спасибо, нашёл и попробовал. Исходник конечно не фонтан, но лучше чем ничего. Самое неприятное, что в результате сборки бесселя получается dll-ка размером 132КБ, поскольку приходится фортрановский рантайм прилинковывать. При таких размерах она скорее всего не влезет в страницу кэша и будет тормозить.

=GM=
Согласен, решение логичное, и в моём случае удобное, поскольку мне нужно считать I_0, I_n и I_2n для одного аргумента. Но смущает то, что его не используют (как мне показалось) ни в netlib.org/amos ни в GSL. Может при рекурсивном вычислении ошибка накапливается? Там ведь дробно-полиномиальный множитель будет...
А не помните, какая точность получалась в сравнении с таблицами? Мне очень желательно 12 знаков получить.