program testmatek
implicit none
integer :: Nc=1000 ,k
integer,parameter :: N=2
REAL :: kx ,a0=1.0 ,t0=0.25
DOUBLE PRECISION :: pi=4*ATAN(1.)
COMPLEX , PARAMETER :: i=(0,1)
complex :: ek1 ,ek2
complex :: MATRIX_ek(N,N)
open(1,file='matek.txt')
MATRIX_ek(N,N)=(0,0)
DO k=-Nc,+Nc...