@brief Subroutine taken from "Numerical Recipes" that calculates the modified Bessel function of
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=SP), | intent(in) | :: | x | |||
real(kind=SP), | intent(in) | :: | xnu | |||
real(kind=SP), | intent(out) | :: | ri | |||
real(kind=SP), | intent(out) | :: | rk | |||
real(kind=SP), | intent(out) | :: | rip | |||
real(kind=SP), | intent(out) | :: | rkp |
SUBROUTINE bessik(x,xnu,ri,rk,rip,rkp)
USE nrtype
REAL(SP), INTENT(IN) :: x,xnu
REAL(SP), INTENT(OUT) :: ri,rk,rip,rkp
END SUBROUTINE bessik