function delta (i, j)
USE Numeric_Kinds_Module
USE FileUnits_Module
IMPLICIT NONE
!=======================
REAL (dp) :: delta
INTEGER :: i, j
!======================
IF (i.eq.j) THEN
   delta = 1.d0
ELSE
   delta = 0.d0
ENDIF
RETURN
END function delta