Fonction de bessel modifiée en C ?

Bonjour à tous :),

Je suis en train de développer un code en C dans lequel je dois calculer des fonctions de Bessel modifiées. J'ai beau farfouiller un peu partout sur la toile, il me semble que les librairies classiques du C ne proposent pas cette fonction. Est ce que je me trompe ? Si non, auriez vous un lien à m'adresser pour récupérer une routine de calcul des fonctions de Bessel ?
NB : j'ai regardé dans le Numerical Recipes mais ils ne traitent que les fonctions de Bessel d'ordre n>1, or j'ai un ordre n=-1 :?

merci d'avance!!