file bessfin.f90 for modified Bessel function In(x) of first kind, integer order, real by Brian Geelen lang Fortran90 seealso toms/597