-
double precision
BESSI(N, X)
from http://perso.orange.fr/jean-pierre.moreau/Fortran/tbessi_f90.txt
This subroutine calculates the first kind modified Bessel function
of integer order N, for any REAL X. We use here the classical
recursion formula, when X > N. For X < N, the Miller's algorithm
is used to avoid overflows.
REFERENCE:
C.W.CLENSHAW, CHEBYSHEV SERIES FOR MATHEMATICAL FUNCTIONS,
MATHEMATICAL TABLES, VOL.5, 1962.
-
integer intent(in) :: N
-
double precision :: X
double precision
BESSI0(X)
double precision :: X
double precision
BESSI1(X)
double precision :: X
real(dl)
bessj(x, n)
real(dl) :: x
integer :: n
real(dl)
bessj0(x)
real(dl) :: x