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