SUBROUTINE Kin_Radial_s(n, toe, kinetic,nd_daf,kmax)
USE Numeric_Kinds_Module
!     CALL kin_radial(nmax+1, toeplitz, kinetic, rvals, kmax)
!-----------------------------------------------------------------------
!this subroutine determines the radial-kinetic energy matrix
!-----------------------------------------------------------------------
IMPLICIT NONE
LOGICAL, PARAMETER :: idbug=.false.
INTEGER  :: i, j, n, index, kmax,nd_daf
REAL(dp) :: toe(0:nd_daf), kinetic(n,n)

DO i=1,n
   DO j=1,n
      index=iabs(i-j)
      IF(index.lt.kmax)THEN
         kinetic(i,j)=kinetic(i,j)+toe(index)
      ENDIF
   ENDDO
ENDDO

IF(idbug)THEN
   WRITE(*,*)'Radial-Kinetic Energy Matrix'
   CALL matrix_out(kinetic, n, n, "Radial-Kinetic.txt")
ENDIF

RETURN
END SUBROUTINE Kin_Radial_s