File: NUC1.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1

(Source file text) 

---

C	LIEFERT TABELLE FUER COBALTTHERMOMETER
C	VERSION FUER TELETYPE
	REAL NORM
	KHILF=-1
	DELTA=0.00592811
	C1=1.17669692
	C2=1.13227703
	C3=0.113227703
	C4=0.408575275
	C5=0.016343011
	WRITE(4,37) DELTA
37	FORMAT(2X,'DELTA           =',E20.8)
	WRITE(4,30)
30	FORMAT(2X,'KOEFFIZIENT UF2 =',T30,$)
	READ(4,10) UF2
10	FORMAT(E15.8)
C	WRITE(4,36) UF2
36	FORMAT(2X,'KOEFFIZIENT UF2 =',E20.8)
	WRITE(4,31)
31	FORMAT(2X,'KOEFFIZIENT UF4 =',T30,$)
	READ(4,10) UF4
C	WRITE(4,35) UF4
35	FORMAT(2X,'KOEFFIZIENT UF4 =',E20.8//)
	T=4.2
48	M=0
	P2=0.
	P4=0.
	K=1

	YZ1=1.
	Y1=1.
	SUMZ=1.
	X=DELTA/T
	Y=EXP(X)
21	YZ1=YZ1*Y
	SUMZ=SUMZ+YZ1+1./YZ1
	IF (K-5) 22,23,23
22	K=K+1
	GOTO 21
23	Z=SUMZ
1	M=M+1
	M1=M*M
	Y1=Y1*Y
	P2=P2+M1*(Y1+1./Y1)
	P4=P4+M1*M1*(Y1+1./Y1)
	IF (M-5) 1,2,2
2	P2=P2/Z
	P4=P4/Z
	W1=1.+C1*UF4-C2*UF2
	W2=(C3*UF2-C4*UF4)*P2
	W3=C5*UF4*P4
	W4=W1+W2+W3
	IF (T-1.) 4,4,6
6	NORM=W4
	WRITE(4,33) T,NORM
33	FORMAT(2X,'T =  ',F6.2,'   W(4.2K,0) =',E14.8/)
	WRITE(4,13)
13	FORMAT(2X,' INTENSITAET BEI 4.2K =',T30,$)
	READ(4,14) CONST
14	FORMAT(E15.8)
	WRITE(4,38) CONST
38	FORMAT(2X,'INTENSITAET BEI 4.2K =',E15.5//)
	WRITE(4,15)
15	FORMAT(2X,'BEGINN DER TABELLE BEI T =',T30,$)
	READ(4,14) T
	WRITE(4,34)
34	FORMAT(2X,'     T =          W(T,0)/W(4.2K,0) =       DT/DW ='/)
4	W=W4*CONST/NORM
	IF (KHILF) 5,7,8
7	WRITE(4,11) T,W
11	FORMAT(3X,F8.5,E20.5)
	GOTO 5
8	X1=(T-TM1)/(W-WM1)	
	WRITE(4,12) T,W,X1
12	FORMAT(3X,F8.5,E20.5,E25.5)
5	TM1=T
	WM1=W
	KHILF=KHILF+1
	T=T-0.0005
	IF (T-0.0011) 46,46,48
46	STOP
	END

---