File: FERMEN.FC of Tape: Various/ETH/fc4

(.FC,.FO Focal format converted to 'WRITE' listing) 

---

C-OS/8 FOCAL, 1972

01.01 C FERMIENERGIE UND DEREN VOLUMENABHAENGIGKEIT
01.10 A !"GITTERKONSTANTE IN A. U.   "A0;S CU=2*PI/A0
01.11 A !"S   ",S
01.12 A !"Z   ",Z
01.20 A !"ANZAHL VERSCHIEDENE BRAGGEBENEN   ",II
01.21 F I=1,II;A !"(Q(I))^2 IN C.U.   ",Q(I);A !"V(I)   ",V(I);A !"WIE OFT ALS BRAGGEBENE   ",NI(I)
01.22 F I=1,II;S Q(I)=FSQT(Q(I))*CU
01.30 A !"EF   ",EF;S KF=FSQT(2*EF);T !"KF=",%,KF

02.01 C FERMIEENERGIE DES FREIEN ELEKTRONS
02.10 S K0=FEXP(FLOG(3*PI^2*S*Z)/3)/A0;S E0=K0^2/2
02.20 T !!"EF0=",%,E0,!"KF0=",%,K0,!

03.01 C KORREKTUR DER FERMIENERGIE WEGEN BRAGGEBENEN
03.15 S EF=E0;S EB=0;S DE=0
03.20 F I=1,II;D 3.30;D 3.40;D 4;S EF=EF+DE(I);S EB=EB+DB(I)
03.25 T !!"EF= ",EF;T ",  EFB=",EF-EB
03.26 R
03.30 S DE(I)=-NI(I)*(V(I)^2/(2*Q(I)*K0))*FLOG(FABS((Q(I)+2*K0)/(Q(I)-2*K0)));T !"DEF=",%,DE(I)
03.40 S DB(I)=-NI(I)*(2*V(I)^2)/(Q(I)^2);T "  DEB=",DB(I)

04.01 C D(DEF)/D(V(I))
04.10 S EV=2*DE(I)/FABS(V(I))
04.12 S BV=-4*NI(I)*FABS(V(I))/(Q(I)^2)
04.15 T ",   DEF/DV=",%,EV,",  DEB/DV=",BV,",  DEFB/DV=",EV-BV

09.70 

---