SUBROUTINE kin_rho (n, toe, kinetic, nd_daf,kmax)
USE Numeric_Kinds_Module
USE FileUnits_Module
IMPLICIT NONE
LOGICAL :: debug
INTEGER :: i, j, n, index, nd_daf, kmax
REAL(dp) :: toe (0:nd_daf), kinetic (n, n)

debug=.false.

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 (debug) THEN
   WRITE(Output_Unit, * ) 'Rho-Kinetic Energy Matrix'
   CALL mxout (kinetic, n, n)
ENDIF

RETURN
END SUBROUTINE kin_rho