- #1

- 290

- 0

## Main Question or Discussion Point

I'm new in programing. Does enybody knows why this don't work.

implicit none

integer i,N,Ntot,istag

real*8 t,U,tpin,k,eg,mstag

real*8 rho,Umstag,Urho

real*8 eaf,eaftot,lambdaminus

write(6,*) 'N,Ntot,t,U'

read (5,*) N,Ntot,t,U

write(36,*) Ntot/2+1

tpin=8.d0*datan(1.do)/dfloat(N)

rho=dfloat(Ntot)/dfloat(N)

Urho=U*rho/2.do

do 1000 istag,Ntot,2

mstag=dfloat(istag)/dfloat(N)

Umstag=U*mstag/2.do

eaftot=0.do

do 200,i=-Ntot/4+1,Ntot/4

k=tpin*dfloat(i)

ek=-2.do*t*dcos(k)

lambdaminus=-dsqrt(ek*ek+Umstag*Umstag)

lambdaminus=lambdaminus+Urho

eaftot=eaftot+lambdaminus

200 continue

eaftot=2.do*eaftot/dfloat(N)-U*(rho*rho-mstag*mstag)/4.do

write(36,990)

990 format(i6,f16.6)

1000 continue

end

implicit none

integer i,N,Ntot,istag

real*8 t,U,tpin,k,eg,mstag

real*8 rho,Umstag,Urho

real*8 eaf,eaftot,lambdaminus

write(6,*) 'N,Ntot,t,U'

read (5,*) N,Ntot,t,U

write(36,*) Ntot/2+1

tpin=8.d0*datan(1.do)/dfloat(N)

rho=dfloat(Ntot)/dfloat(N)

Urho=U*rho/2.do

do 1000 istag,Ntot,2

mstag=dfloat(istag)/dfloat(N)

Umstag=U*mstag/2.do

eaftot=0.do

do 200,i=-Ntot/4+1,Ntot/4

k=tpin*dfloat(i)

ek=-2.do*t*dcos(k)

lambdaminus=-dsqrt(ek*ek+Umstag*Umstag)

lambdaminus=lambdaminus+Urho

eaftot=eaftot+lambdaminus

200 continue

eaftot=2.do*eaftot/dfloat(N)-U*(rho*rho-mstag*mstag)/4.do

write(36,990)

990 format(i6,f16.6)

1000 continue

end