Bessel function in quad or multiple precision 1

[vit25]

Hello,
I need some help with Bessel functions in Fortran. I'm interesting in non-stationary wave propagation in solids and I need to evaluate arbitrary Bessel function (also modified) of arbitrary order and for arbitrary argument in quad precision (or in multiple precision, if possible) for my computations. I didn't find any Fortran library for this. Please, could somebody help me with this problem? Thank you very much.
Regards
vit25

 

[xwb]

I thought Fortran already had Bessel functions or maybe that is C - j0 and j1. Anyway, you could try

http://jin.ece.uiuc.edu/routines/routines.html.

You'll need to modify the real precision to quad precision.

 

[vit25]

Thank you xwb for valuable help, I will try it.
vit25

 

[xwb]

On IVF, they are called besj0, besj1, besjn, besy0, besy1, besyn but it is only for real(4).