File: NUC.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1
(Source file text)
C LIEFERT TABELLE FUER COBALTTHERMOMETER REAL NORM KHILF=-1 DELTA=0.00592811 C1=1.17669692 C2=1.13227703 C3=0.113227703 C4=0.408575275 C5=0.016343011 WRITE(3,37) DELTA 37 FORMAT(2X,'DELTA =',E20.8) WRITE(4,30) 30 FORMAT(2X,'KOEFFIZIENT UF2 =',T60,$) READ(4,10) UF2 10 FORMAT(E15.8) WRITE(3,36) UF2 36 FORMAT(2X,'KOEFFIZIENT UF2 =',E20.8) WRITE(4,31) 31 FORMAT(2X,'KOEFFIZIENT UF4 =',T60,$) READ(4,10) UF4 WRITE(3,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(3,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 =',T60,$) READ(4,14) CONST 14 FORMAT(E15.8) WRITE(3,38) CONST 38 FORMAT(2X,'INTENSITAET BEI 4.2K =',E15.5//) WRITE(4,15) 15 FORMAT(2X,'BEGINN DER TABELLE BEI T =',T60,$) READ(4,14) T WRITE(3,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(3,11) T,W 11 FORMAT(3X,F8.5,E20.5) GOTO 5 8 X1=(T-TM1)/(W-WM1) WRITE(3,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