C
      SUBROUTINE HANK(WR,H0,H1,HR,HD)
C
      IMPLICIT REAL*8(A-H,O-Z)
      COMPLEX*16 CPX,H0,H1,HR,HD
C
      CALL BESJN0(WR,BJ0,BN0)
      CALL BESJN1(WR,BJ1,BN1)
C
      CPX=(0.D0,1.D0)
      H0=BJ0-CPX*BN0
      H1=BJ1-CPX*BN1
      HR=H1/WR
      HD=H0-HR
C
      RETURN
      END
