S***y
发帖数: 186 | 1 SUBROUTINE SB(LMAX,Z,JL)
INTEGER::LMAX
COMPLEX::JL(0:LMAX),Z
INTEGER::L
COMPLEX(KIND=8)::DJL(0:LMAX),NORM
DJL(LMAX)=(0.0D0,0.0D0)
DJL(LMAX-1)=(1.0D0,0.0D0)
DO L=LMAX-2,0,-1
DJL(L)=(L+L+3)*DJL(L+1)/Z-DJL(L+2)
END DO
NORM=DJL(0)
DJL(0)=SIN(Z)/Z
NORM=DJL(0)/NORM
DO L=1,LMAX
DJL(L)=DJL(L)*NORM
END DO
JL=DJL
END SUBROUTINE SB
This subroutine tabulates Spherical Bessel Function in JL(0 |
|
