1 SUBROUTINE FIG26 ( RN,MACH, CF ) 2C 3C********* SUBROUTINE TO COMPUTE FIGURE 4.1.5.1-26 4C 5 REAL MACH 6 DIMENSION A(9),B(9),C(9),D(9),E(9),AMACH(9) 7 DATA A/4.12963E-6,3.92725E-6,4.55853E-6,4.49735E-6,4.11442E-6, 8 1 4.51587E-6,4.61971E-6,4.53836E-6,3.86772E-6/ 9 DATA B/-1.36204E-4,-1.3037E-4,-1.48715E-4,-1.47407E-4,-1.36505E-4, 10 1 -1.47222E-4,-1.48773E-4,-1.44676E-4,-1.23287E-4/ 11 DATA C/1.7162E-3,1.65388E-3,1.85005E-3,1.83955E-3,1.72383E-3, 12 1 1.8252E-3,1.81973E-3,1.74944E-3,1.49335E-3/ 13 DATA D/-9.88935E-3,-9.59519E-3,-1.0503E-2,-1.04587E-2,-9.91294E-3, 14 1 -1.02911E-2,-1.01084E-2,-9.59421E-3,-8.22227E-3/ 15 DATA E/2.23641E-2,2.18366E-2,2.33437E-2,2.32398E-2,2.22626E-2, 16 1 2.2622E-2,2.18584E-2,2.04571E-2,1.76472E-2/ 17 DATA AMACH/0.0,0.3,0.7,0.9,1.0,1.5,2.0,2.5,3.0/ 18C 19 X=ALOG10(RN) 20 DO 1000 M=1,8 21 IF(MACH.GE.AMACH(M).AND.MACH.LT.AMACH(M+1)) GO TO 1020 22 1000 CONTINUE 23 M=9 24 1010 ASSIGN 1050 TO NT 25 GO TO 1060 26 1020 ASSIGN 1030 TO NT 27 GO TO 1060 28 1030 IF(ABS(MACH-AMACH(M)).LE.0.02)GO TO 1050 29 YM=TM 30 N= M+1 31 ASSIGN 1040 TO NT 32 L=N 33 GO TO 1070 34 1040 YN=TM 35 CF= YM+((MACH-AMACH(M))/(AMACH(N)-AMACH(M))) * (YN-YM) 36 RETURN 37 1050 CF=TM 38 RETURN 39 1060 L=M 40 1070 TM=X*(E(L)+X*(D(L)+X*(C(L)+X*(B(L)+X*A(L))))) 41 GO TO NT,(1030,1040,1050) 42 END 43