program bessel_test
  use m_precision, only : wp 
  use m_bessel,    only  :  BesselJn
  implicit none
  !integer     :: order
  !complex(wp) :: arg



  ! --- Test some Bessel function evaluations
  write(*,*) "BesselJn(  0,  ( 0._wp, 0._wp) ) : ", &
       BesselJn(  0,  ( 0._wp, 0._wp) )
  
  write(*,*) "BesselJn(  2,  ( 1._wp, 0._wp) ) : ", &
       BesselJn(  2,  ( 1._wp, 0._wp) )
  
  write(*,*) "BesselJn( -2,  ( 0.5_wp,0._wp) ) : ", &
       BesselJn( -2,  ( 0.5_wp,0._wp) )
  
  write(*,*) "BesselJn( -3,  (-3._wp, 0._wp) ) : ", &
       BesselJn( -3,  (-3._wp, 0._wp) )
  
  write(*,*) "BesselJn( -3,  (-3._wp, 2._wp) ) : ", &
       BesselJn( -3,  (-3._wp, 2._wp) )
  

end program bessel_test