| 1 $JOB 99015,RATAPP RAT00000 $EXECUTE IBJOB IBJOB VERSION 5 HAS CONTROL. | $IBJOB DLOGIC,MAP,NODECK RAT00010 | $IBFTC INPUT RAT00020 1 12/01/64 PAGE 1 | INPUT - EFN SOURCE STATEMENT - IFN(S) - | 0| COMMON/LAES/U,N1,N2,M1 RAT00030 | COMMON/RECO/ITEST(80),N4 RAT00040 0COMMON X(1000),Y(1000),W(1000),JX,ZERO(25),FZERO(25), RAT00050 | 1COEF(25),ACC,BMPAR,CPAR,M,N,INDER,INDX0,INDXK,INDPOT,INDMN, RAT00060 | 2INDCHA,INDACC,MN2,MPTEST,YMEAN(5),ITERA,MNTEST,AM(25),AMPAR RAT00070 | 0DIMENSION LT(15),ARG(1000),FUNC(1000),IOPW(110),IAD1W(110), RAT00080 | 1IAD2W(110),IAD3W(110),IOPX(110),IAD1X(110),IAD2X(110),IAD3X(110), RAT00090 2IOPY(110),IAD1Y(110),IAD2Y(110),IAD3Y(110),CSX(137),CSY(137), RAT00100 | 3CSW(137),XP(1000) RAT00110 | 0DATA LT/5HTITLE,5HTABLE,6HDEGREE,6HACCURA,6HZEROES,6HMPARAM, RAT00120 | 16HCPARAM,6HINDICA,5HSTART,5HPRINT,6HWEIGHT,6HARGUME,6HFUNCTI, RAT00130 | 25HERROR,6HAPPROX/ RAT00140 5 N4=0 RAT00150 | 1 CALL RECORD (U,N1,N2,M1) RAT00160 3 | IF(M1.NE.64) GO TO 1 RAT00170 | 7 DO 2 I=1,13 RAT00180 | IF(LT(I).EQ.N1) GO TO 3 RAT00190 2 CONTINUE RAT00200 | IFEL=1 RAT00210 | 6 GO TO (401,402,403,404,405,406,407,408,409,410, RAT00220 | 1411,411,411,411,411,411,411,411,411,420,421,1,423,424),IFEL RAT00230 | 401 WRITE(6,441) RAT00240 20 441 FORMAT(16H0WRONG DIRECTIVE) RAT00250 | GO TO 1 RAT00260 | 402 WRITE(6,442) RAT00270 22 | 442 FORMAT(20H0TABLE WITHOUT TITLE) RAT00280 | GO TO 1 RAT00290 403 WRITE(6,443)JY,I RAT00300 24 | 443 FORMAT(26H0TABLE OR PRINT DATA WRONG,2I6) RAT00310 | GO TO 8 RAT00320 | 404 WRITE(6,444) RAT00330 26 | 444 FORMAT(19H0DEGREES DATA WRONG) RAT00340 GO TO 8 RAT00350 | 405 WRITE(6,445) RAT00360 28 | 445 FORMAT(20H0ACCURACY DATA WRONG) RAT00370 | GO TO 8 RAT00380 | 406 WRITE(6,446) RAT00390 30 446 FORMAT(18H0ZEROES DATA WRONG) RAT00400 | GO TO 8 RAT00410 | 407 WRITE(6,447) RAT00420 32 | 447 FORMAT(22H0MPARAMETER DATA WRONG) RAT00430 | GO TO 8 RAT00440 408 WRITE(6,448) RAT00450 34 | 448 FORMAT(22H0CPARAMETER DATA WRONG) RAT00460 | GO TO 8 RAT00470 | 409 WRITE(6,449) RAT00480 36 | 449 FORMAT(21H0INDICATOR DATA WRONG) RAT00490 GO TO 8 RAT00500 | 410 GO TO (490,491,492),JJJ RAT00510 | 490 WRITE(6,493) RAT00520 39 | 493 FORMAT(22H0WEIGHT CONSTANT WRONG) RAT00530 | GO TO 8 RAT00540 491 WRITE(6,494) RAT00550 41 | 494 FORMAT(24H0ARGUMENT CONSTANT WRONG) RAT00560 | GO TO 8 RAT00570 | 1 12/01/64 PAGE 2 | INPUT - EFN SOURCE STATEMENT - IFN(S) - 0 492 WRITE(6,495) RAT00580 43 | 495 FORMAT(24H0FUNCTION CONSTANT WRONG) RAT00590 | GO TO 8 RAT00600 | 411 IFEL=IFEL-10 RAT00610 | GO TO (480,481,482),JJJ RAT00620 480 WRITE(6,483)IFEL RAT00630 47 | 483 FORMAT(6H0FAULT,I2,18H IN WEIGHT FORMULA) RAT00640 | GO TO 1 RAT00650 | 481 WRITE(6,484)IFEL RAT00660 49 | 484 FORMAT(6H0FAULT,I2,20H IN ARGUMENT FORMULA) RAT00670 GO TO 1 RAT00680 | 482 WRITE(6,485)IFEL RAT00690 51 | 485 FORMAT(6H0FAULT,I2,20H IN FUNCTION FORMULA) RAT00700 | GO TO 1 RAT00710 | 420 WRITE(6,460) RAT00720 53 460 FORMAT(25H0ARGUMENTS NOT MONOTONOUS) RAT00730 | GO TO 1 RAT00740 | 421 WRITE(6,461) RAT00750 55 | 461 FORMAT(36H0START BEFORE TABLE OR FORMULA-FAULT) RAT00760 | GO TO 1 RAT00770 423 WRITE(6,463) RAT00780 57 | 463 FORMAT(19H0FAULTS AFTER PRINT) RAT00790 | GO TO 8 RAT00800 | 424 WRITE(6,464) RAT00810 59 | 464 FORMAT(17H0NO PRINT ALLOWED) RAT00820 GO TO 1 RAT00830 | 8 IF(M1.EQ.64)GO TO 7 RAT00840 | GO TO 1 RAT00850 | 3 IF(I.NE.10)IPRES=0 RAT00860 | GO TO (10,20,30,40,50,40,40,80,90,100,110,120,130),I RAT00870 C RAT00880 | C TITLE RAT00890 | 10 READ(5,11)ITEST RAT00900 69 | 11 FORMAT (80A1) RAT00910 | WRITE (6,12)ITEST RAT00920 71 12 FORMAT (1H1,80A1) RAT00930 | MDEG=0 RAT00940 | NDEG=0 RAT00950 | ACCIN=0. RAT00960 | AMPAR=2. RAT00970 CPAR=1.05 RAT00980 | INDER=0 RAT00990 | INDX0=0 RAT01000 | INDXK=0 RAT01010 | INDPOT=0 RAT01020 INDCHA=0 RAT01030 | INDACC=0 RAT01040 | K2X=0 RAT01050 | K2Y=0 RAT01060 | K2W=0 RAT01070 IPWEI=0 RAT01080 | JX=0 RAT01090 | IPSTAB=0 RAT01100 | IPZER=0 RAT01110 | GO TO 5 RAT01120 C RAT01130 | 1 12/01/64 PAGE 3 | INPUT - EFN SOURCE STATEMENT - IFN(S) - | 0C TABLE RAT01140 | 20 IFEL=2 RAT01150 IF(JX.NE.0) GO TO 6 RAT01160 | N4=0 RAT01170 | IFEL=3 RAT01180 | I=0 RAT01190 | IF(ITEST(10).EQ.1891830832) GO TO 25 RAT01200 22 CALL REC(ARG,JY,K) RAT01210 100 | JY=JY+1 RAT01220 | GO TO (29,6),K RAT01230 | 29 I=I+1 RAT01240 | IF(I.GT.1000) GO TO 6 RAT01250 CALL RECORD (FUNC(I),N1,N2,M1) RAT01260 109 | 24 IF((M1.NE.48.AND.M1.NE.64).OR.N2.NE.48)GO TO 6 RAT01270 | IF(M1.NE.64) GO TO 29 RAT01280 | IF(JY.NE.I) GO TO 6 RAT01290 | GO TO 220 RAT01300 25 JY=0 RAT01310 | 23 JY=JY+1 RAT01320 | IF(JY.GT.1000)GO TO 6 RAT01330 | CALL RECORD (ARG(JY),N1,N2,M1) RAT01340 128 | IF((M1.NE.48.AND.M1.NE.64).OR.N2.NE.48)GO TO 6 RAT01350 I=JY RAT01360 | IF(M1.EQ.64) GO TO 220 RAT01370 | CALL RECORD (FUNC(I),N1,N2,M1) RAT01380 139 | IF(M1.NE.48.OR.N2.NE.48) GO TO 6 RAT01390 | GO TO 23 RAT01400 220 IF(I.LT.7)GO TO 6 RAT01410 | JX=I-1 RAT01420 | GO TO 7 RAT01430 | C RAT01440 | C DEGREES RAT01450 30 IFEL=4 RAT01460 | CALL RECORD (V,N1,N2,M1) RAT01470 151 | IF(N2.NE.48.OR.M1.NE.48) GO TO 6 RAT01480 | 239 CALL RECORD (U,N1,N2,M1) RAT01490 157 | IF(N2.NE.48.OR.M1.NE.48) GO TO 6 RAT01500 39 CALL RECORD (UU,N1,N2,M1) RAT01510 163 | IF(M1.NE.64) GO TO 6 RAT01520 | GO TO (6,6,31,32,33,34,35,7),I RAT01530 | 31 MDEG=V RAT01540 | NDEG=U RAT01550 GO TO 7 RAT01560 | 32 ACCIN=U RAT01570 | GO TO 7 RAT01580 | 33 IPZER=30*MDEG+NDEG+1 RAT01590 | GO TO 7 RAT01600 34 AMPAR=U RAT01610 | GO TO 7 RAT01620 | 35 CPAR=U RAT01630 | GO TO 7 RAT01640 | C RAT01650 C ACCURACY,MPARAMETER,CPARAMETER RAT01660 | 40 IFEL=I+1 RAT01670 | GO TO 239 RAT01680 | C RAT01690 | 1 12/01/64 PAGE 4 INPUT - EFN SOURCE STATEMENT - IFN(S) - | 0C ZEROES RAT01700 | 50 J=MDEG+NDEG+2 RAT01710 | IFEL=6 RAT01720 | DO 51 K=2,J RAT01730 CALL RECORD (ZERO(K),N1,N2,M1) RAT01740 189 | IF(N2.NE.48.OR.M1.NE.48) GO TO 6 RAT01750 | 51 CONTINUE RAT01760 | GO TO 39 RAT01770 | C RAT01780 C INDICATOR RAT01790 | 80 READ(5,81)INDER,INDX0,INDXK,INDPOT,INDMN,INDCHA,INDACC RAT01800 198 | 81 FORMAT(7I1) RAT01810 | N4=21 RAT01820 | IFEL=I+1 RAT01830 GO TO 39 RAT01840 | C RAT01850 | C START RAT01860 | 90 K=1 RAT01870 | 299 IFEL=21 RAT01880 IF(JX.EQ.0.OR.IPSTAB.NE.0)GO TO 6 RAT01890 | IF(K2X.EQ.0)GO TO 91 RAT01900 | DO 92 I=1,JX RAT01910 | CSX(1)=ARG(I) RAT01920 | CSX(2)=FUNC(I) RAT01930 DO 93 J=1,K2X RAT01940 | 93 CALL OPERA(IOPX(J),IAD1X(J),IAD2X(J),IAD3X(J),CSX) RAT01950 229 | J=IAD3X(K2X) RAT01960 | 92 X(I)=CSX(J) RAT01970 | GO TO 94 RAT01980 91 DO 95 I=1,JX RAT01990 | 95 X(I)=ARG(I) RAT02000 | 94 IF(K2Y.EQ.0) GO TO 96 RAT02010 | DO 97 I=1,JX RAT02020 | CSY(1)=ARG(I) RAT02030 CSY(2)=FUNC(I) RAT02040 | DO 98 J=1,K2Y RAT02050 | 98 CALL OPERA(IOPY(J),IAD1Y(J),IAD2Y(J),IAD3Y(J),CSY) RAT02060 263 | J=IAD3Y(K2Y) RAT02070 | 97 Y(I)=CSY(J) RAT02080 GO TO 99 RAT02090 | 96 DO 290 I=1,JX RAT02100 | 290 Y(I)=FUNC(I) RAT02110 | 99 IF(K2W.EQ.0)GO TO 291 RAT02120 | DO 292 I=1,JX RAT02130 CSW(1)=ARG(I) RAT02140 | CSW(2)=FUNC(I) RAT02150 | DO 293 J=1,K2W RAT02160 | 293 CALL OPERA(IOPW(J),IAD1W(J),IAD2W(J),IAD3W(J),CSW) RAT02170 297 | J=IAD3W(K2W) RAT02180 292 W(I)=CSW(J) RAT02190 | GO TO 294 RAT02200 | 291 IF(IPWEI.EQ.1)GO TO 294 RAT02210 | DO 295 I=1,JX RAT02220 | 295 W(I)=1. RAT02230 294 Z1=X(1)-X(2) RAT02240 | IFEL=20 RAT02250 | 1 12/01/64 PAGE 5 | INPUT - EFN SOURCE STATEMENT - IFN(S) - | 0 DO 298 I=3,JX RAT02260 IF(Z1*(X(I-1)-X(I)).LE.0.)GO TO 6 RAT02270 | 298 CONTINUE RAT02280 | GO TO (297,309),K RAT02290 | 297 ACC=ACCIN RAT02300 | CALL CHECK RAT02310 332 M=MDEG RAT02320 | N=NDEG RAT02330 | I=30*M+N+1 RAT02340 | L=1 RAT02350 | IF(IPZER.EQ.I) L=2 RAT02360 CALL COMPUT(L,I) RAT02370 340 | IF(I.NE.3)IPRES=1 RAT02380 | GO TO 5 RAT02390 | C RAT02400 | C PRINT RAT02410 100 CALL RECORD(U,N1,N2,M1) RAT02420 346 | DO 101 I1=1,3 RAT02430 | IF(N1.EQ.LT(I1+12)) GO TO 102 RAT02440 | 101 CONTINUE RAT02450 | IFEL=23 RAT02460 GO TO 6 RAT02470 | 102 IFEL=24 RAT02480 | IF(IPRES.EQ.1)GO TO 309 RAT02490 | IF(I1.NE.1)GO TO 6 RAT02500 | K=2 RAT02510 GO TO 299 RAT02520 | 309 CALL REC(XP,JY,K) RAT02530 369 | GO TO (103,6),K RAT02540 | 103 SPAR1=ZERO(2) RAT02550 | SPAR 2=FZERO(2) RAT02560 GO TO (300,301,302),I1 RAT02570 | 300 WRITE(6,303) RAT02580 373 | 303 FORMAT(9H-FUNCTION) RAT02590 | GO TO 306 RAT02600 | 301 WRITE(6,304) RAT02610 375 304 FORMAT(6H-ERROR) RAT02620 | GO TO 306 RAT02630 | 302 WRITE(6,305) RAT02640 377 | 305 FORMAT(13H-APPROXIMATOR) RAT02650 | 306 K=1 RAT02660 109 L=MIN0(K+6,JY) RAT02670 | DO 107K1=K,L RAT02680 | K2=K1-K RAT02690 | ZERO(2)=XP(K1) RAT02700 | GO TO (104,105,106),I1 RAT02710 104 CALL INTERP(M1,2) RAT02720 388 | AM(K2+1)=FZERO(2) RAT02730 | GO TO 107 RAT02740 | 105 CALL INTERP(M1,2) RAT02750 393 | AM(K2+1)=FZERO(2)-FAPP(ZERO(2)) RAT02760 395 GO TO 107 RAT02770 | 106 AM(K2+1)=FAPP(ZERO(2)) RAT02780 399 | 107 CONTINUE RAT02790 | K3=L-K+1 RAT02800 | WRITE(6,108)XP(K),(AM(K2),K2=1,K3) RAT02810 404 1 12/01/64 PAGE 6 | INPUT - EFN SOURCE STATEMENT - IFN(S) - | 0 108 FORMAT(1H ,1PE12.5,2X,7E15.7) RAT02820 | IF(L.EQ.JY) GO TO 307 RAT02830 | K=L+1 RAT02840 GO TO 109 RAT02850 | 307 ZERO(2)=SPAR1 RAT02860 | FZERO(2)=SPAR2 RAT02870 | GO TO 1 RAT02880 | C RAT02890 C WEIGHT RAT02900 | 110 N4=10 RAT02910 | J=J+1 RAT02920 | IF(ITEST(10).EQ.1891830832) GO TO 119 RAT02930 | 112 CALL STABLE (3467,K2W,IFEL,IOPW,IAD1W,IAD2W,IAD3W) RAT02940 425 JJJ=1 RAT02950 | 118 M1=0 RAT02960 | IPSTAB=IPSTAB+1 RAT02970 | IF(IFEL.NE.0) GO TO 6 RAT02980 | IFEL=10 RAT02990 I=2 RAT03000 | 113 I=I+1 RAT03010 | GO TO (114,115,116),JJJ RAT03020 | 114 CALL RECORD (CSW(I),N1,N2,M1) RAT03030 438 | GO TO 117 RAT03040 115 CALL RECORD (CSX(I),N1,N2,M1) RAT03050 442 | GO TO 117 RAT03060 | 116 CALL RECORD(CSY(I),N1,N2,M1) RAT03070 446 | 117 IF(N2.NE.48.OR.(M1.NE.48.AND.M1.NE.64).OR.I.GT.28)GO TO 6 RAT03080 | IF(M1.EQ.48) GO TO 113 RAT03090 GO TO 311 RAT03100 | 119 IFEL=10 RAT03110 | IF(JX.EQ.0) GO TO 6 RAT03120 | I=0 RAT03130 | 310 I=I+1 RAT03140 CALL RECORD(W(I),N1,N2,M1) RAT03150 463 | IF(N2.EQ.48.AND.M1.EQ.48)GO TO 310 RAT03160 | IF(M1.NE.64.OR.N2.NE.48) GO TO 6 RAT03170 | IF(I-1.NE.JX)GO TO 6 RAT03180 | IPWEI=1 RAT03190 311 IPSTAB=IPSTAB-1 RAT03200 | GO TO 7 RAT03210 | C RAT03220 | C ARGUMENT RAT03230 | 120 CALL STABLE (1099,K2X,IFEL,IOPX,IAD1X,IAD2X,IAD3X) RAT03240 477 JJJ=2 RAT03250 | GO TO 118 RAT03260 | C RAT03270 | C FUNCTION RAT03280 | 130 CALL STABLE (1419,K2Y,IFEL,IOPY,IAD1Y,IAD2Y,IAD3Y) RAT03290 481 JJJ=3 RAT03300 | GO TO 118 RAT03310 | END RAT03320 | 1 12/01/64 PAGE 7 | 0| -| $IBFTC RECRD1 RAT03330 | 1 12/01/64 PAGE 8 | RECRD1 - EFN SOURCE STATEMENT - IFN(S) - 0| SUBROUTINE REC(ARG,I,K) RAT03340 | COMMON/LAES/U,N1,N2,M1 RAT03350 | DIMENSION ARG(1),ISIG(4) RAT03360 | DATA ISIG/48,60,28,44/ RAT03370 I=0 RAT03380 | 22 MG=48 RAT03390 | M=48 RAT03400 | 28 CALL RECORD(U,N1,N2,M1) RAT03410 5 | DO 1 L=1,4 RAT03420 IF(ISIG(L).EQ.M1) GO TO 2 RAT03430 | 1 CONTINUE RAT03440 | GO TO 6 RAT03450 | 2 IF(M.NE.48.AND.N2.NE.48) GO TO 6 RAT03460 | IF(M.NE.48) GO TO 25 RAT03470 M=N2 RAT03480 | 25 IF((I.EQ.0.AND.M.NE.48).OR.(I.EQ.1.AND.M.EQ.28))GO TO 6 RAT03490 | IF(MG.NE.60) GO TO 229 RAT03500 | IF(M.NE.28)GO TO 6 RAT03510 | DELT=ARG(I) RAT03520 BEG=ARG(I-1) RAT03530 | J=(U-BEG)/DELT-0.999999 RAT03540 | IF(J.LT.1) GO TO 26 RAT03550 | I=I-1 RAT03560 | DO 27 K=1,J RAT03570 I=I+1 RAT03580 | 27 ARG(I)=BEG+FLOAT(K)*DELT RAT03590 | 26 I=I+1 RAT03600 | IF(I.GT.1000) GO TO 6 RAT03610 | ARG(I)=U RAT03620 IF(M1.EQ.44) GO TO 228 RAT03630 | MG=M RAT03640 | M=M1 RAT03650 | GO TO 28 RAT03660 | 229 IF(M.NE.28) GO TO 26 RAT03670 6 K=2 RAT03680 | RETURN RAT03690 | 228 IF(M.EQ.60) GO TO 6 RAT03700 | K=1 RAT03710 | RETURN RAT03720 END RAT03730 | 1 12/01/64 PAGE 9 | | 0| - $IBFTC RECRD2 RAT03740 | 1 12/01/64 PAGE 10 | RECRD2 - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE RECORD (X,N1,N2,M1) RAT03750 COMMON/RECO/ITEST (80),N4 RAT03760 | DIMENSION L(5) RAT03770 | N2=48 RAT03780 | DO 71 I=1,5 RAT03790 | 71 L(I)=I RAT03800 I5=-1 RAT03810 | I10=1 RAT03820 | 150 I8=1 RAT03830 | 1 X=0. RAT03840 | Y=1. RAT03850 I9=1 RAT03860 | 2 GO TO 110 RAT03870 | 100 IF(M1.GE.10) GO TO 21 RAT03880 | X=10.*X+Y*FLOAT (M1) RAT03890 | GO TO (73,74),I9 RAT03900 73 I5=0 RAT03910 | 75 GO TO (76,110),I8 RAT03920 | 74 I5=I5+1 RAT03930 | 76 I8=2 RAT03940 | 28 L(1)=7 RAT03950 L(2)=7 RAT03960 | 29 L(4)=6 RAT03970 | GO TO 110 RAT03980 | 31 L(3)=7 RAT03990 | I9=2 RAT04000 GO TO 30 RAT04010 | 96 Y=-Y RAT04020 | 30 I5=0 RAT04030 | GO TO 28 RAT04040 | 41 IF(I8.EQ.1) GO TO 13 RAT04050 X1=X RAT04060 | I51=I5 RAT04070 | 17 L(1)=1 RAT04080 | L(2)=2 RAT04090 | L(3)=7 RAT04100 L(4)=4 RAT04110 | L(5)=7 RAT04120 | GO TO 150 RAT04130 | 21 M2=-1 RAT04140 | IF(M1.EQ.16)M2=1 RAT04150 IF(M1.EQ.32)M2=2 RAT04160 | IF(M1.EQ.27)M2=3 RAT04170 | IF(M1.EQ.48)M2=4 RAT04180 | IF(M1.EQ.21.AND.I5.NE.(-1)) M2=5 RAT04190 | I=L(M2) RAT04200 IF(M2.GE.0) GO TO (30,96,31,2,41,40,13),I RAT04210 | IF((M1.GE.17.AND.M1.LE.25).OR.(M1.GE.33.AND.M1.LE.41). RAT04220 | 1OR.(M1.GE.50.AND.M1.LE.57)) GO TO 154 RAT04230 | IF(I5.NE.(-1)) GO TO 40 RAT04240 | 24 IF(N2.NE.48) GO TO 13 RAT04250 N2=M1 RAT04260 | GO TO 110 RAT04270 | 154 IF(I5.NE.(-1)) GO TO 13 RAT04280 | N1=M1 RAT04290 | 1 12/01/64 PAGE 11 RECRD2 - EFN SOURCE STATEMENT - IFN(S) - | 0 I10=2 RAT04300 | DO 79I=1,5 RAT04310 | IF(M1.NE.48) GO TO 110 RAT04320 | 101 IF(N1.LT.536870912) GO TO 79 RAT04330 N1=64*(536870912-N1)-M1 RAT04340 | GO TO 80 RAT04350 | 79 N1=64*N1+M1 RAT04360 | 80 I10=3 RAT04370 | 102 IF(M1.NE.48) GO TO 110 RAT04380 160 M1=64 RAT04390 | RETURN RAT04400 | 40 IF(I8.NE.2) GO TO 13 RAT04410 | IF(L(5).EQ.5) GO TO 45 RAT04420 | I5=I51-INT(X) RAT04430 X=X1 RAT04440 | 45 X=X/10.**I5 RAT04450 97 | 6 RETURN RAT04460 | 13 M1=-M1 RAT04470 | RETURN RAT04480 110 IF(N4.GE.1.AND.N4.LT.80) GO TO 111 RAT04490 | READ (5,112)(ITEST(I),I=1,80) RAT04500 103 | 112 FORMAT (80A1) RAT04510 | N4=0 RAT04520 | 111 N4=N4+1 RAT04530 M1=ITEST (N4)/1073741824 RAT04540 | IF(ITEST(N4).LT.0)M1=32-M1 RAT04550 | GO TO (100,101,102),I10 RAT04560 | END RAT04570 | 0ERROR MESSAGE NUMBER 1 1 12/01/64 PAGE 12 | | 0| -| $IBFTC STABL RAT04580 1 12/01/64 PAGE 13 | STABL - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE STABLE (K1,K2,K3,IOP,IAD1,IAD2,IAD3) RAT04590 | DIMENSION ITE(19),IOP(1),IAD1(1),IAD2(1),IAD3(1),INTAB(320) RAT04600 DATA ITE/16,32,44,49,2860,60,28,6H000ABS,6H00AINT,6H000EXP,6H00ALORAT04610 | 2G,6HALOG10,6H00ATAN,6H000SIN,6H000COS,6H00TANH,6H00SQRT,6H0ARSIN, RAT04620 | 36H0ARCOS/ RAT04630 | K3=100 RAT04640 | DO 72 I1=1,4 RAT04650 J=79*I1-78 RAT04660 | K=J+79 RAT04670 | READ (5,1)(INTAB(I),I=J,K) RAT04680 7 | 1 FORMAT (80A1) RAT04690 | WRITE (6,2)(INTAB(I),I=J,K) RAT04700 12 2 FORMAT (1H0,80A1) RAT04710 | IF(INTAB(K).EQ.(-17997958192))GO TO 71 RAT04720 | 72 CONTINUE RAT04730 | GO TO 100 RAT04740 | 71 JSIST=0 RAT04750 DO 3 I=1,K RAT04760 | J=INTAB(I)/1073741824 RAT04770 | IF(INTAB(I).LT.0) J=32-J RAT04780 | IF(J.EQ.48) GO TO 3 RAT04790 | IF(JSIST.EQ.0.OR.J.NE.ITE(3)) GO TO 25 RAT04800 IF(INTAB(JSIST).NE.J) GO TO 25 RAT04810 | INTAB (JSIST)=ITE(5) RAT04820 | GO TO 3 RAT04830 | 25 JSIST=JSIST+1 RAT04840 | INTAB (JSIST)=J RAT04850 3 CONTINUE RAT04860 | IF(64*INTAB(1)+INTAB(2).NE.K1.OR.JSIST.LE.2) GO TO 100 RAT04870 | IPAR=64 RAT04880 | IIN=2 RAT04890 | IPLATS=1 RAT04900 INTAB(1)=-64 RAT04910 | 16 ISUM=0 RAT04920 | K=0 RAT04930 | 6 IIN=IIN+1 RAT04940 | IF(IIN.GT.JSIST) GO TO 50 RAT04950 DO 4 IS=1,7 RAT04960 | 4 IF(INTAB(IIN).EQ.ITE(IS)) GO TO 5 RAT04970 | ISUM=64*ISUM+INTAB(IIN) RAT04980 | K=K+1 RAT04990 | K3=101 RAT05000 IF(K.EQ.7) GO TO 100 RAT05010 | GO TO 6 RAT05020 | 5 IF(K.EQ.0) GO TO 7 RAT05030 | IUT=1 RAT05040 | 17 IF(ISUM.EQ.56) GO TO 18 RAT05050 IF(ISUM.NE.55) GO TO 19 RAT05060 | 18 I=54-ISUM RAT05070 | GO TO 9 RAT05080 | 19 DO 8 I=8,19 RAT05090 | 8 IF(ISUM.EQ.ITE(I)) GO TO 9 RAT05100 K3=102 RAT05110 | IF(K-3)26,10,100 RAT05120 | 26 I=ISUM/64 RAT05130 | 1 12/01/64 PAGE 14 | STABL - EFN SOURCE STATEMENT - IFN(S) - 0 ISUM=ISUM-64*I RAT05140 | GO TO 11 RAT05150 | 10 I=ISUM/4096 RAT05160 | ISUM=ISUM-4096*I RAT05170 | ISUM=ISUM-ISUM/64*54 RAT05180 11 IF(I.NE.19.OR.ISUM.GT.25) GO TO 100 RAT05190 | I=-2-ISUM RAT05200 | 9 IPLATS=IPLATS+1 RAT05210 | INTAB(IPLATS)=-I RAT05220 | GO TO (7,52),IUT RAT05230 7 K3=103 RAT05240 | I=INTAB(IPLATS) RAT05250 | GO TO (90,90,91,91,91,92,91),IS RAT05260 | 90 IF(I.GT.(-64).AND.I.LT.0) GO TO 100 RAT05270 | GO TO 13 RAT05280 91 IF(INTAB(IIN-1).NE.28.AND.I.LT.0) GO TO 100 RAT05290 | GO TO (100,100,13,13,13,100,93),IS RAT05300 | 92 IF(INTAB(IIN-1).EQ.28.OR.I.GT.0) GO TO 100 RAT05310 | IPAR=IPAR+1 RAT05320 | IS=IPAR RAT05330 GO TO 15 RAT05340 | 93 IS=IPAR RAT05350 | IPAR=IPAR-1 RAT05360 | 15 K3=104 RAT05370 | IF(IPAR.LT.64) GO TO 100 RAT05380 13 IPLATS=IPLATS+1 RAT05390 | INTAB(IPLATS)=-IS RAT05400 | GO TO 16 RAT05410 | 50 IUT=2 RAT05420 | IF(K.NE.0) GO TO 17 RAT05430 K3=105 RAT05440 | IF(INTAB(IIN-1).NE.28) GO TO 100 RAT05450 | GO TO 94 RAT05460 | 52 K3=105 RAT05470 | IF(I.GT.0) GO TO 100 RAT05480 94 IF(IPAR.NE.64) GO TO 100 RAT05490 | IPLATS=IPLATS+1 RAT05500 | INTAB(IPLATS)=-64 RAT05510 | K2=0 RAT05520 | 70 IL=1 RAT05530 DO 53 I=2,IPLATS RAT05540 | IF(INTAB(I).LT.INTAB(IL))IL=I RAT05550 | 53 CONTINUE RAT05560 | IL=IL+1 RAT05570 | DO 55 IS=IL,IPLATS RAT05580 IF(INTAB(IS).EQ.INTAB(IL-1)) GO TO 56 RAT05590 | 55 CONTINUE RAT05600 | 56 IS=IS-1 RAT05610 | ILL=IL RAT05620 | IF(INTAB(IL).LT.0) ILL=IL+1 RAT05630 K3=106 RAT05640 | IF(ILL.GT.IS.OR.INTAB(ILL).LT.0) GO TO 100 RAT05650 | DO 57 I=ILL,IS RAT05660 | IF(I.EQ.ILL) GO TO 57 RAT05670 | IF(INTAB(I-1)*INTAB(I)) 57,80,100 RAT05680 80 IF(INTAB(I-1).LT.0) GO TO 100 RAT05690 | 1 12/01/64 PAGE 15 | STABL - EFN SOURCE STATEMENT - IFN(S) - | 0 57 CONTINUE RAT05700 | IF(IL.EQ.ILL.OR.INTAB(IL).EQ.(-1)) GO TO 85 RAT05710 K3=108 RAT05720 | IF(INTAB(IL).NE.(-2)) GO TO 100 RAT05730 | K2=K2+1 RAT05740 | IOP(K2)=-2 RAT05750 | IAD 1 (K2)=INTAB(ILL) RAT05760 IAD 2 (K2)=0 RAT05770 | IAD 3 (K2)=K2+27 RAT05780 | INTAB (IL)=-1 RAT05790 | INTAB (ILL)=IAD 3 (K2) RAT05800 | 85 DO 58 I=1,5,2 RAT05810 62 DO 59 J=ILL,IS RAT05820 | IF(INTAB(J).GT.(I-6)) GO TO 59 RAT05830 | K3=107 RAT05840 | IF(INTAB(J).LT.(-5)) GO TO 100 RAT05850 | K2=K2+1 RAT05860 IOP(K2)=INTAB(J) RAT05870 | JD=J-1 RAT05880 | DO 63 JM=ILL,JD RAT05890 | IF(INTAB(JM).NE.0) JA=JM RAT05900 | 63 CONTINUE RAT05910 IAD1(K2)=INTAB(JA) RAT05920 | JD=J+1 RAT05930 | DO 60 JJ=JD,IS RAT05940 | IF(INTAB(JJ).NE.0) GO TO 61 RAT05950 | 60 CONTINUE RAT05960 61 IAD2(K2)=INTAB(JJ) RAT05970 | IAD3(K2)=K2+27 RAT05980 | INTAB(JA)=IAD3(K2) RAT05990 | INTAB(J)=0 RAT06000 | INTAB(JJ)=0 RAT06010 GO TO 62 RAT06020 | 59 CONTINUE RAT06030 | 58 CONTINUE RAT06040 | IF(ILL.EQ.2.OR.(ILL.EQ.3.AND.INTAB(2).EQ.(-1))) GO TO 101 RAT06050 | IF(IL.EQ.2) GO TO 64 RAT06060 IF(INTAB(IL-2).GE.(-7).OR.INTAB(IL-2).LE.(-64)) GO TO 64 RAT06070 | I=INTAB(IL-2) RAT06080 | 68 K2=K2+1 RAT06090 | IAD2(K2)=I RAT06100 | IAD1(K2)=INTAB(ILL) RAT06110 IOP(K2)=-1 RAT06120 | IF(ILL.EQ.IL) GO TO 65 RAT06130 | 65 IAD3(K2)=K2+27 RAT06140 | INTAB(IL)=0 RAT06150 | IF(I.NE.0)IL=IL-1 RAT06160 INTAB(IL-1)=IAD3(K2) RAT06170 | INTAB(IL)=0 RAT06180 | 66 INTAB(ILL)=0 RAT06190 | INTAB(IS+1)=0 RAT06200 | IF(IL.EQ.2) GO TO 101 RAT06210 GO TO 70 RAT06220 | 64 IF(ILL.NE.IL) INTAB(IL)=0 RAT06230 | 76 INTAB(IL-1)=INTAB(ILL) RAT06240 | GO TO 66 RAT06250 | 1 12/01/64 PAGE 16 STABL - EFN SOURCE STATEMENT - IFN(S) - | 0 100 K3=K3-89 RAT06260 | K2=0 RAT06270 | RETURN RAT06280 | 101 K3=0 RAT06290 RETURN RAT06300 | END RAT06310 | 1 12/01/64 PAGE 17 | | 0 -| $IBFTC OPER RAT06320 | 1 12/01/64 PAGE 18 | OPER - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE OPERA (IOP,IAD1,IAD2,IAD3,CS) RAT06330 | DIMENSION CS(1) RAT06340 | JOP=-IOP RAT06350 | IF(IAD2)4,1,50 RAT06360 | 50 GO TO (5,6,7,8,9),JOP RAT06370 5 CS(IAD3)=CS(IAD1)+CS(IAD2) RAT06380 | RETURN RAT06390 | 6 CS(IAD3)=CS(IAD1)-CS(IAD2) RAT06400 | RETURN RAT06410 | 7 CS(IAD3)=CS(IAD1)*CS(IAD2) RAT06420 RETURN RAT06430 | 8 CS(IAD3)=CS(IAD1)/CS(IAD2) RAT06440 | RETURN RAT06450 | 9 CS(IAD3)=CS(IAD1)**CS(IAD2) RAT06460 24 | RETURN RAT06470 1 GO TO (2,3),JOP RAT06480 | 2 CS(IAD3)=CS(IAD1) RAT06490 | RETURN RAT06500 | 3 CS(IAD3)=-CS(IAD1) RAT06510 | RETURN RAT06520 4 J=-7-IAD2 RAT06530 | X=(3.-FLOAT(2*JOP))*CS(IAD1) RAT06540 | GO TO (10,11,12,13,14,15,16,17,18,19,20,21),J RAT06550 | 10 CS(IAD3)=ABS(X) RAT06560 | RETURN RAT06570 11 CS(IAD3)=AINT(X) RAT06580 | RETURN RAT06590 | 12 CS(IAD3)=EXP(X) RAT06600 42 | RETURN RAT06610 | 13 CS(IAD3)=ALOG(X) RAT06620 45 RETURN RAT06630 | 14 CS(IAD3)=ALOG10(X) RAT06640 48 | RETURN RAT06650 | 15 CS(IAD3)=ATAN(X) RAT06660 51 | RETURN RAT06670 16 CS(IAD3)=SIN(X) RAT06680 54 | RETURN RAT06690 | 17 CS(IAD3)=COS(X) RAT06700 57 | RETURN RAT06710 | 18 CS(IAD3)=TANH(X) RAT06720 60 RETURN RAT06730 | 19 CS(IAD3)=SQRT(X) RAT06740 63 | RETURN RAT06750 | 20 CS(IAD3)=ARSIN(X) RAT06760 66 | RETURN RAT06770 21 CS(IAD3)=ARCOS(X) RAT06780 69 | RETURN RAT06790 | END RAT06800 | 1 12/01/64 PAGE 19 | 0| -| $IBFTC CHEC RAT06810 | 1 12/01/64 PAGE 20 | CHEC - EFN SOURCE STATEMENT - IFN(S) - 0| SUBROUTINE CHECK RAT06820 | 0COMMON X(1000),Y(1000),W(1000),JX,ZERO(25),FZERO(25), RAT06830 | 1COEF(25),ACC,BMPAR,CPAR,M,N,INDER,INDX0,INDXK,INDPOT,INDMN, RAT06840 | 2INDCHA,INDACC,MN2,MPTEST,YMEAN(5),ITERA,MNTEST,AM(25),AMPAR RAT06850 REAL KN1,KP1,KP3 RAT06860 | D=0. RAT06870 | EA=0. RAT06880 | ET=0. RAT06890 | I5=0 RAT06900 DO 1 K=6,JX RAT06910 | KN1=(X(K-3)-X(K-2))*(X(K-3)-X(K))/((X(K-4)-X(K-2))*(X(K-4)-X(K))) RAT06920 | KP1=(X(K-3)-X(K-4))*(X(K-3)-X(K))/((X(K-2)-X(K-4))*(X(K-2)-X(K))) RAT06930 | KP3=1.-KN1-KP1 RAT06940 | 0CN2=(X(K-3)-X(K-4))*(X(K-3)-X(K-2))*(X(K-3)-X(K-1))/ RAT06950 1((X(K-5)-X(K-4))*(X(K-5)-X(K-2))*(X(K-5)-X(K-1))) RAT06960 | 0CN1=(X(K-3)-X(K-5))*(X(K-3)-X(K-2))*(X(K-3)-X(K-1))/ RAT06970 | 1((X(K-4)-X(K-5))*(X(K-4)-X(K-2))*(X(K-4)-X(K-1))) RAT06980 | 0CP1=(X(K-3)-X(K-5))*(X(K-3)-X(K-4))*(X(K-3)-X(K-1))/ RAT06990 | 1((X(K-2)-X(K-5))*(X(K-2)-X(K-4))*(X(K-2)-X(K-1))) RAT07000 CP2=1.-CN2-CN1-CP1 RAT07010 | E3=W(K-3)*(Y(K-3)-KN1*Y(K-4)-KP1*Y(K-2)-KP3*Y(K)) RAT07020 | E4=W(K-3)*(Y(K-3)-CN2*Y(K-5)-CN1*Y(K-4)-CP1*Y(K-2)-CP2*Y(K-1)) RAT07030 | SIGMA=SQRT(CN2*CN2+CN1*CN1+CP1*CP1+CP2*CP2+1.) RAT07040 69 | AMY=CN1*KN1+1.+CP1*KP1 RAT07050 ETK=E4/SIGMA RAT07060 | EAK=E3-AMY*ETK/SIGMA RAT07070 | DK=2.*ABS(ETK)+0.25*ABS(EAK) RAT07080 | IF(D.GE.DK) GO TO 2 RAT07090 | D=DK RAT07100 2 IF(ABS(EA).GE.ABS(EAK)) GO TO 3 RAT07110 | EA=EAK RAT07120 | 3 IF(ABS(ET).GE.ABS(ETK)) GO TO 4 RAT07130 | ET=ETK RAT07140 | 4 IF(ACC.GT.0..AND.DK.GT.ACC) GO TO 9 RAT07150 1 CONTINUE RAT07160 | WRITE (6,6)D,EA,ET RAT07170 84 | 6 FORMAT (15H0MAXDIVERGENCES,20X,1P3E11.2) RAT07180 | IF((ACC.GT.0..AND.INDACC.NE.0).OR.ACC.GT.D)GO TO 7 RAT07190 | ACC=D RAT07200 7 WRITE(6,10)ACC RAT07210 89 | 10 FORMAT(9H0ACCURACY,1PE11.2) RAT07220 | RETURN RAT07230 | 9 IF(I5.EQ.0) WRITE(6,8) RAT07240 91 | 8 FORMAT (1H0,8X,1HX,16X,1HY,12X,3HDIV,9X,2HEA,9X,2HET) RAT07250 I5=1 RAT07260 | WRITE (6,5)X(K-3),Y(K-3),DK,EAK,ETK RAT07270 93 | 5 FORMAT (1H ,1P2E17.8,3E11.2) RAT07280 | GO TO 1 RAT07290 | END RAT07300 0ERROR MESSAGE NUMBER 1 | 1 12/01/64 PAGE 21 | | 0| - $IBFTC COMPU RAT07310 | 1 12/01/64 PAGE 22 | COMPU - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE COMPUT(L5,L) RAT07320 0COMMON X(1000),Y(1000),W(1000),JX,ZERO(25),FZERO(25), RAT07330 | 1COEF(25),ACC,BMPAR,CPAR,M,N,INDER,INDX0,INDXK,INDPOT,INDMN, RAT07340 | 2INDCHA,INDACC,MN2,MPTEST,YMEAN(5),ITERA,MNTEST,AM(25),AMPAR RAT07350 | YMEAN(4)=ACC RAT07360 | MNTEST=1 RAT07370 MN2=M+N+2 RAT07380 | GO TO (24,296),L5 RAT07390 | 24 MN2=M+N+2 RAT07400 | DO 297 I=2,MN2 RAT07410 | 297 ZERO(I)=0.5*((X(JX)-X(1))*COS(3.14159265*(1.-FLOAT(I-1)/ RAT07420 1FLOAT(MN2)))+X(JX)+X(1)) RAT07430 11 | 296 ZERO(1)=X(1) RAT07440 | ZERO(MN2+1)=X(JX) RAT07450 | IF(MN2.LT.4)GO TO 90 RAT07460 | IF(INDX0.NE.0)ZERO(2)=X(1) RAT07470 IF(INDXK.NE.0)ZERO(MN2)=X(JX) RAT07480 | 90 BMPAR=AMPAR RAT07490 | MPTEST=0 RAT07500 | ITERA=0 RAT07510 | IF(MN2.GE.24)GO TO 11 RAT07520 2 IF(INDER.EQ.1)WRITE(6,20)(ZERO(I),I=2,MN2) RAT07530 35 | 20 FORMAT (7H0ZEROES/( 1H ,1P8E15.7)) RAT07540 | CALL INTERP (L,MN2) RAT07550 43 | GO TO (1,100),L RAT07560 | 1 CALL EQUAT(I5) RAT07570 46 IF(I5.LT.0) GO TO 101 RAT07580 | CALL ERROR (L) RAT07590 51 | IF (INDER.EQ.0)GOTO23 RAT07600 | WRITE(6,21)(AM(I),I=1,MN2) RAT07610 55 | 21 FORMAT(7H ERRORS/( 1H ,1P8E15.7)) RAT07620 WRITE(6,22)YMEAN(3) RAT07630 62 | 22 FORMAT(11H MEANERROR ,1PE15.7) RAT07640 | 23 GO TO (2,3,15,16,4),L RAT07650 | 3 WRITE(6,25) RAT07660 64 | 25 FORMAT(32H0CONVERGENCE WITHIN THE ACCURACY) RAT07670 GO TO 26 RAT07680 | 15 WRITE(6,27) RAT07690 66 | 27 FORMAT(15H0NO CONVERGENCE) RAT07700 | GO TO 26 RAT07710 | 16 WRITE(6,28) RAT07720 68 28 FORMAT(40H0CONVERGENCE BUT NOT WITHIN THE ACCURACY) RAT07730 | 26 WRITE(6,29)X(1),X(JX),ITERA RAT07740 69 | 29 FORMAT(24H0APPROXIMATION INTERVAL ,1P2E15.7/ RAT07750 | 121H0NUMBER OF ITERATIONS,I5) RAT07760 | IF(INDPOT.EQ.0)GO TO 42 RAT07770 WRITE(6,41) RAT07780 73 | 41 FORMAT(22H0CHEBYSHEV POLYNOMIALS) RAT07790 | GO TO 40 RAT07800 | 42 WRITE(6,43) RAT07810 75 | 43 FORMAT(18H0POWER POLYNOMIALS) RAT07820 40 WRITE(6,20)(ZERO(I),I=2,MN2) RAT07830 76 | WRITE(6,21)(AM(I),I=1,MN2) RAT07840 83 | WRITE(6,22)YMEAN(3) RAT07850 90 | M2=M+1 RAT07855 | 1 12/01/64 PAGE 23 COMPU - EFN SOURCE STATEMENT - IFN(S) - | 0 WRITE(6,30)M,N,(COEF(I),I=1,M) RAT07860 91 | 30 FORMAT(8H0DEGREES,2I5/13H0COEFFICIENTS/( 1H 1P8E16.8)) RAT07870 | J=M+2 RAT07880 | WRITE(6,31)(COEF(I),I=J,MN2) RAT07890 98 31 FORMAT(1H 1P8E16.8) RAT07900 | IF(INDCHA.EQ.1.OR.L.LE.3)GO TO 10 RAT07910 | 14 GO TO (5,6,7),MNTEST RAT07920 | 5 M=M+1 RAT07930 | 8 MNTEST=MNTEST+1 RAT07940 GO TO (10,24,90), MNTEST RAT07950 | 6 YMEAN(5)=YMEAN(3) RAT07960 | M=M-1 RAT07970 | N=N+1 RAT07980 | GO TO 8 RAT07990 7 IF(YMEAN(5).GT.YMEAN(3))GO TO 19 RAT08000 | 140 M=M+1 RAT08010 | N=N-1 RAT08020 | 19 MNTEST=1 RAT08030 | GO TO 5 RAT08040 11 WRITE(6,12) RAT08050 122 | 12 FORMAT(17H0DEGREES TOO HIGH) RAT08060 | 100 L=3 RAT08070 | RETURN RAT08080 | 4 WRITE(6,13)M,N RAT08090 125 13 FORMAT(8H0DEGREES,2I4,19H GIVE WORSE RESULT) RAT08100 | GO TO 140 RAT08110 | 101 WRITE(6,102) M,N RAT08120 127 | 102 FORMAT(28H0SINGULAR MATRIX FOR DEGREES,2I5) RAT08130 | 10 IF(INDMN.EQ.0.OR.N.EQ.0) RETURN RAT08140 M=M+1 RAT08150 | N=N-1 RAT08160 | GO TO 90 RAT08170 | END RAT08180 | 1 12/01/64 PAGE 24 | 0| -| $IBFTC INTER RAT08190 | 1 12/01/64 PAGE 25 INTER - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE INTERP(L,L1) RAT08200 | 0COMMON X(1000),Y(1000),W(1000),JX,ZERO(25),FZERO(25), RAT08210 | 1COEF(25),ACC,BMPAR,CPAR,M,N,INDER,INDX0,INDXK,INDPOT,INDMN, RAT08220 2INDCHA,INDACC,MN2,MPTEST,YMEAN(5),ITERA,MNTEST,AM(25),AMPAR RAT08230 | DIMENSION IPLATS(25) RAT08240 | L=1 RAT08250 | Z1=X(1)-X(2) RAT08260 | K=3 RAT08270 DO1I=2,L1 RAT08280 | DO2J=K,JX RAT08290 | IF(Z1*(ZERO(I)-X(J-1)).GT.0.)GO TO 3 RAT08300 | 2 CONTINUE RAT08310 | J=JX RAT08320 3 IPLATS(I)=J RAT08330 | K=J RAT08340 | IF(I.EQ.2)GO TO 1 RAT08350 | IF(IPLATS(I).NE.IPLATS(I-1))GO TO 1 RAT08360 | WRITE(6,5)ZERO(I-1),ZERO(I) RAT08370 26 5 FORMAT(7H0ZEROES,1P2E15.4,4X,21H BETWEEN GIVEN POINTS) RAT08380 | L=2 RAT08390 | 1 CONTINUE RAT08400 | GO TO (7,8),L RAT08410 | 7 DO10I=2,L1 RAT08420 J=IPLATS(I) RAT08430 | 0Y01=(Y(J-2)*(X(J-1)-ZERO(I))-Y(J-1)*(X(J-2)-ZERO(I)))/ RAT08440 | 1(X(J-1)-X(J-2)) RAT08450 | 0Y02=(Y(J-2)*(X(J)-ZERO(I))-Y(J)*(X(J-2)-ZERO(I)))/ RAT08460 | 1(X(J)-X(J-2)) RAT08470 100FZERO(I)=(Y01*(X(J)-ZERO(I))-Y02*(X(J-1)-ZERO(I)))/ RAT08480 | 1(X(J)-X(J-1)) RAT08490 | 8 RETURN RAT08500 | END RAT08510 | 1 12/01/64 PAGE 26 | 0| -| $IBFTC ERRROR RAT08520 | 1 12/01/64 PAGE 27 ERRROR - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE ERROR(L) RAT08530 | 0COMMON X(1000),Y(1000),W(1000),JX,ZERO(25),FZERO(25), RAT08540 | 1COEF(25),ACC,BMPAR,CPAR,M,N,INDER,INDX0,INDXK,INDPOT,INDMN, RAT08550 2INDCHA,INDACC,MN2,MPTEST,YMEAN(5),ITERA,MNTEST,AM(25),AMPAR RAT08560 | K=1 RAT08570 | L1=1 RAT08580 | L2=MN2 RAT08590 | AM(1)=0. RAT08600 AM(MN2)=0. RAT08610 | IF(MN2.LT.4)GO TO 36 RAT08620 | IF(INDX0.NE.0)L1=2 RAT08630 | IF(INDXK.NE.0)L2=MN2-1 RAT08640 | 36 DO 1 I=L1,L2 RAT08650 Y2=0. RAT08660 | DO 2 J=K,JX RAT08670 | IF(ZERO(I+1).LT.X(J))GO TO 3 RAT08680 | F=FAPP(X(J))-Y(J) RAT08690 25 | IF(ABS(F).LT.ABS(Y2)) GO TO 2 RAT08700 JSPAR=J RAT08710 | Y2=F RAT08720 | 2 CONTINUE RAT08730 | 3 K=J RAT08740 | IF(JSPAR.EQ.1.OR.JSPAR.EQ.JX) GO TO 16 RAT08750 X1=X(JSPAR-1)-X(JSPAR) RAT08760 | Y1=FAPP(X(JSPAR-1))-Y(JSPAR-1)-Y2 RAT08770 40 | X3=X(JSPAR+1)-X(JSPAR) RAT08780 | Y3=FAPP(X(JSPAR+1))-Y(JSPAR+1)-Y2 RAT08790 45 | 0AM(I)=ABS(Y2+(X1*X1*Y3-X3*X3*Y1)**2/ RAT08800 1(4.*(X1*Y3-X3*Y1)*X1*X3*(X1-X3))) RAT08810 | 1 CONTINUE RAT08820 | ITERA=ITERA+1 RAT08830 | MPTEST=MPTEST+1 RAT08840 | YMEAN(1)=YMEAN(2) RAT08850 YMEAN(2)=YMEAN(3) RAT08860 | IS=0 RAT08870 | 23 YS=YMEAN(4) RAT08880 | 20 SI=0. RAT08890 | DO 18 I=L1,L2 RAT08900 18 SI=SI+(ZERO(I+1)-ZERO(I))/((BMPAR+1.)*AM(I)+YS) RAT08910 | YMEAN(3)=((ZERO(MN2+1)-ZERO(1))/SI-YS)/(BMPAR+1.) RAT08920 | IF(ABS(YMEAN(3)/YS-1.).LT.0.001) GO TO 19 RAT08930 | YS=YMEAN(3) RAT08940 | GO TO 20 RAT08950 190IF(MPTEST.LT.3.OR.ABS(YMEAN(3)-YMEAN(2)).LT. RAT08960 | 1ABS(YMEAN(2)-YMEAN(1)).OR.IS.EQ.1) GO TO 21 RAT08970 | IS=1 RAT08980 | BMPAR=0.5*BMPAR RAT08990 | MPTEST=1 RAT09000 WRITE(6,22)BMPAR,ITERA RAT09010 74 | 220FORMAT(27H0MOVING PARAMETER HALVED TO,1PE12.3,3X, RAT09020 | 15HAFTER,I5,3X,10HITERATIONS) RAT09030 | GO TO 23 RAT09040 | 21 IF(MNTEST.NE.3) GO TO 30 RAT09050 DO 29 I=L1,L2 RAT09060 | IF(AM(I).LT.YMEAN(5)) GO TO 30 RAT09070 | 1 12/01/64 PAGE 28 | ERRROR - EFN SOURCE STATEMENT - IFN(S) - | 0 29 CONTINUE RAT09080 C BYT GRADTAL RAT09090 | L=5 RAT09100 | GO TO 28 RAT09110 | 30 DO 24 I=L1,L2 RAT09120 | IF(AM(I).GT.ACC) GO TO 25 RAT09130 24 CONTINUE RAT09140 | C SLUTKONVERGENS RAT09150 | L=2 RAT09160 | GO TO 28 RAT09170 | 25 DO 26 I=L1,L2 RAT09180 IF((CPAR*YMEAN(3)+0.125*ACC).LT.AM(I))GO TO 27 RAT09190 | 26 CONTINUE RAT09200 | C KONVERGENS RAT09210 | L=4 RAT09220 | GO TO 28 RAT09230 27 IF(ITERA.LE.100)GO TO 40 RAT09240 | C OEVER 100 ITERATIONER RAT09250 | L=3 RAT09260 | GO TO 28 RAT09270 | C EJ KONVERGENS RAT09280 40 L=1 RAT09290 | L1=L1+1 RAT09300 | SPAR=ZERO(1) RAT09310 | DO 35 I=L1,L2 RAT09320 | 0CL=(AM(I-1)+YS+BMPAR*YMEAN(3))/ RAT09330 1((BMPAR+1.)*AM(I-1)+YS)*(ZERO(I)-SPAR) RAT09340 | SPAR=ZERO(I) RAT09350 | 35 ZERO(I)=CL+ZERO(I-1) RAT09360 | YMEAN(4)=YMEAN(3) RAT09370 | 28 RETURN RAT09380 16 AM(I)=ABS(Y2) RAT09390 | GO TO 1 RAT09400 | END RAT09410 | 0ERROR MESSAGE NUMBER 1 | 1 12/01/64 PAGE 29 | 0| -| $IBFTC FAP RAT09420 | 1 12/01/64 PAGE 30 FAP - EFN SOURCE STATEMENT - IFN(S) - | 0| FUNCTION FAPP(XARG) RAT09430 | 0COMMON X(1000),Y(1000),W(1000),JX,ZERO(25),FZERO(25), RAT09440 | 1COEF(25),ACC,BMPAR,CPAR,M,N,INDER,INDX0,INDXK,INDPOT,INDMN, RAT09450 2INDCHA,INDACC,MN2,MPTEST,YMEAN(5),ITERA,MNTEST,AM(25),AMPAR RAT09460 | DIMENSION B(26) RAT09470 | I=INDPOT+1 RAT09480 | GO TO (4,5),I RAT09490 | 4 AT=0. RAT09500 M1=M-1 | DO 6 II=1,M1 RAT09510 | IJ=M-(II-1) RAT09520 | 6 AT=XARG*AT+COEF(IJ+1) RAT09530 | AN=0. RAT09540 N1=N-1 | DO 7 II=1,N1 RAT09550 | IJ=MN2-(II-1) RAT09560 | 7 AN=XARG*AN+COEF(IJ) RAT09570 | GO TO 8 RAT09580 5 Z=(2.*XARG-X(JX)-X(1))/(X(JX)-X(1)) RAT09590 | IF(M.EQ.0) GO TO 9 RAT09600 | B(M+2)=0. RAT09610 | B(M+1)=0. RAT09620 | DO 10 II=1,M RAT09630 IJ=M-II RAT09640 | 10 B(IJ+1)=2.*Z*B(IJ+2)-B(IJ+3)+COEF(IJ+2) RAT09650 | AT=Z*B(1)-B(2)+COEF(1) RAT09660 | 19 IF(N.EQ.0) GO TO 11 RAT09670 | B(MN2+2)=0. RAT09680 B(MN2+1)=0. RAT09690 | DO 12 II=1,N RAT09700 | IJ=MN2-II RAT09710 | 12 B(IJ+1)=2.*Z*B(IJ+2)-B(IJ+3)+COEF(IJ+1) RAT09720 | AN=Z*B(M+3)-B(M+4)+COEF(M+2) RAT09730 8 FAPP=AT/AN RAT09740 | RETURN RAT09750 | 9 AT=COEF(1) RAT09760 | GO TO 19 RAT09770 | 11 FAPP=AT RAT09780 RETURN RAT09790 | END RAT09800 | 1 12/01/64 PAGE 31 | | 0 -| $IBFTC EQU RAT09810 | 1 12/01/64 PAGE 32 | EQU - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE EQUAT(IERR) RAT09820 | 0COMMON X(1000),Y(1000),W(1000),JX,ZERO(25),FZERO(25), RAT09830 | 1COEF(25),ACC,BMPAR,CPAR,M,N,INDER,INDX0,INDXK,INDPOT,INDMN, RAT09840 | 2INDCHA,INDACC,MN2,MPTEST,YMEAN(5),ITERA,MNTEST,AM(25),AMPAR RAT09850 | DIMENSION CMAT(25,25),CF(25,1) RAT09860 EQUIVALENCE (COEF(1),CF(1,1)) RAT09870 | IPOT=INDPOT+1 RAT09880 | MNMAX=MAX0(M,N) RAT09890 | IF(MNMAX.EQ.0) GO TO 11 RAT09900 | DO 1 I=2,MN2 RAT09910 CMAT(I-1,1)=1. RAT09920 | F1=1. RAT09930 | GO TO (14,15),IPOT RAT09940 | 14 F2=ZERO(I) RAT09950 | GO TO 16 RAT09960 15 Z=(2.*ZERO(I)-X(JX)-X(1))/(X(JX)-X(1)) RAT09970 | F2=Z RAT09980 | 16 DO 2 J=1,MNMAX RAT09990 | IF(J.EQ.1) GO TO 6 RAT10000 | GO TO (4,5), IPOT RAT10010 4 F2=F2*ZERO(I) RAT10020 | GO TO 6 RAT10030 | 5 F3=2.*Z*F2-F1 RAT10040 | F1=F2 RAT10050 | F2=F3 RAT10060 6 IF(J.GT.M) GO TO 7 RAT10070 | CMAT(I-1,J+1)=F2 RAT10080 | 7 IF(J.GT.N) GO TO 2 RAT10090 | JM=J+M RAT10100 | CMAT(I-1,JM+1)=-F2*FZERO(I) RAT10110 2 CONTINUE RAT10120 | 1 CMAT(I-1,MN2)=FZERO(I) RAT10130 | CALL MATS (MN2-1,1,CMAT,25,25,CF,25,1,IERR) RAT10140 47 | IF(N.EQ.0.OR.IERR.EQ.(-1)) GO TO 8 RAT10150 | DO 9 J=1,N RAT10160 JJ=MN2-J RAT10170 | 9 COEF(JJ+1)=COEF(JJ) RAT10180 | 8 COEF(M+2)=1. RAT10190 | RETURN RAT10200 | 11 IERR=1 RAT10210 10 COEF(1)=FZERO(2) RAT10220 | GO TO 8 RAT10230 | END RAT10240 | 1 12/01/64 PAGE 33 | 0| -| $IBFTC MATRS RAT10250 | 1 12/01/64 PAGE 34 | MATRS - EFN SOURCE STATEMENT - IFN(S) - 0| SUBROUTINE MATS(N1,M1,A,ARAD,AKOL,X,XRAD,XKOL,IERR) RAT10260 | INTEGER N1,M1,ARAD,AKOL,XRAD,XKOL,IERR RAT10270 | DIMENSION A(ARAD,AKOL),X(XRAD,XKOL) RAT10280 | N=MAX0(1,MIN0(ARAD,XRAD,N1)) RAT10290 M=MAX0(1,MIN0(AKOL-N,XKOL,M1)) RAT10300 | IF(N.EQ.N1.AND.M.EQ.M1) GO TO 1965 RAT10310 | WRITE(6,1966) N1,M1,ARAD,AKOL,XRAD,XKOL RAT10320 6 | 1966 FORMAT(27H ILLEGAL ARGUMENTS TO MATS. /3H N=I3,4H, M=I3,4H, A(I3, RAT10330 | 11H,I3,5H), X(I3,1H,I3,1H)) RAT10340 GO TO 1967 RAT10350 | 1965 MM=N+M RAT10360 | DO 15 I=2,N RAT10370 | 70 II=I-1 RAT10380 | 7 DO 15 J=1,II RAT10390 8 IF (A(I,J))9,15,9 RAT10400 | 9 IF (ABS (A(J,J))-ABS (A(I,J)))11,10,10 RAT10410 | 10 R=A(I,J)/A(J,J) RAT10420 | GO TO 130 RAT10430 | 11 R=A(J,J)/A(I,J) RAT10440 DO 12 K=1,MM RAT10450 | B=A(J,K) RAT10460 | A(J,K)=A(I,K) RAT10470 | 12 A(I,K)=B RAT10480 | 130 JJ=J+1 RAT10490 13 DO 14 K=JJ,MM RAT10500 | 14 A(I,K)=A(I,K)-R*A(J,K) RAT10510 | 15 CONTINUE RAT10520 | IF (ABS (A(N,N))-1.0E-10)16,16,17 RAT10530 | 16 WRITE(6,1) RAT10540 54 1 FORMAT(34H ERROR IN 'MATS'. MATRIX SINGULAR. ) RAT10550 | 1967 IERR=-1 RAT10560 | GO TO 1968 RAT10570 | 17 DO28J=1,M RAT10580 | KK=N+J RAT10590 X(N,J)=A(N,KK)/A(N,N) RAT10600 | DO28I=2,N RAT10610 | JJ=N-I+1 RAT10620 | B=0. RAT10630 | II=N-I+2 RAT10640 DO 25 K=II,N RAT10650 | 25 B=B+A(JJ,K)*X(K,J) RAT10660 | IF (ABS (A(JJ,JJ))-1.0E-10)16,16,28 RAT10670 | 28 X(JJ,J)=(A(JJ,KK)-B)/A(JJ,JJ) RAT10680 | IERR=1 RAT10690 1968 RETURN RAT10700 | END RAT10710 | 1 12/01/64 PAGE 35 | INPUT STORAGE MAP | 0 MAIN PROGRAM - COMMON VARIABLES | 0 COMMON BLOCK LAES ORIGIN 00001 LENGTH 00004 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | U 00000 R N1 00001 I N2 00002 I | M1 00003 I 0 COMMON BLOCK RECO ORIGIN 00005 LENGTH 00121 | ITEST 00000 I N4 00120 I | 0 COMMON BLOCK // ORIGIN 00126 LENGTH 06063 | X 00000 R Y 01750 R W 03720 R | JX 05670 I ZERO 05671 R FZERO 05722 R COEF 05753 R ACC 06004 R BMPAR 06005 R | CPAR 06006 R M 06007 I N 06010 I | INDER 06011 I INDX0 06012 I INDXK 06013 I | INDPOT 06014 I INDMN 06015 I INDCHA 06016 I | INDACC 06017 I MN2 06020 I MPTEST 06021 I YMEAN 06022 R ITERA 06027 I MNTEST 06030 I | AM 06031 R AMPAR 06062 R | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | LT 06211 I ARG 06230 R FUNC 10200 R IOPW 12150 I IAD1W 12326 I IAD2W 12504 I | IAD3W 12662 I IOPX 13040 I IAD1X 13216 I | IAD2X 13374 I IAD3X 13552 I IOPY 13730 I | IAD1Y 14106 I IAD2Y 14264 I IAD3Y 14442 I | CSX 14620 R CSY 15031 R CSW 15242 R XP 15453 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | I 17423 I IFEL 17424 I JY 17425 I | JJJ 17426 I IPRES 17427 I MDEG 17430 I NDEG 17431 I ACCIN 17432 R K2X 17433 I | K2Y 17434 I K2W 17435 I IPWEI 17436 I | IPSTAB 17437 I IPZER 17440 I K 17441 I | V 17442 R UU 17443 R J 17444 I | Z1 17445 R L 17446 I I1 17447 I SPAR1 17450 R SPAR2 17451 R K1 17452 I | K2 17453 I K3 17454 I | - ENTRY POINTS | 0 ...... SECTION 5 | 1 12/01/64 PAGE 36 INPUT STORAGE MAP | 0 SUBROUTINES CALLED | 0 RECORD SECTION 6 .FWRD. SECTION 7 .FRDD. SECTION 8 | .FSLI. SECTION 9 .FSLO. SECTION 10 REC SECTION 11 | OPERA SECTION 12 CHECK SECTION 13 COMPUT SECTION 14 INTERP SECTION 15 FAPP SECTION 16 STABLE SECTION 17 | .FXEM. SECTION 18 .UN06. SECTION 19 .FFIL. SECTION 20 | .FCNV. SECTION 21 .UN05. SECTION 22 .FRTN. SECTION 23 | CC.1 SECTION 24 CC.2 SECTION 25 CC.3 SECTION 26 | CC.4 SECTION 27 SYSLOC SECTION 28 - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 5 1A 17677 1 2A 17700 7 8A 17713 | 2 16A 17724 3 65A 20234 6 19A 17730 | 401 20A 17765 402 22A 17774 403 24A 20003 404 26A 20016 405 28A 20025 406 30A 20034 | 407 32A 20043 408 34A 20052 409 36A 20061 | 410 38A 20070 411 45A 20125 420 53A 20173 | 421 55A 20202 423 57A 20211 424 59A 20220 | 441 FORMAT 17515 442 FORMAT 17521 443 FORMAT 17526 8 61A 20227 444 FORMAT 17534 445 FORMAT 17540 | 446 FORMAT 17545 447 FORMAT 17551 448 FORMAT 17556 | 449 FORMAT 17563 490 39A 20100 491 41A 20107 | 492 43A 20116 493 FORMAT 17570 494 FORMAT 17575 | 495 FORMAT 17602 480 47A 20140 481 49A 20151 482 51A 20162 483 FORMAT 17607 484 FORMAT 17615 | 485 FORMAT 17624 460 FORMAT 17633 461 FORMAT 17640 | 463 FORMAT 17647 464 FORMAT 17653 10 69A 20262 | 20 88A 20335 30 149A 20541 40 180A 20662 | 50 182A 20666 80 198A 20730 90 208A 20762 100 345A 21375 110 418A 21663 120 476A 22103 | 130 480A 22120 11 FORMAT 17657 12 FORMAT 17660 | 25 121A 20441 22 99A 20353 29 103A 20373 | 24 110A 20416 220 145A 20530 23 122A 20442 | 239 156A 20561 39 162A 20577 31 169A 20627 32 172A 20642 33 174A 20645 34 176A 20654 | 35 178A 20657 51 194A 20725 81 FORMAT 17662 | 299 209A 20764 91 239A 21064 92 233A 21054 | 93 224A 21024 94 248A 21074 95 243A 21070 | 96 273A 21163 97 267A 21153 98 258A 21123 99 282A 21173 290 277A 21167 291 307A 21262 | 292 301A 21252 293 292A 21222 294 317A 21276 | 295 313A 21272 298 327A 21323 297 331A 21334 | 309 368A 21437 101 355A 21415 102 359A 21422 | 103 371A 21454 300 373A 21470 301 375A 21477 302 377A 21506 303 FORMAT 17663 306 378A 21514 | 304 FORMAT 17666 305 FORMAT 17670 109 379A 21516 | 107 401A 21617 104 387A 21554 105 392A 21567 | 106 398A 21607 108 FORMAT 17673 307 416A 21656 | 119 455A 22025 112 424A 21674 118 427A 21710 113 434A 21724 114 436A 21737 115 440A 21753 | 116 444A 21767 117 447A 22002 311 474A 22077 | 310 460A 22034 | 1 12/01/64 PAGE 37 | INPUT STORAGE MAP 0 THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 22162.| 1 12/01/64 PAGE 38 | RECRD1 STORAGE MAP | 0 SUBROUTINE REC | - COMMON VARIABLES 0 COMMON BLOCK LAES ORIGIN 00001 LENGTH 00004 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | U 00000 R N1 00001 I N2 00002 I | M1 00003 I | - DIMENSIONED PROGRAM VARIABLES 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | ISIG 00005 I | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | MG 00011 I M 00012 I L 00013 I DELT 00014 R BEG 00015 R J 00016 I | - ENTRY POINTS | 0 REC SECTION 3 | - SUBROUTINES CALLED | 0 RECORD SECTION 4 E.1 SECTION 5 E.2 SECTION 6 E.3 SECTION 7 E.4 SECTION 8 CC.1 SECTION 9 | CC.2 SECTION 10 CC.3 SECTION 11 CC.4 SECTION 12 | SYSLOC SECTION 13 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 22 2A 00034 28 4A 00040 1 14A 00060 | 2 18A 00063 6 62A 00236 25 24A 00100 | 229 59A 00232 26 48A 00206 27 44A 00166 | 228 64A 00241 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00322. 1 12/01/64 PAGE 39 | RECRD2 STORAGE MAP | 0 SUBROUTINE RECORD | - COMMON VARIABLES | 0 COMMON BLOCK RECO ORIGIN 00001 LENGTH 00121 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | ITEST 00000 I N4 00120 I | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | L 00122 I - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | I 00127 I I5 00130 I I10 00131 I | I8 00132 I Y 00133 R I9 00134 I | X1 00135 R I51 00136 I M2 00137 I - ENTRY POINTS | 0 RECORD SECTION 3 | - SUBROUTINES CALLED | 0 .XP2. SECTION 4 .FRDD. SECTION 5 .FXEM. SECTION 6 | .UN05. SECTION 7 .FRTN. SECTION 8 .FCNV. SECTION 9 E.1 SECTION 10 E.2 SECTION 11 E.3 SECTION 12 | E.4 SECTION 13 CC.1 SECTION 14 CC.2 SECTION 15 | CC.3 SECTION 16 CC.4 SECTION 17 SYSLOC SECTION 18 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 71 6A 00211 150 11A 00221 1 12A 00223 | 2 13A 00230 110 101A 00605 100 14A 00231 | 21 35A 00343 73 18A 00263 74 20A 00273 | 75 19A 00264 76 21A 00276 28 22A 00300 | 29 23A 00304 31 25A 00307 30 28A 00316 96 27A 00314 41 30A 00320 13 99A 00602 | 17 33A 00330 40 89A 00546 154 64A 00470 | 24 60A 00461 79 79A 00526 101 75A 00511 | 80 83A 00535 102 84A 00537 160 87A 00543 | 45 96A 00570 6 98A 00601 111 110A 00632 112 FORMAT 00201 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00775.| 1 12/01/64 PAGE 40 | STABL STORAGE MAP | 0 SUBROUTINE STABLE - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | ITE 00001 I INTAB 00024 I | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE I1 00524 I J 00525 I K 00526 I | I 00527 I JSIST 00530 I IPAR 00531 I | IIN 00532 I IPLATS 00533 I ISUM 00534 I | IS 00535 I IUT 00536 I IL 00537 I | ILL 00540 I JD 00541 I JM 00542 I JA 00543 I JJ 00544 I | - ENTRY POINTS | 0 STABLE SECTION 2 | - SUBROUTINES CALLED | 0 .FRDD. SECTION 3 .FWRD. SECTION 4 .FXEM. SECTION 5 .UN05. SECTION 6 .FRTN. SECTION 7 .FCNV. SECTION 8 | .UN06. SECTION 9 .FFIL. SECTION 10 SYSLOC SECTION 11 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 72 21A 00677 1 FORMAT 00612 2 FORMAT 00613 71 25A 00702 100 336A 02061 3 50A 00763 | 25 47A 00753 16 59A 01010 6 61A 01012 | 50 150A 01335 4 69A 01032 5 84A 01060 | 7 120A 01213 17 87A 01066 18 93A 01076 | 19 95A 01102 9 116A 01176 8 99A 01107 26 107A 01123 10 110A 01137 11 113A 01163 | 52 161A 01353 90 124A 01234 91 128A 01246 | 92 133A 01273 13 146A 01326 93 139A 01312 | 15 141A 01317 94 166A 01362 70 171A 01375 | 53 180A 01420 55 192A 01444 56 195A 01446 57 219A 01526 80 215A 01521 85 240A 01571 | 58 289A 01741 62 244A 01604 59 287A 01736 | 63 265A 01656 60 276A 01704 61 278A 01706 | 101 339A 02066 64 328A 02047 68 304A 01776 | 65 313A 02015 66 321A 02036 76 332A 02056 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 02215.| 1 12/01/64 PAGE 41 | OPER STORAGE MAP | 0 SUBROUTINE OPERA - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | JOP 00001 I J 00002 I X 00003 R | - ENTRY POINTS | 0 OPERA SECTION 2 - SUBROUTINES CALLED | 0 .XP3. SECTION 3 EXP SECTION 4 ALOG SECTION 5 | ALOG10 SECTION 6 ATAN SECTION 7 SIN SECTION 8 | COS SECTION 9 TANH SECTION 10 SQRT SECTION 11 | ARSIN SECTION 12 ARCOS SECTION 13 .FXEM. SECTION 14 E.1 SECTION 15 E.2 SECTION 16 E.3 SECTION 17 | E.4 SECTION 18 CC.1 SECTION 19 CC.2 SECTION 20 | CC.3 SECTION 21 CC.4 SECTION 22 SYSLOC SECTION 23 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 4 33A 00131 1 26A 00110 50 4A 00022 | 5 5A 00034 6 9A 00043 7 13A 00052 | 8 17A 00061 9 21A 00070 2 27A 00117 | 3 30A 00124 10 37A 00175 11 39A 00202 | 12 41A 00210 13 44A 00217 14 47A 00226 15 50A 00235 16 53A 00244 17 56A 00253 | 18 59A 00262 19 62A 00271 20 65A 00300 | 21 68A 00307 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00455.| 1 12/01/64 PAGE 42 CHEC STORAGE MAP | 0 SUBROUTINE CHECK | - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 06063 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE X 00000 R Y 01750 R W 03720 R | JX 05670 I ZERO 05671 R FZERO 05722 R | COEF 05753 R ACC 06004 R BMPAR 06005 R | CPAR 06006 R M 06007 I N 06010 I | INDER 06011 I INDX0 06012 I INDXK 06013 I INDPOT 06014 I INDMN 06015 I INDCHA 06016 I | INDACC 06017 I MN2 06020 I MPTEST 06021 I | YMEAN 06022 R ITERA 06027 I MNTEST 06030 I | AM 06031 R AMPAR 06062 R | - UNDIMENSIONED PROGRAM VARIABLES 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | KN1 06064 R KP1 06065 R KP3 06066 R | D 06067 R EA 06070 R ET 06071 R | I5 06072 I K 06073 I CN2 06074 R | CN1 06075 R CP1 06076 R CP2 06077 R E3 06100 R E4 06101 R SIGMA 06102 R | AMY 06103 R ETK 06104 R EAK 06105 R | DK 06106 R | - ENTRY POINTS | 0 CHECK SECTION 3 - SUBROUTINES CALLED | 0 SQRT SECTION 4 .FWRD. SECTION 5 .UN06. SECTION 6 | .FFIL. SECTION 7 .FCNV. SECTION 8 SYSLOC SECTION 9 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 1 82A 06565 2 73A 06526 3 76A 06541 | 4 79A 06554 9 90A 06632 6 FORMAT 06127 | 7 89A 06621 10 FORMAT 06135 8 FORMAT 06141 | 5 FORMAT 06151 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 06715. 1 12/01/64 PAGE 43 | COMPU STORAGE MAP | 0 SUBROUTINE COMPUT | - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 06063 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | X 00000 R Y 01750 R W 03720 R | JX 05670 I ZERO 05671 R FZERO 05722 R | COEF 05753 R ACC 06004 R BMPAR 06005 R | CPAR 06006 R M 06007 I N 06010 I INDER 06011 I INDX0 06012 I INDXK 06013 I | INDPOT 06014 I INDMN 06015 I INDCHA 06016 I | INDACC 06017 I MN2 06020 I MPTEST 06021 I | YMEAN 06022 R ITERA 06027 I MNTEST 06030 I | AM 06031 R AMPAR 06062 R - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | I 06064 I I5 06065 I M2 06066 I | J 06067 I | - ENTRY POINTS 0 COMPUT SECTION 3 | - SUBROUTINES CALLED | 0 COS SECTION 4 .FWRD. SECTION 5 INTERP SECTION 6 | EQUAT SECTION 7 ERROR SECTION 8 .FXEM. SECTION 9 | .UN06. SECTION 10 .FFIL. SECTION 11 .FCNV. SECTION 12 E.1 SECTION 13 E.2 SECTION 14 E.3 SECTION 15 | E.4 SECTION 16 CC.1 SECTION 17 CC.2 SECTION 18 | CC.3 SECTION 19 CC.4 SECTION 20 SYSLOC SECTION 21 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 24 4A 06257 296 16A 06352 297 9A 06277 | 90 28A 06401 11 122A 07011 2 34A 06412 | 20 FORMAT 06112 1 45A 06447 100 123A 07017 | 101 127A 07035 23 63A 06516 21 FORMAT 06117 | 22 FORMAT 06124 3 64A 06530 15 66A 06537 16 68A 06546 4 125A 07022 25 FORMAT 06130 | 26 69A 06554 27 FORMAT 06137 28 FORMAT 06143 | 29 FORMAT 06153 42 75A 06604 41 FORMAT 06172 | 40 76A 06612 43 FORMAT 06177 30 FORMAT 06203 | 31 FORMAT 06213 10 128A 07047 14 108A 06734 5 109A 06744 6 112A 06762 7 115A 06773 | 8 110A 06747 19 120A 07006 140 118A 07000 | 1 12/01/64 PAGE 44 | COMPU STORAGE MAP | 012 FORMAT 06216 13 FORMAT 06222 102 FORMAT 06231 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 07123.| 1 12/01/64 PAGE 45 | INTER STORAGE MAP | 0 SUBROUTINE INTERP - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 06063 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | X 00000 R Y 01750 R W 03720 R | JX 05670 I ZERO 05671 R FZERO 05722 R COEF 05753 R ACC 06004 R BMPAR 06005 R | CPAR 06006 R M 06007 I N 06010 I | INDER 06011 I INDX0 06012 I INDXK 06013 I | INDPOT 06014 I INDMN 06015 I INDCHA 06016 I | INDACC 06017 I MN2 06020 I MPTEST 06021 I YMEAN 06022 R ITERA 06027 I MNTEST 06030 I | AM 06031 R AMPAR 06062 R | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | IPLATS 06064 I - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | Z1 06115 R K 06116 I I 06117 I | J 06120 I Y01 06121 R Y02 06122 R | - ENTRY POINTS 0 INTERP SECTION 3 | - SUBROUTINES CALLED | 0 .FWRD. SECTION 4 .FXEM. SECTION 5 .UN06. SECTION 6 | .FFIL. SECTION 7 .FCNV. SECTION 8 SYSLOC SECTION 9 | - EFN IFN CORRESPONDENCE 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 1 29A 06244 2 14A 06210 3 16A 06214 | 5 FORMAT 06136 7 33A 06255 8 65A 06406 | 10 55A 06353 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 06443. 1 12/01/64 PAGE 46 | ERRROR STORAGE MAP | 0 SUBROUTINE ERROR | - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 06063 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | X 00000 R Y 01750 R W 03720 R | JX 05670 I ZERO 05671 R FZERO 05722 R | COEF 05753 R ACC 06004 R BMPAR 06005 R | CPAR 06006 R M 06007 I N 06010 I INDER 06011 I INDX0 06012 I INDXK 06013 I | INDPOT 06014 I INDMN 06015 I INDCHA 06016 I | INDACC 06017 I MN2 06020 I MPTEST 06021 I | YMEAN 06022 R ITERA 06027 I MNTEST 06030 I | AM 06031 R AMPAR 06062 R - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | K 06064 I L1 06065 I L2 06066 I | I 06067 I Y2 06070 R J 06071 I | F 06072 R JSPAR 06073 I X1 06074 R Y1 06075 R X3 06076 R Y3 06077 R | IS 06100 I YS 06101 R SI 06102 R | SPAR 06103 R CL 06104 R | - ENTRY POINTS | 0 ERROR SECTION 3 - SUBROUTINES CALLED | 0 FAPP SECTION 4 .FWRD. SECTION 5 .UN06. SECTION 6 | .FFIL. SECTION 7 .FCNV. SECTION 8 SYSLOC SECTION 9 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 36 13A 06201 1 49A 06405 2 31A 06260 | 3 33A 06262 16 134A 06732 23 54A 06422 | 20 55A 06424 18 59A 06431 19 70A 06504 | 21 76A 06554 22 FORMAT 06131 30 90A 06602 | 29 86A 06575 28 133A 06731 24 98A 06617 25 102A 06624 26 110A 06646 27 114A 06653 | 40 118A 06663 35 128A 06722 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 06777.| 1 12/01/64 PAGE 47 | FAP STORAGE MAP 0 FUNCTION FAPP TYPE R | - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 06063 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | X 00000 R Y 01750 R W 03720 R JX 05670 I ZERO 05671 R FZERO 05722 R | COEF 05753 R ACC 06004 R BMPAR 06005 R | CPAR 06006 R M 06007 I N 06010 I | INDER 06011 I INDX0 06012 I INDXK 06013 I | INDPOT 06014 I INDMN 06015 I INDCHA 06016 I INDACC 06017 I MN2 06020 I MPTEST 06021 I | YMEAN 06022 R ITERA 06027 I MNTEST 06030 I | AM 06031 R AMPAR 06062 R | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE B 06064 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | F.0000 06116 R I 06117 I AT 06120 R | M1 06121 I II 06122 I IJ 06123 I AN 06124 R N1 06125 I Z 06126 R | - ENTRY POINTS | 0 FAPP SECTION 3 | - SUBROUTINES CALLED | 0 .FXEM. SECTION 4 SYSLOC SECTION 5 - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 4 3A 06152 5 21A 06237 6 8A 06173 | 7 16A 06225 8 57A 06400 9 58A 06405 | 10 32A 06275 19 39A 06325 11 60A 06410 12 47A 06347 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 06443.| 1 12/01/64 PAGE 48 | EQU STORAGE MAP | 0 SUBROUTINE EQUAT - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 06063 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | X 00000 R Y 01750 R W 03720 R | JX 05670 I ZERO 05671 R FZERO 05722 R COEF 05753 R ACC 06004 R BMPAR 06005 R | CPAR 06006 R M 06007 I N 06010 I | INDER 06011 I INDX0 06012 I INDXK 06013 I | INDPOT 06014 I INDMN 06015 I INDCHA 06016 I | INDACC 06017 I MN2 06020 I MPTEST 06021 I YMEAN 06022 R ITERA 06027 I MNTEST 06030 I | AM 06031 R AMPAR 06062 R CF 05753 R | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | CMAT 06064 R - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | IPOT 07245 I MNMAX 07246 I I 07247 I | F1 07250 R F2 07251 R Z 07252 R | J 07253 I F3 07254 R JM 07255 I JJ 07256 I | - ENTRY POINTS | 0 EQUAT SECTION 3 | - SUBROUTINES CALLED | 0 MATS SECTION 4 .FXEM. SECTION 5 SYSLOC SECTION 6 - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 11 62A 07555 1 42A 07473 14 12A 07346 | 15 15A 07351 16 19A 07366 2 40A 07467 | 6 30A 07436 4 26A 07420 5 29A 07424 7 34A 07445 8 60A 07551 9 55A 07540 | 10 63A 07557 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 07613.| 1 12/01/64 PAGE 49 | MATRS STORAGE MAP 0 SUBROUTINE MATS | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | N 00001 I M 00002 I MM 00003 I | I 00004 I II 00005 I J 00006 I R 00007 R B 00010 R JJ 00011 I | KK 00012 I | - ENTRY POINTS | 0 MATS SECTION 2 | - SUBROUTINES CALLED 0 .FWRD. SECTION 3 .UN06. SECTION 4 .FFIL. SECTION 5 | .FCNV. SECTION 6 SYSLOC SECTION 7 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 1965 8A 00150 1966 FORMAT 00031 1967 55A 00364 15 48A 00332 70 13A 00201 7 14A 00204 | 8 17A 00231 9 20A 00233 11 28A 00256 | 10 24A 00251 130 39A 00304 12 36A 00275 | 13 40A 00307 14 43A 00320 16 54A 00356 | 17 57A 00367 1 FORMAT 00050 1968 84A 00563 28 77A 00522 25 71A 00475 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00734.| 1 12/01/64 PAGE 50 | COMPILATION RECRD2 | 0 DIAGNOSTIC MESSAGES -1 SOURCE ERROR 127 LEVEL 1 - WARNING ONLY| TRANSFER INTO DO NEST TO IFN 75 FROM OUTSIDE ITS RANGE| 1 12/01/64 PAGE 51 | COMPILATION CHEC | 0 DIAGNOSTIC MESSAGES -1 SOURCE ERROR 127 LEVEL 1 - WARNING ONLY| TRANSFER INTO DO NEST TO IFN 82 FROM OUTSIDE ITS RANGE| 1 12/01/64 PAGE 52 | COMPILATION ERRROR | 0 DIAGNOSTIC MESSAGES -1 SOURCE ERROR 127 LEVEL 1 - WARNING ONLY| TRANSFER INTO DO NEST TO IFN 49 FROM OUTSIDE ITS RANGE 1 12/01/64 PAGE 53 | | 0$ENTRY RAT10720 | $DATA 0LEVEL=2 UNRECOGNIZABLE PARAMETER ON $IBJOB CD... NODECK IS IGNORED.| 1 IBLDR 12/01/64 PAGE 54 | * CROSS REFERENCE LOGIC * | 0 | 0DECK INPUT 03004 SIZE = 14077 REAL LAES 03005 SIZE = 00004 REF'D BY - RECRD1| REAL RECO 03011 SIZE = 00121 REF'D BY - RECRD2| REAL // DELETED REFERS TO // AT 71714| REAL ...... 17065| VIRT RECORD REFERS TO RECRD2 AT 20203 VIRT .FWRD. REFERS TO FRDD AT 33412| VIRT .FRDD. REFERS TO FRDD AT 33407| VIRT .FSLI. REFERS TO FSLDI AT 42006| VIRT .FSLO. REFERS TO FSLDO AT 42102| VIRT REC REFERS TO RECRD1 AT 17351 VIRT OPERA REFERS TO OPER AT 23037| VIRT CHECK REFERS TO CHEC AT 24005| VIRT COMPUT REFERS TO COMPU AT 25032| VIRT INTERP REFERS TO INTER AT 25414| VIRT FAPP REFERS TO FAP AT 26715 VIRT STABLE REFERS TO STABL AT 22376| VIRT .FXEM. REFERS TO FXEM AT 32713| VIRT .UN06. REFERS TO UN06 AT 40404| VIRT .FFIL. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762 VIRT .UN05. REFERS TO UN05 AT 40403| VIRT .FRTN. REFERS TO FIOH AT 35057| VIRT CC.1 REFERS TO .XCC. AT 32707| VIRT CC.2 REFERS TO .XCC. AT 32710| VIRT CC.3 REFERS TO .XCC. AT 32711 VIRT CC.4 REFERS TO .XCC. AT 32712| VIRT SYSLOC REFERS TO .JBCON AT 02652| 0DECK RECRD1 17103 SIZE = 00316| REAL LAES DELETED REFERS TO INPUT AT 03005| REAL REC 17351 REF'D BY - INPUT VIRT RECORD REFERS TO RECRD2 AT 20203| VIRT E.1 REFERS TO .ERAS. AT 32703| VIRT E.2 REFERS TO .ERAS. AT 32704| VIRT E.3 REFERS TO .ERAS. AT 32705| VIRT E.4 REFERS TO .ERAS. AT 32706 VIRT CC.1 REFERS TO .XCC. AT 32707| VIRT CC.2 REFERS TO .XCC. AT 32710| VIRT CC.3 REFERS TO .XCC. AT 32711| VIRT CC.4 REFERS TO .XCC. AT 32712| VIRT SYSLOC REFERS TO .JBCON AT 02652 0DECK RECRD2 17421 SIZE = 00654| REAL RECO DELETED REFERS TO INPUT AT 03011| REAL RECORD 20203 REF'D BY - RECRD1, INPUT | VIRT .XP2. REFERS TO FXP2 AT 41525| VIRT .FRDD. REFERS TO FRDD AT 33407 VIRT .FXEM. REFERS TO FXEM AT 32713| VIRT .UN05. REFERS TO UN05 AT 40403| VIRT .FRTN. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762| VIRT E.1 REFERS TO .ERAS. AT 32703 VIRT E.2 REFERS TO .ERAS. AT 32704| 1 IBLDR 12/01/64 PAGE 55 | * CROSS REFERENCE LOGIC * | 0 VIRT E.3 REFERS TO .ERAS. AT 32705| VIRT E.4 REFERS TO .ERAS. AT 32706 VIRT CC.1 REFERS TO .XCC. AT 32707| VIRT CC.2 REFERS TO .XCC. AT 32710| VIRT CC.3 REFERS TO .XCC. AT 32711| VIRT CC.4 REFERS TO .XCC. AT 32712| VIRT SYSLOC REFERS TO .JBCON AT 02652 0DECK STABL 20275 SIZE = 02215| REAL STABLE 22376 REF'D BY - INPUT | VIRT .FRDD. REFERS TO FRDD AT 33407| VIRT .FWRD. REFERS TO FRDD AT 33412| VIRT .FXEM. REFERS TO FXEM AT 32713 VIRT .UN05. REFERS TO UN05 AT 40403| VIRT .FRTN. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762| VIRT .UN06. REFERS TO UN06 AT 40404| VIRT .FFIL. REFERS TO FIOH AT 35057 VIRT SYSLOC REFERS TO .JBCON AT 02652| 0DECK OPER 22512 SIZE = 00455| REAL OPERA 23037 REF'D BY - INPUT | VIRT .XP3. REFERS TO FXP3 AT 41643| VIRT EXP REFERS TO FXPF AT 40611 VIRT ALOG REFERS TO FLOG AT 40406| VIRT ALOG10 REFERS TO FLOG AT 40405| VIRT ATAN REFERS TO FATN AT 41274| VIRT SIN REFERS TO FSCN AT 40733| VIRT COS REFERS TO FSCN AT 40732 VIRT TANH REFERS TO FTNH AT 41201| VIRT SQRT REFERS TO FSQR AT 41126| VIRT ARSIN REFERS TO FASC AT 42161| VIRT ARCOS REFERS TO FASC AT 42160| VIRT .FXEM. REFERS TO FXEM AT 32713 VIRT E.1 REFERS TO .ERAS. AT 32703| VIRT E.2 REFERS TO .ERAS. AT 32704| VIRT E.3 REFERS TO .ERAS. AT 32705| VIRT E.4 REFERS TO .ERAS. AT 32706| VIRT CC.1 REFERS TO .XCC. AT 32707 VIRT CC.2 REFERS TO .XCC. AT 32710| VIRT CC.3 REFERS TO .XCC. AT 32711| VIRT CC.4 REFERS TO .XCC. AT 32712| VIRT SYSLOC REFERS TO .JBCON AT 02652| 0DECK CHEC 23167 SIZE = 00632 REAL // DELETED REFERS TO // AT 71714| REAL CHECK 24005 REF'D BY - INPUT | VIRT SQRT REFERS TO FSQR AT 41126| VIRT .FWRD. REFERS TO FRDD AT 33412| VIRT .UN06. REFERS TO UN06 AT 40404 VIRT .FFIL. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762| VIRT SYSLOC REFERS TO .JBCON AT 02652| 0DECK COMPU 24021 SIZE = 01040| REAL // DELETED REFERS TO // AT 71714 1 IBLDR 12/01/64 PAGE 56 | * CROSS REFERENCE LOGIC * | 0 REAL COMPUT 25032 REF'D BY - INPUT | VIRT COS REFERS TO FSCN AT 40732| VIRT .FWRD. REFERS TO FRDD AT 33412 VIRT INTERP REFERS TO INTER AT 25414| VIRT EQUAT REFERS TO EQU AT 30443| VIRT ERROR REFERS TO ERRROR AT 26331| VIRT .FXEM. REFERS TO FXEM AT 32713| VIRT .UN06. REFERS TO UN06 AT 40404 VIRT .FFIL. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762| VIRT E.1 REFERS TO .ERAS. AT 32703| VIRT E.2 REFERS TO .ERAS. AT 32704| VIRT E.3 REFERS TO .ERAS. AT 32705 VIRT E.4 REFERS TO .ERAS. AT 32706| VIRT CC.1 REFERS TO .XCC. AT 32707| VIRT CC.2 REFERS TO .XCC. AT 32710| VIRT CC.3 REFERS TO .XCC. AT 32711| VIRT CC.4 REFERS TO .XCC. AT 32712 VIRT SYSLOC REFERS TO .JBCON AT 02652| 0DECK INTER 25061 SIZE = 00360| REAL // DELETED REFERS TO // AT 71714| REAL INTERP 25414 REF'D BY - COMPU , INPUT | VIRT .FWRD. REFERS TO FRDD AT 33412 VIRT .FXEM. REFERS TO FXEM AT 32713| VIRT .UN06. REFERS TO UN06 AT 40404| VIRT .FFIL. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762| VIRT SYSLOC REFERS TO .JBCON AT 02652 0DECK ERRROR 25441 SIZE = 00714| REAL // DELETED REFERS TO // AT 71714| REAL ERROR 26331 REF'D BY - COMPU | VIRT FAPP REFERS TO FAP AT 26715| VIRT .FWRD. REFERS TO FRDD AT 33412 VIRT .UN06. REFERS TO UN06 AT 40404| VIRT .FFIL. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762| VIRT SYSLOC REFERS TO .JBCON AT 02652| 0DECK FAP 26355 SIZE = 00360 REAL // DELETED REFERS TO // AT 71714| REAL FAPP 26715 REF'D BY - ERRROR, INPUT | VIRT .FXEM. REFERS TO FXEM AT 32713| VIRT SYSLOC REFERS TO .JBCON AT 02652| 0DECK EQU 26735 SIZE = 01530 REAL // DELETED REFERS TO // AT 71714| REAL EQUAT 30443 REF'D BY - COMPU | VIRT MATS REFERS TO MATRS AT 31337| VIRT .FXEM. REFERS TO FXEM AT 32713| VIRT SYSLOC REFERS TO .JBCON AT 02652 0DECK MATRS 30465 SIZE = 00734| REAL MATS 31337 REF'D BY - EQU | VIRT .FWRD. REFERS TO FRDD AT 33412| 1 IBLDR 12/01/64 PAGE 57 | * CROSS REFERENCE LOGIC * 0 VIRT .UN06. REFERS TO UN06 AT 40404| VIRT .FFIL. REFERS TO FIOH AT 35057| VIRT .FCNV. REFERS TO FIOH AT 34762| VIRT SYSLOC REFERS TO .JBCON AT 02652| 1 IBLDR 12/01/64 PAGE 58 * IOCS CONTROLS * | 0 | 0FILE LIST | LOC | 0 02750 GROUP RESERVED 00001 BUFFERS, MAXIMUM FILE OPEN COUNT = 00001 02751 UNIT05 | 0 02752 GROUP RESERVED 00001 BUFFERS, MAXIMUM FILE OPEN COUNT = 00001| 02753 UNIT06 | 1 IBLDR 12/01/64 PAGE 59 | * MEMORY MAP * 0 | 0SYSTEM 00000 THRU 02717| FILE BLOCK ORIGIN 02720| FILES 1. UNIT05 | 2. UNIT06 FILE LIST ORIGIN 02750| PRE-EXECUTION INITIALIZATION 02754| CALL ON OBJECT PROGRAM 02777| OBJECT PROGRAM 03004 THRU 46123| 0 DECK ORIGIN CONTROL SECTIONS (/NAME/=NON 0 LENGTH, (LOC)=DELETED, 'LOC'=MOVED, *=NOT REFERENCED) | 1. INPUT 03004 /LAES / 03005 /RECO / 03011 /// /(71714) ...... 17065 * | 2. RECRD1 17103 /LAES /(03005) REC 17351 | 3. RECRD2 17421 /RECO /(03011) RECORD 20203 | 4. STABL 20275 STABLE 22376 5. OPER 22512 OPERA 23037 | 6. CHEC 23167 /// /(71714) CHECK 24005 | 7. COMPU 24021 /// /(71714) COMPUT 25032 | 8. INTER 25061 /// /(71714) INTERP 25414 | 9. ERRROR 25441 /// /(71714) ERROR 26331 10. FAP 26355 /// /(71714) FAPP 26715 | 11. EQU 26735 /// /(71714) EQUAT 30443 | 12. MATRS 30465 MATS 31337 | 13. .LXCON 31421 .LXSTR 31421 * .LXSTP 31424 .LXOUT 31472 * .LXERR 31501 .LXCAL 31504 * | .LXRTN 31504 IBEXIT 31504 * .DBCLS 31666 * .LXARG 32035 * .L0 32060 * .CLSE 32066 .LFBL 32067 * .LUNB 32070 .DFOUT 32071 | 14. .IODEF 32075 .DEFIN 32075 .ATTAC 32101 * .CLOSE 32103 .OPEN 32105 .READ 32107 | .WRITE 32111 .BSR 32121 * .READR 32131 .RELES 32133 * .LAREA 32144 | .LFBLK 32162 .LTSX 32165 * .RLHLD 32175 * .AREA1 32200 .LUNBL 32206 | .ENTRY 32212 .GOA 32247 .GO 32253 .DERR 32267 .NOPXI 32270 .COMXI 32272 .EX34 32314 | 15. .LXSL 32321 .LXSLB 32321 * .LXSEL 32322 .LXSL1 32323 .LXTST 32326 * .LXOVL 32371 * | .LXMOD 32433 * .LXIND 32457 * .LXDIS 32462 * .LXFLG 32463 * .LTCH 32464 | 16. .FPTRP 32473 .FFPT. 32473 * .FPOUT 32622 * .FPARG 32630 * /.COUNT/ 32632 * OVFLOW 32676 * | 17. .ERAS. 32703 E.1 32703 E.2 32704 E.3 32705 E.4 32706 18. .XCC. 32707 CC.1 32707 CC.2 32710 CC.3 32711 CC.4 32712 | 19. FXEM 32713 .FXEM. 32713 TRACE 32721 * /.OPTW./ 33256 * | 20. FCOM 33272 .FCOM. 33272 /FORCOM/ 33375 * | 21. FRDD 33407 .FRDD. 33407 .FWRD. 33412 | 22. FIOS 33524 .FIOS. 33524 .FSEL. 33713 .FRTB. 33730 .FRTD. 33745 /.FBF1./ 34177 * 23. FIOH 34203 .FIOH. 34203 .FCNV. 34762 .FFIL. 35057 .FRTN. 35057 | 24. FCNV 35176 .FCON. 35176 | 25. FIOT 40264 .FIOT. 40264 .FIOE. 40355 | 26. UN05 40403 .UN05. 40403 | 27. UN06 40404 .UN06. 40404 ..UN06 40404 * 28. FLOG 40405 ALOG10 40405 ALOG 40406 | 29. FXPF 40611 EXP 40611 | 30. FSCN 40732 COS 40732 SIN 40733 | 31. FSQR 41126 SQRT 41126 | 32. FTNH 41201 TANH 41201 33. FATN 41273 ATAN2 41273 * ATAN 41274 | 34. FXP2 41525 .XP2. 41525 | 35. FXP3 41643 .XP3. 41643 | 36. FSLDI 41770 .FSLI. 42006 .FSDI. 42014 * | 1 IBLDR 12/01/64 PAGE 60 * MEMORY MAP * | 0 37. FSLI 42025 .SLI. 42025 .SLI1. 42032 .SDI. 42040 .SDI1. 42046 | 38. FSLDO 42064 .FSLO. 42102 .FSDO. 42110 * | 39. FSLO 42121 .SLO. 42121 .SLO2. 42127 .SDO. 42134 .SDO2. 42146 | 40. FASC 42160 ARCOS 42160 ARSIN 42161 41. .IOCS 42311 .L(0) 42311 .MONSW 42331 .TEOR 42400 .DEFI. 42460 .JOINX 42524 * | .CLOS. 42543 .ATTC. 42556 .SH1 42770 * .SH9 43032 * .OPEN. 43053 | .OP4 43101 * .OP7 43132 * .OP9.2 43146 * .RLSE. 43220 .RER2. 43220 | .READ. 43221 .RER1. 43244 .WRIT. 43246 .MNT1A 43436 * .EOFEX 43517 * | .FEEIT 43567 .GTIOX 43610 .RW7 43726 * .RE7 44345 * .ENDTR 45006 .SEL59 45010 * .BSR. 45427 .EOTOF 45554 .ETOF3 45562 * .SWITC 45611 | .TCHEX 46116 .BASIO 46121 * | 42. .IOCSM 46124 | 43. // 71714 | 0I/O BUFFERS 46124 THRU 71671 0UNUSED CORE 71672 THRU 71713| 0 1ARCTG(X)/X 7.7.65 0A=X*X 0F=ATAN(X)/X 0W=C1*X/ATAN(X) 0MAXDIVERGENCES 1.10E-06 -3.97E-06 5.59E-08 0ACCURACY 1.10E-06 0CONVERGENCE BUT NOT WITHIN THE ACCURACY 0APPROXIMATION INTERVAL 9.9999998E-17 1.0000000E 00 0NUMBER OF ITERATIONS 4 0POWER POLYNOMIALS 0ZEROES 4.0849488E-01 ERRORS 1.1020483E-01 1.0439701E-01 MEANERROR 1.0671264E-01 0DEGREES 0 0 0COEFFICIENTS 8.89795160E-01 1.00000000E 00 0MOVING PARAMETER HALVED TO 1.000E 00 AFTER 3 ITERATIONS 0MOVING PARAMETER HALVED TO 5.000E-01 AFTER 5 ITERATIONS 0MOVING PARAMETER HALVED TO 2.500E-01 AFTER 7 ITERATIONS 0MOVING PARAMETER HALVED TO 1.250E-01 AFTER 10 ITERATIONS 0MOVING PARAMETER HALVED TO 6.250E-02 AFTER 13 ITERATIONS 0MOVING PARAMETER HALVED TO 3.125E-02 AFTER 15 ITERATIONS 0MOVING PARAMETER HALVED TO 1.562E-02 AFTER 19 ITERATIONS 0MOVING PARAMETER HALVED TO 7.812E-03 AFTER 22 ITERATIONS 0MOVING PARAMETER HALVED TO 3.906E-03 AFTER 25 ITERATIONS 0MOVING PARAMETER HALVED TO 1.953E-03 AFTER 28 ITERATIONS 0MOVING PARAMETER HALVED TO 9.766E-04 AFTER 31 ITERATIONS 0MOVING PARAMETER HALVED TO 4.883E-04 AFTER 34 ITERATIONS 0MOVING PARAMETER HALVED TO 2.441E-04 AFTER 37 ITERATIONS 0MOVING PARAMETER HALVED TO 1.221E-04 AFTER 40 ITERATIONS 0MOVING PARAMETER HALVED TO 6.104E-05 AFTER 42 ITERATIONS 0MOVING PARAMETER HALVED TO 3.052E-05 AFTER 45 ITERATIONS 0MOVING PARAMETER HALVED TO 1.526E-05 AFTER 48 ITERATIONS 0MOVING PARAMETER HALVED TO 7.629E-06 AFTER 50 ITERATIONS 0MOVING PARAMETER HALVED TO 3.815E-06 AFTER 52 ITERATIONS 0MOVING PARAMETER HALVED TO 1.907E-06 AFTER 54 ITERATIONS 0MOVING PARAMETER HALVED TO 9.537E-07 AFTER 56 ITERATIONS 0MOVING PARAMETER HALVED TO 4.768E-07 AFTER 61 ITERATIONS 0MOVING PARAMETER HALVED TO 2.384E-07 AFTER 63 ITERATIONS 0MOVING PARAMETER HALVED TO 1.192E-07 AFTER 65 ITERATIONS 0MOVING PARAMETER HALVED TO 5.960E-08 AFTER 67 ITERATIONS 0MOVING PARAMETER HALVED TO 2.980E-08 AFTER 70 ITERATIONS 0MOVING PARAMETER HALVED TO 1.490E-08 AFTER 73 ITERATIONS 0MOVING PARAMETER HALVED TO 7.451E-09 AFTER 75 ITERATIONS 0MOVING PARAMETER HALVED TO 3.725E-09 AFTER 78 ITERATIONS 0MOVING PARAMETER HALVED TO 1.863E-09 AFTER 80 ITERATIONS 0MOVING PARAMETER HALVED TO 9.313E-10 AFTER 82 ITERATIONS 0MOVING PARAMETER HALVED TO 4.657E-10 AFTER 85 ITERATIONS 0MOVING PARAMETER HALVED TO 2.328E-10 AFTER 88 ITERATIONS 0MOVING PARAMETER HALVED TO 1.164E-10 AFTER 91 ITERATIONS 0MOVING PARAMETER HALVED TO 5.821E-11 AFTER 94 ITERATIONS 0MOVING PARAMETER HALVED TO 2.910E-11 AFTER 97 ITERATIONS 0MOVING PARAMETER HALVED TO 1.455E-11 AFTER 99 ITERATIONS 0MOVING PARAMETER HALVED TO 7.276E-12 AFTER 101 ITERATIONS 0NO CONVERGENCE 0APPROXIMATION INTERVAL 9.9999998E-17 1.0000000E 00 0NUMBER OF ITERATIONS 101 0POWER POLYNOMIALS 0ZEROES 1.5718377E-01 6.6636829E-01 ERRORS 1.2226941E 00 1.0402451E 00 9.3663055E-01 MEANERROR 1.0302590E 00 0DEGREES 1 0 0COEFFICIENTS 9.87056065E-01 1.00000000E 00 0NO PRINT ALLOWED 1 2000 LINES OUTPUT. $IBSYS $STOP | 1 PERIPHERAL UNIT POSITIONS AT END OF JOBS | 1 SYSPP1 IS B2 REC. 00391, FILE 00000 SYSOU1 IS B1 REC. 00465, FILE 00000 SYSIN1 IS A2 REC. 00002, FILE 00001| 1 END OF JOBS | 1