testbesselkv Subroutine

  SUBROUTINE testbesselkv()
    IMPLICIT NONE
    REAL(rp) :: v
    REAL(rp), DIMENSION(:), ALLOCATABLE :: x
    REAL(rp), DIMENSION(:), ALLOCATABLE :: R
    INTEGER :: nx
    INTEGER :: ii
    REAL(4) :: xnu,ri,rk,rip,rkp

    nx = 1000
    v = 1.0_rp/3.0_rp

    ALLOCATE(x(nx))
    ALLOCATE(R(nx))

    do ii=1_idef,nx
       x(ii) = REAL(ii,rp)*0.01_rp
    end do

    do ii=1_idef,nx
       call bessik(REAL(x(ii),4),REAL(v,4),ri,rk,rip,rkp)
       write(6,'(F25.16)') rk
    end do

    DEALLOCATE(x)
    DEALLOCATE(R)
  END SUBROUTINE testbesselkv