$EXECUTE FORTRAN 1 * ID MPR3 1 * XEQ * LIST 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 1 C MPR30002 C L1 = CLASS OF A (OPERAND) MPR30003 C L2 = INDEX OF A MPR30004 C C1 = CONSTANT TO BE ADDED TO A MPR30005 C L3 = OPERATION CODE MPR30006 C L4 = NO. OF TIMES OPERATION IS TO BE PERFORMED MPR30007 C L5 = CLASS OF B (OPERATOR) MPR30008 C L6 = INDEX OF B MPR30009 C C2 = CONSTANT TO BE ADDED TO B MPR30010 C L7 = CLASS OF RESULT (C) MPR30011 C L8 = INDEX OF RESULT MPR30012 C MPR30013 C Z = INPUT VARIABLE (130) MPR30014 C Y = REGRESSION DEPENDENT VARIABLE (20) MPR30015 C X = REGRESSION INDEPENDENT VARIABLE (130) MPR30016 C W = INTERMEDIATE VARIABLE (70) MPR30017 C MPR30018 C 1FWT =1,THEN ALL WHTS =1.0 MPR30019 C 1FSTEP =1, DO NOT PRINT EACH STEP MPR30020 C 1FRAW =1, DO NO PRINT RAW SUMS AND SQUARES MPR30021 C 1FAVE =1, DO NOT PRINT AVERAGES MPR30022 C 1FRESD =1, DO NOT PRINT RESIDUAL SUMS OF SQUARES MPR30023 C 1FCOEN =1, DO NOT PRINT PARTIAL COEFFICIENTS MPR30024 C 1FPRED =1, DO NOT CALC PREDICTED VALUES MPR30025 C 1FCNST =1, DO NOT HAVE CNST TERM IN EQUATION MPR30026 C IFTRAN = 1, DO NOT PRINT TABLE OF VARIABLE TRANSFORMATIONS MPR30027 C MPR30028 C INPUT DATA SPECIFICATIONS MPR30029 C MAX. OF 25 OPERATION COMMAND CARDS EXCLUDING 5 CONTROL AT END MPR30030 C MPR30031 C MPR30032 EQUIVALENCE (DATA,X) MPR30033 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30034 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30035 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30036 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30037 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1,NAMVAR,NAMEOP MPR30038 5, NODATA MPR30039 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30040 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30041 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30042 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30043 4,NAMEOP(8),NAMVAR(4), FMT(12) MPR30044 C MPR30045 C MPR30046 CALL PERM(6H Z , 6H W , 6H X , 6H Y , 6H + , MPR30047 1 6H - , 6H * , 6H / , 6HLOGE , 6HE ** , 6HARCTAN, MPR30048 2 6H** , NAMVAR(1),NAMVAR(2),NAMVAR(3),NAMVAR(4),NAMEOP(1), MPR30049 3NAMEOP(2),NAMEOP(3),NAMEOP(4),NAMEOP(5),NAMEOP(6),NAMEOP(7), MPR30050 4NAMEOP(8) ) MPR30051 READ INPUT TAPE5, 5,TOL,EFIN,EFOUT,NOPROB,INVAR,NODATA,IFWT, MPR30052 1IFSTEP,IFRAW,IFAVE, MPR30053 2IFRESD,IFCOEN,IFPRED,IFCNST,N,IFTRAN, IFTVAR,IFFORM MPR30054 REWIND 3 MPR30055 NOIN = 0 MPR30056 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 2 VAR = 0 MPR30057 K = 0 MPR30058 FLEVEL =0 MPR30059 NOENT = 0 MPR30060 NOMIN = 0 MPR30061 NOMAX = 0 MPR30062 NOVAR = INVAR MPR30063 NVP1 = NOVAR+1 MPR30064 NIX = 1 MPR30065 IF (INVAR - 130) 100, 100, 1800 MPR30066 100 CALL SETUP MPR30067 IF (IFTRAN) 900, 102, 105 MPR30068 102 CALL PRTTRF MPR30069 105 DO 110 I = 1,NVP1 MPR30070 DO 110 J = 1,NVP1 MPR30071 110 VECTOR (I,J) = 0.0 MPR30072 IF (IFWT) 900,114,115 MPR30073 114 ASSIGN 130 TO IFWGT MPR30074 GO TO 120 MPR30075 115 ASSIGN 140 TO IFWGT MPR30076 120 IF (IFFORM) 900, 124,122 MPR30077 122 READ INPUT TAPE 5, 125, (FMT (I), I = 1,12) MPR30078 ASSIGN 127 TO IFMT MPR30079 GO TO 126 MPR30080 124 ASSIGN 128 TO IFMT MPR30081 126 DO 150 K = 1, NODATA MPR30082 GO TO IFMT, (127, 128) MPR30083 127 READ INPUT TAPE 5, FMT, RUN, (Z(N1), N1 = 1,N), WHT MPR30084 GO TO 129 MPR30085 128 READ INPUT TAPE 5, 1000, RUN, (Z(N1), N1 = 1,N) ,WHT MPR30086 129 CALL SWITCH MPR30087 GO TO IFWGT, (130,140) MPR30088 130 CALL WGTSUM MPR30089 GO TO 150 MPR30090 140 CALL SUM MPR30091 150 CONTINUE MPR30092 C MPR30093 C HAVE CALCULATED SUMS OF SQUARES AND CROSS PRODUCTS - PRINT RESULTSMPR30094 C IF 1FRAW SAY O.K. MPR30095 565 NOVMI = NOVAR - 1 MPR30096 566 NOVPL = NOVAR + 1 MPR30097 567 WRITE OUTPUT TAPE 6, 90, NOPROB,NODATA, NOVAR, VECTOR(NOVPL, MPR30098 1NOVPL),EFIN, EFOUT MPR30099 568 PRINT 90,NOPROB,NODATA,NOVAR,VECTOR(NOVPL,NOVPL), MPR30100 1EFIN, EFOUT MPR30101 570 IF (IFRAW) 900, 580, 650 MPR30102 580 WRITE OUTPUT TAPE 6, 15 MPR30103 590 WRITE OUTPUT TAPE 6, 20,(I,VECTOR(I,NOVPL),I=1,NOVMI) MPR30104 600 WRITE OUTPUT TAPE 6, 25, VECTOR (NOVAR, NOVPL) MPR30105 610 WRITE OUTPUT TAPE 6, 30 MPR30106 620 WRITE OUTPUT TAPE 6, 35,((I,J,VECTOR(I,J),J=I,NOVMI),I=1,NOVMI) MPR30107 630 WRITE OUTPUT TAPE 6, 40,(I,VECTOR(I,NOVAR),I=1,NOVMI) MPR30108 640 WRITE OUTPUT TAPE 6, 45, VECTOR(NOVAR,NOVAR) MPR30109 C MPR30110 C CALCULATE RESIDUAL SUMS OF SQUARES AND CROSS PRODUCTS MPR30111 650 CALL RESID MPR30112 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 3 CALL NORM MPR30113 CALL CALC MPR30114 C MPR30115 C PRINT FINAL RESULTS AND RETURN MPR30116 1580 WRITE OUTPUT TAPE 6, 1586,(L,VECTOR(L,L),L=1,NOVMI ) MPR30117 1581 IF( IFPRED) 900,1582,1900 MPR30118 1582 REWIND 3 MPR30119 1583 WRITE OUTPUT TAPE 6, 85 MPR30120 1590 DO 1660 N = 1, NODATA MPR30121 1600 READ TAPE 3,(DATA(L),L= 1, NOVAR), RUN MPR30122 1610 YPRED = CNST MPR30123 1620 DO 1630 I = 1,NOIN MPR30124 1640 K = INDEX(I) MPR30125 1630 YPRED= YPRED + COEN(I) * DATA(K) MPR30126 1650 DEV = DATA(NOVAR) - YPRED MPR30127 1660 WRITE OUTPUT TAPE 6, 80, RUN, DATA(NOVAR), YPRED, DEV MPR30128 C MPR30129 C TRANSFORM RESULTS TO FORM OF ORIGINAL INPUT MPR30130 1900 IF (IFTVAR) 900,1905, 910 MPR30131 1905 CALL PRINTX MPR30132 CALL OUTPUT MPR30133 C MPR30134 910 REWIND 3 MPR30135 121 CALL EXIT MPR30136 900 PRINT 905 MPR30137 GO TO 910 MPR30138 1800 WRITE OUTPUT TAPE 6, 93, INVAR MPR30139 GO TO 910 MPR30140 C MPR30141 5 FORMAT (3F10.5,3I5,1H 8I2, I4, 3I2) MPR30142 10 FORMAT (6(F12.5)) MPR30143 15 FORMAT (1H 49H SUM OF VARIABLES//)MPR30144 20 FORMAT (1H 11H SUM X( I2,3H) = F12.4,8H SUM X( I2,3H) =F12.4,MPR30145 18H SUM X( I2,3H) =F12.4,8H SUM X(I2,3H) =F12.4 ) MPR30146 25 FORMAT (17H SUM Y =F12.4) MPR30147 30 FORMAT(1H0 70H RAW SUM OF SQUARES AMPR30148 1ND CROSS PRODUCTS// ) MPR30149 35 FORMAT (1H 7H X(I2,7H) VS X(I2,3H) = F15.6, MPR30150 1 6H X(I2,7H) VS X(I2,3H) = F15.6, MPR30151 2 6H X(I2,7H) VS X(I2,3H) = F15.6 ) MPR30152 40 FORMAT (1H 7H X(I2,12H) VS Y =F15.6, MPR30153 1 6H X(I2,12H) VS Y =F15.6, MPR30154 2 6H X(I2,12H) VS Y =F15.6 ) MPR30155 45 FORMAT (1H 21H Y VS Y =F15.6) MPR30156 80 FORMAT (8H F6.0,2H F12.5,3H F12.5,2H F12.5) MPR30157 85 FORMAT (1H043H PREDICTED VS ACTUAL RESULTS /59H MPR30158 1 RUN NO. ACTUAL PREDICTED DEVIATION //) MPR30159 90 FORMAT (22H1STEPWISE REGRESSION //12H PROBLEM NO I10 //13H NO OF MPR30160 1DATA = I5 //18H NO OF VARIABLES = I10 //30H WEIGHTED DEGREES OF FRMPR30161 2EEDOM = F12.2 //28H F LEVEL TO ENTER VARIABLE = F10.3 //29H F LEVEMPR30162 3L TO REMOVE VARIABLE = F9.3 /// ) MPR30163 93 FORMAT (45H1INVAR SHOULD NOT BE GREATER THAN 150. IT IS I5 ) MPR30164 905 FORMAT (42H ERROR IN CONTROL CARD, PROBLEM TERMINATED) MPR30165 1000 FORMAT (6F12.6) MPR30166 1586 FORMAT (24H0 DIAGONAL ELEMENTS //20H VAR.NO. VALUE// MPR30167 1(1H I 7, F16.6)) MPR30168 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 4 125 FORMAT (12A6) MPR30169 END(1,1,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 5 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 1006 01756 14076 33374 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NODATA 14077 33375 NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 | NOSTEP 31539 75463 NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 | VECTOR 31537 75461 WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 | Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT A 1005 01755 B 995 01743 C 985 01731 CTEMP 920 01630 FMT 910 01616 | 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | DEV 898 01602 IFFORM 897 01601 IFMT 896 01600 IFPRED 895 01577 IFRAW 894 01576 | IFTRAN 893 01575 IFTVAR 892 01574 IFWGT 891 01573 IFWT 890 01572 INVAR 889 01571 I 888 01570 J 887 01567 L 886 01566 NOPROB 885 01565 N 884 01564 | YPRED 883 01563 | 0 SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT STATEMENTS | 0 EFN LOC EFN LOC EFN LOC EFN LOC EFN LOC | 8)5 5 01546 8)A 10 01541 8)F 15 01537 8)K 20 01525 8)P 25 01502 8)U 30 01475 8)13 35 01457 8)18 40 01422 8)1D 45 01366 8)2G 80 01360 | 8)2L 85 01350 8)2Q 90 01324 8)2T 93 01256 8)3T 125 01215 8)S9 905 01244 | 8)V8 1000 01234 8)1HI 1586 01232 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT 1) 871 01547 2) 614 01146 3) 645 01205 4) 32767 77777 5) 610 01142 | 6) 646 01206 A)101 571 01073 A)103 584 01110 A)105 597 01125 B) 621 01155 | C)G4 875 01553 C)G5 876 01554 C)101 877 01555 C)103 878 01556 C)105 879 01557 | C)200 880 01560 C)203 881 01561 C)206 882 01562 D)215 446 00676 D)219 482 00742 | D)419 484 00744 D)41F 533 01025 D)41I 549 01045 D)61F 532 01024 D)81H 540 01034 D)81I 551 01047 D)C1I 548 01044 D)E1I 547 01043 E)2 160 00240 E)1J 558 01056 | E)23 166 00246 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 6 | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | CALC 15 00017 EXIT 20 00024 NORM 14 00016 OUTPUT 19 00023 PERM 1 00001 | PRINTX 18 00022 PRTTRF 6 00006 RESID 13 00015 SETUP 5 00005 SUM 9 00011 | SWITCH 7 00007 WGTSUM 8 00010 (FIL) 11 00013 (FPT) 0 00000 (RLR) 17 00021 | (RTN) 3 00003 (RWT) 4 00004 (SPH) 12 00014 (STH) 10 00012 (TSB) 16 00020 (TSH) 2 00002 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | | CALC EXIT NORM OUTPUT PERM PRINTX PRTTRF RESID SETUP SUM | SWITCH WGTSUM (FIL) (FPT) (RLR) (RTN) (RWT) (SPH) (STH) (TSB) (TSH) | 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 100 37 00231 102 39 00242 105 40 00250 110 42 00263 114 44 00276 | 115 46 00301 120 47 00303 122 48 00307 124 55 00323 126 56 00325 127 58 00331 128 65 00355 129 71 00400 130 73 00404 140 75 00410 | 150 76 00413 565 77 00416 566 78 00421 567 79 00430 568 81 00453 | 570 83 00475 580 84 00501 590 85 00507 600 90 00545 610 92 00557 | 620 93 00565 630 101 00634 640 106 00666 650 108 00677 1580 111 00702 | 1581 116 00724 1582 117 00730 1583 118 00732 1590 119 00736 1600 120 00745 1610 125 00761 1620 126 00763 1640 127 00766 1630 128 00772 1650 129 01000 | 1660 130 01003 1900 132 01026 1905 133 01032 910 135 01035 121 136 01040 | 900 137 01050 1800 139 01060 | 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 7 | 00000 (FPT) BCD 1(FPT) 00070 TSX NAMEOP 00160 LXD 6)+4,4 | 00001 PERM BCD 1PERM 00071 TSX NAMEOP-1 00161 26A CLA 2)+2 | 00002 (TSH) BCD 1(TSH) 00072 TSX NAMEOP-2 00162 STO NOIN | 00003 (RTN) BCD 1(RTN) 00073 TSX NAMEOP-3 00163 27A CLA 2)+2 | 00004 (RWT) BCD 1(RWT) 00074 TSX NAMEOP-4 00164 LRS 18 00005 SETUP BCD 1SETUP 00075 TSX NAMEOP-5 00165 ORA 6) | 00006 PRTTRF BCD 1PRTTRF 00076 TSX NAMEOP-6 00166 FAD 6) | 00007 SWITCH BCD 1SWITCH 00077 TSX NAMEOP-7 00167 STO VAR | 00010 WGTSUM BCD 1WGTSUM 00100 LXD 6)+4,4 00170 28A CLA 2)+2 | 00011 SUM BCD 1SUM 00101 23A CAL 2) 00171 STO K 00012 (STH) BCD 1(STH) 00102 SXD 6)+4,4 00172 LXD K,2 | 00013 (FIL) BCD 1(FIL) 00103 TSX (TSH),4 00173 29A CLA 2)+2 | 00014 (SPH) BCD 1(SPH) 00104 PZE 8)5 00174 LRS 18 | 00015 RESID BCD 1RESID 00105 LXD 6)+4,4 00175 ORA 6) | 00016 NORM BCD 1NORM 00106 24A STR 00176 FAD 6) 00017 CALC BCD 1CALC 00107 STQ TOL 00177 STO FLEVEL | 00020 (TSB) BCD 1(TSB) 00110 STR 00200 30A CLA 2)+2 | 00021 (RLR) BCD 1(RLR) 00111 STQ EFIN 00201 STO NOENT | 00022 PRINTX BCD 1PRINTX 00112 STR 00202 31A CLA 2)+2 | 00023 OUTPUT BCD 1OUTPUT 00113 STQ EFOUT 00203 STO NOMIN 00024 EXIT BCD 1EXIT 00114 STR 00204 32A CLA 2)+2 | 00025 $$ CLA 2 00115 STQ NOPROB 00205 STO NOMAX | 00026 STO 6)+5 00116 STR 00206 33A CLA INVAR | 00027 CLA (FPT) 00117 STQ INVAR 00207 STO NOVAR | 00030 STO 8 00120 STR 00210 SXD 6)+4,4 00031 STZ 4)-205 00121 STQ NODATA 00211 TSX A)103,4 | 00032 TSX A)101,4 00122 STR 00212 LXD C)103,4 | 00033 LXD C)101,1 00123 STQ IFWT 00213 SXD 6)+4,4 | 00034 TSX A)103,4 00124 STR 00214 TSX A)105,4 | 00035 LXD C)103,4 00125 STQ IFSTEP 00215 LXD C)105,1 00036 SXD 6)+4,4 00126 STR 00216 LXD 6)+4,4 | 00037 TSX A)105,4 00127 STQ IFRAW 00217 LXD NOVAR,1 | 00040 LXD C)105,2 00130 STR 00220 34A CLA NOVAR | 00041 LXD 6)+4,4 00131 STQ IFAVE 00221 ADD 2)+3 | 00042 LXD K,4 00132 STR 00222 STO NVP1 00043 SXD C)G4,4 00133 STQ IFRESD 00223 35A CLA 2)+3 | 00044 LXD NOVAR,4 00134 STR 00224 STO NIX | 00045 SXD C)G5,4 00135 STQ IFCOEN 00225 36A CLA INVAR | 00046 22A SXD 6)+4,4 00136 STR 00226 SUB 2)+4 | 00047 TSX PERM,4 00137 STQ IFPRED 00227 36A1 TZE 37A 00050 TSX B) 00140 STR 00230 TPL E)1J | 00051 TSX B)+2 00141 STQ IFCNST 00231 37A SXD 6)+4,4 | 00052 TSX B)+4 00142 STR 00232 TSX SETUP,4 | 00053 TSX B)+6 00143 STQ N 00233 LXD 6)+4,4 | 00054 TSX B)+8 00144 STR 00234 38A CLA IFTRAN 00055 TSX B)+10 00145 STQ IFTRAN 00235 38A1 TZE E)2 | 00056 TSX B)+12 00146 STR 00236 TPL E)23 | 00057 TSX B)+14 00147 STQ IFTVAR 00237 TRA 137A | 00060 TSX B)+16 00150 STR 00240 E)2 SXD C)G4,2 | 00061 TSX B)+18 00151 STQ IFFORM 00241 SXD C)G5,1 00062 TSX B)+20 00152 SXD 6)+4,4 00242 39A SXD 6)+4,4 | 00063 TSX B)+22 00153 TSX (RTN),4 00243 TSX PRTTRF,4 | 00064 TSX NAMVAR 00154 LXD 6)+4,4 00244 LXD 6)+4,4 | 00065 TSX NAMVAR-1 00155 25A CAL 2)+1 00245 TRA 40A | 00066 TSX NAMVAR-2 00156 SXD 6)+4,4 00246 E)23 SXD C)G4,2 00067 TSX NAMVAR-3 00157 TSX (RWT),4 00247 SXD C)G5,1 | 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 8 | | 00250 40A LXD 2)+3,4 00336 59A STR 00425 LXD C)101,2 | 00251 CLA NVP1 00337 STQ RUN 00426 TSX A)103,4 00252 STD 42A4 00340 60A LXD 2)+3,2 00427 LXD C)103,4 | 00253 LDQ NVP1 00341 CLA N 00430 79A CAL 2)+5 | 00254 MPY 2)+6 00342 STD 61A2 00431 SXD 6)+4,4 | 00255 ALS 17 00343 61A STR 00432 TSX (STH),4 | 00256 STD 42A2 00344 STQ Z+1,2 00433 PZE 8)2Q 00257 LDQ NVP1 00345 61A1 TXI *+1,2,1 00434 LXD 6)+4,4 | 00260 MPY 2)+6 00346 61A2 TXL 61A,2 00435 80A LDQ NOPROB | 00261 ALS 17 00347 63A STR 00436 STR | 00262 STD 42A3 00350 STQ WHT 00437 LDQ NODATA | 41A BSS 00351 SXD 6)+4,4 00440 STR 00263 42A CLA 3) 00352 TSX (RTN),4 00441 LDQ NOVAR | 00264 STO VECTOR+1,4 00353 LXD 6)+4,4 00442 STR | 00265 42A1 TXI *+1,4,131 00354 64A TRA 71A 00443 LDQ VECTOR+1,2 | 00266 42A2 TXL 42A,4 00355 65A CAL 2) 00444 STR | 00267 42A3 TIX *+1,4 00356 SXD 6)+4,4 00445 LDQ EFIN 00270 TXI *+1,4,1 00357 TSX (TSH),4 00446 STR | 00271 42A4 TXL 41A,4 00360 PZE 8)V8 00447 LDQ EFOUT | 00272 43A CLA IFWT 00361 LXD 6)+4,4 00450 STR | 00273 43A1 TZE 44A 00362 66A STR 00451 SXD 6)+4,4 | 00274 TPL 46A 00363 STQ RUN 00452 TSX (FIL),4 00275 TRA D)81I 00364 67A LXD 2)+3,2 81A BSS | 00276 44A CLA 5) 00365 CLA N 00453 TSX (SPH),4 | 00277 STO IFWGT 00366 STD 68A2 00454 PZE 8)2Q | 00300 45A TRA 47A 00367 68A STR 00455 LXD 6)+4,4 | 00301 46A CLA 5)+1 00370 STQ Z+1,2 00456 82A LDQ NOPROB 00302 STO IFWGT 00371 68A1 TXI *+1,2,1 00457 STR | 00303 47A CLA IFFORM 00372 68A2 TXL 68A,2 00460 LDQ NODATA | 00304 47A1 TZE 55A 00373 70A STR 00461 STR | 00305 TPL 48A 00374 STQ WHT 00462 LDQ NOVAR | 00306 TRA D)81I 00375 SXD 6)+4,4 00463 STR 00307 48A CAL 2) 00376 TSX (RTN),4 00464 LDQ VECTOR+1,2 | 00310 TSX (TSH),4 00377 LXD 6)+4,4 00465 STR | 00311 PZE 8)3T 00400 71A SXD 6)+4,4 00466 LDQ EFIN | 00312 49A LXD 2)+3,4 00401 TSX SWITCH,4 00467 STR | 00313 50A STR 00402 LXD 6)+4,4 00470 LDQ EFOUT 00314 STQ FMT+1,4 00403 72A TRA IFWGT 00471 STR | 00315 50A1 TXI *+1,4,1 00404 73A SXD 6)+4,4 00472 SXD 6)+4,4 | 00316 50A2 TXL 50A,4,12 00405 TSX WGTSUM,4 00473 TSX (FIL),4 | 52A BSS 00406 LXD 6)+4,4 00474 LXD 6)+4,4 | 00317 TSX (RTN),4 00407 74A TRA 76A 00475 83A CLA IFRAW 00320 53A CLA 5)+2 00410 75A SXD 6)+4,4 00476 83A1 TZE 84A | 00321 STO IFMT 00411 TSX SUM,4 00477 TPL D)215 | 00322 54A TRA 56A 00412 LXD 6)+4,4 00500 TRA D)41I | 00323 55A CLA 5)+3 76A BSS 00501 84A CAL 2)+5 | 00324 STO IFMT 00413 76A1 TXI *+1,4,1 00502 SXD 6)+4,4 00325 56A LXD 2)+3,4 00414 SXD C)200,4 00503 TSX (STH),4 | 00326 CLA NODATA 00415 76A2 TXL 57A,4 00504 PZE 8)F | 00327 STD 76A2 00416 77A CLA NOVAR 00505 TSX (FIL),4 | 00330 57A TRA IFMT 00417 SUB 2)+3 00506 LXD 6)+4,4 | 00331 58A CAL 2) 00420 STO NOVMI 00507 85A CAL 2)+5 00332 SXD 6)+4,4 00421 78A CLA NOVAR 00510 SXD 6)+4,4 | 00333 TSX (TSH),4 00422 ADD 2)+3 00511 TSX (STH),4 | 00334 PZE FMT 00423 STO NOVPL 00512 PZE 8)K | 00335 LXD 6)+4,4 00424 TSX A)101,4 00513 LXD 6)+4,4 | 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 9 | 00514 86A CLA 6)+3 00603 MPY 2)+6 00670 PZE 8)1D | 00515 STO 1)+3 00604 ALS 17 00671 LXD C)105,1 | 00516 LDQ NOVPL 00605 SUB 2)+6 00672 107A LDQ VECTOR+1,1 | 00517 MPY 2)+6 00606 ADD 1)+3 00673 STR 00520 ALS 17 00607 STO 1)+3 00674 TSX (FIL),4 | 00521 SUB 2)+6 00610 LXD 1)+3,2 00675 TRA 108A | 00522 ADD 1)+3 00611 LXD I,4 00676 D)215 LXD C)105,1 | 00523 STO 1)+3 00612 CLA NOVMI 108A BSS | 00524 LXD 1)+3,1 00613 STD 97A2 00677 TSX RESID,4 00525 LXD 2)+3,2 00614 PXD 0,4 109A BSS | 00526 CLA NOVMI 00615 STO J 00700 TSX NORM,4 | 00527 STD 87A2 00616 97A LDQ I 110A BSS | 00530 PXD 0,2 00617 STR 00701 TSX CALC,4 | 00531 STO I 00620 LDQ J 00702 111A CAL 2)+5 00532 87A LDQ I 00621 STR 00703 TSX (STH),4 | 00533 STR 00622 LDQ VECTOR+1,2 00704 PZE 8)1HI | 00534 LDQ VECTOR+1,1 00623 STR 00705 112A LXD 2)+3,2 | 00535 STR 00624 97A1 TXI *+1,2,131 00706 LXD 2)+3,4 | 00536 87A1 TXI *+1,1,1 00625 TXI *+1,4,1 00707 CLA NOVMI 00537 TXI *+1,2,1 00626 SXD J,4 00710 STD 113A2 | 00540 SXD I,2 00627 97A2 TXL 97A,4 00711 PXD 0,4 | 00541 87A2 TXL 87A,2 98A BSS 00712 STO L | 00542 89A SXD 6)+4,4 00630 98A1 TXI *+1,1,1 00713 113A LDQ L | 00543 TSX (FIL),4 00631 SXD I,1 00714 STR 00544 LXD 6)+4,4 00632 98A2 TXL 95A,1 00715 LDQ VECTOR+1,2 | 00545 90A CAL 2)+5 100A BSS 00716 STR | 00546 SXD 6)+4,4 00633 TSX (FIL),4 00717 113A1 TXI *+1,2,132 | 00547 TSX (STH),4 00634 101A CAL 2)+5 00720 TXI *+1,4,1 | 00550 PZE 8)P 00635 TSX (STH),4 00721 SXD L,4 00551 LXD 6)+4,4 00636 PZE 8)18 00722 113A2 TXL 113A,4 | 00552 91A LDQ VECTOR+1,4 00637 102A CLA 6)+3 115A BSS | 00553 STR 00640 STO 1)+3 00723 TSX (FIL),4 | 00554 SXD 6)+4,4 00641 LDQ NOVAR 00724 116A CLA IFPRED | 00555 TSX (FIL),4 00642 MPY 2)+6 00725 116A1 TZE 117A 00556 LXD 6)+4,4 00643 ALS 17 00726 TPL D)61F | 00557 92A CAL 2)+5 00644 SUB 2)+6 00727 TRA D)E1I | 00560 SXD 6)+4,4 00645 ADD 1)+3 00730 117A CAL 2)+1 | 00561 TSX (STH),4 00646 STO 1)+3 00731 TSX (RWT),4 | 00562 PZE 8)U 00647 LXD 1)+3,1 00732 118A CAL 2)+5 00563 TSX (FIL),4 00650 LXD 2)+3,2 00733 TSX (STH),4 | 00564 LXD 6)+4,4 00651 CLA NOVMI 00734 PZE 8)2L | 00565 93A CAL 2)+5 00652 STD 103A2 00735 TSX (FIL),4 | 00566 SXD 6)+4,4 00653 PXD 0,2 00736 119A LXD 2)+3,2 | 00567 TSX (STH),4 00654 STO I 00737 SXD C)206,2 00570 PZE 8)13 00655 103A LDQ I 00740 CLA NODATA | 00571 LXD 6)+4,4 00656 STR 00741 STD 131A2 | 00572 94A LXD 2)+3,1 00657 LDQ VECTOR+1,1 00742 D)219 LXD C)G5,1 | 00573 SXD C)203,1 00660 STR 00743 TRA 120A | 00574 CLA NOVMI 00661 103A1 TXI *+1,1,1 00744 D)419 LXD C)206,2 00575 STD 98A2 00662 TXI *+1,2,1 00745 120A CAL 2)+1 | 00576 PXD 0,1 00663 SXD I,2 00746 TSX (TSB),4 | 00577 STO I 00664 103A2 TXL 103A,2 00747 121A LXD 2)+3,4 | 95A BSS 105A BSS 00750 CLA NOVAR | 00600 96A CLA I 00665 TSX (FIL),4 00751 STD 122A2 00601 STO 1)+3 00666 106A CAL 2)+5 00752 122A STR | 00602 LDQ I 00667 TSX (STH),4 00753 STQ DATA+1,4 | 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 10 | | 00754 122A1 TXI *+1,4,1 01041 LXD 6)+4,4 01131 ADD 1)+3 00755 122A2 TXL 122A,4 01042 TRA 137A 01132 STO 1)+3 | 00756 124A STR 01043 D)E1I LXD C)G5,1 01133 LDQ NOVAR | 00757 STQ RUN 01044 D)C1I LXD C)103,4 01134 MPY 2)+6 | 00760 TSX (RLR),4 01045 D)41I LXD C)G4,2 01135 ALS 17 | 00761 125A CLA CNST 01046 TRA 137A 01136 SUB 2)+6 00762 STO YPRED 01047 D)81I LXD C)103,4 01137 ADD 1)+3 | 00763 126A LXD 2)+3,4 01050 137A SXD 6)+4,4 01140 STO C)105 | 00764 CLA NOIN 01051 TSX (SPH),4 01141 TRA 1,4 | 00765 STD 128A2 01052 PZE 8)S9 01142 5) TRA 73A | 00766 127A CLA INDEX+1,4 01053 TSX (FIL),4 01143 TRA 75A 00767 STO K 01054 LXD 6)+4,4 01144 TRA 58A | 00770 LXD K,2 01055 138A TRA 135A 01145 TRA 65A | 00771 SXD C)G4,2 01056 E)1J SXD C)G4,2 01146 2) OCT +000005000000 | 00772 128A LDQ COEN+1,4 01057 SXD C)G5,1 01147 OCT +000003000000 | 00773 FMP DATA+1,2 01060 139A CAL 2)+5 01150 OCT +000000000000 00774 FAD YPRED 01061 SXD 6)+4,4 01151 OCT +000001000000 | 00775 STO YPRED 01062 TSX (STH),4 01152 OCT +000202000000 | 00776 128A1 TXI *+1,4,1 01063 PZE 8)2T 01153 OCT +000006000000 | 00777 128A2 TXL 127A,4 01064 LXD 6)+4,4 01154 OCT +000203000000 | 01000 129A CLA DATA+1,1 01065 140A LDQ INVAR 01155 B) BCD 1 Z 01001 FSB YPRED 01066 STR 01156 OCT -377777777777 | 01002 STO DEV 01067 SXD 6)+4,4 01157 BCD 1 W | 01003 130A CAL 2)+5 01070 TSX (FIL),4 01160 OCT -377777777777 | 01004 TSX (STH),4 01071 LXD 6)+4,4 01161 BCD 1 X | 01005 PZE 8)2G 01072 141A TRA 135A 01162 OCT -377777777777 01006 131A LDQ RUN 01073 A)101 CLA 6)+3 01163 BCD 1 Y | 01007 STR 01074 STO 1)+3 01164 OCT -377777777777 | 01010 LDQ DATA+1,1 01075 CLA NOVPL 01165 BCD 1 + | 01011 STR 01076 SUB 6)+3 01166 OCT -377777777777 | 01012 LDQ YPRED 01077 ADD 1)+3 01167 BCD 1 - 01013 STR 01100 STO 1)+3 01170 OCT -377777777777 | 01014 LDQ DEV 01101 LDQ NOVPL 01171 BCD 1 * | 01015 STR 01102 MPY 2)+6 01172 OCT -377777777777 | 01016 TSX (FIL),4 01103 ALS 17 01173 BCD 1 / | 01017 LXD C)206,4 01104 SUB 2)+6 01174 OCT -377777777777 01020 131A1 TXI *+1,4,1 01105 ADD 1)+3 01175 BCD 1LOGE | 01021 SXD C)206,4 01106 STO C)101 01176 OCT -377777777777 | 01022 131A2 TXL D)419,4 01107 TRA 1,4 01177 BCD 1E ** | 01023 TRA 132A 01110 A)103 CLA 6)+3 01200 OCT -377777777777 | 01024 D)61F LXD C)G5,1 01111 STO 1)+3 01201 BCD 1ARCTAN 01025 D)41F LXD C)G4,2 01112 CLA NOVAR 01202 OCT -377777777777 | 01026 132A CLA IFTVAR 01113 SUB 6)+3 01203 BCD 1** | 01027 132A1 TZE 133A 01114 ADD 1)+3 01204 OCT -377777777777 | 01030 TPL D)81H 01115 STO 1)+3 01205 3) OCT +000000000000 | 01031 TRA D)81I 01116 LDQ NOVPL 01206 6) OCT +233000000000 133A BSS 01117 MPY 2)+6 01207 OCT +000000377777 | 01032 TSX PRINTX,4 01120 ALS 17 01210 OCT +000000000000 | 134A BSS 01121 SUB 2)+6 01211 OCT +000001000000 | 01033 TSX OUTPUT,4 01122 ADD 1)+3 01212 OCT +000000000000 | 01034 D)81H LXD C)103,4 01123 STO C)103 01213 OCT +000000000000 01035 135A CAL 2)+1 01124 TRA 1,4 01214 BCD 1) | 01036 SXD 6)+4,4 01125 A)105 CLA 6)+3 01215 8)3T BCD 1 (12A6 | 01037 TSX (RWT),4 01126 STO 1)+3 01216 BCD 16.6)) | 136A BSS 01127 CLA NOVAR 01217 BCD 1 7, F1 | 01040 TSX EXIT,4 01130 SUB 6)+3 01220 BCD 1(1H I 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 11 | | 01221 BCD 1E// 01311 BCD 1 I5 // 01401 BCD 1 | 01222 BCD 1 VALU 01312 BCD 1DATA = 01402 BCD 15.6, | 01223 BCD 1NO. 01313 BCD 1NO OF 01403 BCD 1 =F1 01224 BCD 1 VAR. 01314 BCD 1//13H 01404 BCD 1VS Y | 01225 BCD 1//20H 01315 BCD 1O I10 01405 BCD 1,12H) | 01226 BCD 1MENTS 01316 BCD 1BLEM N 01406 BCD 1 X(I2 | 01227 BCD 1L ELE 01317 BCD 12H PRO 01407 BCD 1 6H | 01230 BCD 1IAGONA 01320 BCD 1ON //1 01410 BCD 1 01231 BCD 1 D 01321 BCD 1GRESSI 01411 BCD 1 | 01232 8)1HI BCD 1 (24H0 01322 BCD 1SE RE 01412 BCD 1 | 01233 BCD 1.6) 01323 BCD 1STEPWI 01413 BCD 1 | 01234 8)V8 BCD 1 (6F12 01324 8)2Q BCD 1 (22H1 01414 BCD 16, | 01235 BCD 1NATED) 01325 BCD 1 //) 01415 BCD 1 =F15. 01236 BCD 1 TERMI 01326 BCD 1IATION 01416 BCD 1 Y | 01237 BCD 1ROBLEM 01327 BCD 1 DEV 01417 BCD 12H) VS | 01240 BCD 1ARD, P 01330 BCD 1CTED 01420 BCD 1X(I2,1 | 01241 BCD 1TROL C 01331 BCD 1 PREDI 01421 BCD 1H | 01242 BCD 1IN CON 01332 BCD 1 01422 8)18 BCD 1 (1H 7 01243 BCD 1ERROR 01333 BCD 1ACTUAL 01423 BCD 1 ) | 01244 8)S9 BCD 1 (42H 01334 BCD 1. 01424 BCD 1 F15.6 | 01245 BCD 15 ) 01335 BCD 1RUN NO 01425 BCD 1,3H) = | 01246 BCD 1 IS I 01336 BCD 1 01426 BCD 1S X(I2 | 01247 BCD 10. IT 01337 BCD 19H 01427 BCD 1,7H) V 01250 BCD 1HAN 15 01340 BCD 1TS /5 01430 BCD 1 X(I2 | 01251 BCD 1ATER T 01341 BCD 1 RESUL 01431 BCD 1 6H | 01252 BCD 1BE GRE 01342 BCD 1ACTUAL 01432 BCD 1 | 01253 BCD 1D NOT 01343 BCD 1ED VS 01433 BCD 1 | 01254 BCD 1 SHOUL 01344 BCD 1REDICT 01434 BCD 1 01255 BCD 11INVAR 01345 BCD 1 P 01435 BCD 1 | 01256 8)2T BCD 1 (45H 01346 BCD 1 01436 BCD 1, | 01257 BCD 1// ) 01347 BCD 13H 01437 BCD 1 F15.6 | 01260 BCD 19.3 / 01350 8)2L BCD 1 (1H04 01440 BCD 1,3H) = | 01261 BCD 1LE = F 01351 BCD 1F12.5) 01441 BCD 1S X(I2 01262 BCD 1VARIAB 01352 BCD 15,2H 01442 BCD 1,7H) V | 01263 BCD 1EMOVE 01353 BCD 1 F12. 01443 BCD 1 X(I2 | 01264 BCD 1L TO R 01354 BCD 1.5,3H 01444 BCD 1 6H | 01265 BCD 1F LEVE 01355 BCD 1H F12 01445 BCD 1 | 01266 BCD 1//29H 01356 BCD 1F6.0,2 01446 BCD 1 01267 BCD 1F10.3 01357 BCD 1 01447 BCD 1 | 01270 BCD 1BLE = 01360 8)2G BCD 1 (8H 01450 BCD 1 | 01271 BCD 1 VARIA 01361 BCD 115.6) 01451 BCD 115.6, | 01272 BCD 1 ENTER 01362 BCD 1Y =F 01452 BCD 1H) = F | 01273 BCD 1VEL TO 01363 BCD 1 VS 01453 BCD 1X(I2,3 01274 BCD 1H F LE 01364 BCD 1 Y 01454 BCD 1H) VS | 01275 BCD 12 //28 01365 BCD 11H 01455 BCD 1X(I2,7 | 01276 BCD 1= F12. 01366 8)1D BCD 1 (1H 2 01456 BCD 1H | 01277 BCD 1EEDOM 01367 BCD 15.6 ) 01457 8)13 BCD 1 (1H 7 | 01300 BCD 1 OF FR 01370 BCD 1 =F1 01460 BCD 1// ) 01301 BCD 1EGREES 01371 BCD 1VS Y 01461 BCD 1ODUCTS | 01302 BCD 1HTED D 01372 BCD 1,12H) 01462 BCD 1OSS PR | 01303 BCD 1H WEIG 01373 BCD 1 X(I2 01463 BCD 1ND CR | 01304 BCD 10 //30 01374 BCD 1 6H 01464 BCD 1RES A | 01305 BCD 1S = I1 01375 BCD 1 01465 BCD 1 SQUA 01306 BCD 1RIABLE 01376 BCD 1 01466 BCD 1UM OF | 01307 BCD 1 OF VA 01377 BCD 1 01467 BCD 1RAW S | 01310 BCD 118H NO 01400 BCD 1 01470 BCD 1 | 1 STEPWISE REGRESSION WITH VARIABLE TRANSFORMATIONS 12/01/64 PAGE 12 | 01471 BCD 1 | 01472 BCD 1 | 01473 BCD 1 | 01474 BCD 10H | 01475 8)U BCD 1(1H0 7 01476 BCD 12.4) | 01477 BCD 1 =F1 | 01500 BCD 1UM Y | 01501 BCD 1 S | 01502 8)P BCD 1 (17H 01503 BCD 1 ) | 01504 BCD 1=F12.4 | 01505 BCD 12,3H) | 01506 BCD 1UM X(I | 01507 BCD 1,8H S 01510 BCD 1=F12.4 | 01511 BCD 12,3H) | 01512 BCD 1M X( I | 01513 BCD 18H SU | 01514 BCD 1F12.4, 01515 BCD 1,3H) = | 01516 BCD 1 X( I2 | 01517 BCD 1H SUM | 01520 BCD 112.4,8 | 01521 BCD 1H) = F 01522 BCD 1( I2,3 | 01523 BCD 1 SUM X | 01524 BCD 11H | 01525 8)K BCD 1 (1H 1 | 01526 BCD 1LES//) 01527 BCD 1VARIAB | 01530 BCD 1 OF | 01531 BCD 1 SUM | 01532 BCD 1 | 01533 BCD 1 01534 BCD 1 | 01535 BCD 1 | 01536 BCD 19H | 01537 8)F BCD 1 (1H 4 | 01540 BCD 12.5)) 01541 8)A BCD 1 (6(F1 | 01542 BCD 1, 3I2) | 01543 BCD 1I2, I4 | 01544 BCD 1,1H 8 | 01545 BCD 1.5,3I5 01546 8)5 BCD 1 (3F10 | 1 1 12/01/64 PAGE 1 SUBROUTINE SETUP MPR30172 EQUIVALENCE (DATA,X) MPR30173 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30174 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30175 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30176 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30177 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30178 5, NAMEOP, NAMVAR MPR30179 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30180 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30181 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30182 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30183 4,NAMEOP(8),NAMVAR(4) MPR30184 DUMMY = 1.0 MPR30185 NN1 = 1 MPR30186 100 READ INPUT TAPE 5, 900, L1(NN1),L2(NN1),C1(NN1),L3(NN1),L4(NN1), MPR30187 1L5(NN1),L6(NN1),C2(NN1), L7(NN1),L8(NN1) MPR30188 IF (L1(NN1) -5) 110,800,850 MPR30189 110 IF (L4(NN1)) 120,120,130 MPR30190 120 L4(NN1)= 1 MPR30191 130 IF (L3(NN1)) 860,860,140 MPR30192 140 IF (L3(NN1)-10) 150,160,160 MPR30193 150 L9(NN1) = L4(NN1) MPR30194 GO TO 180 MPR30195 160 IF (L3(NN1)-20) 170,170,860 MPR30196 170 L9(NN1) =(L4(NN1) * (L4(NN1) +3))/2 MPR30197 C L9 = NO. OF RESULTS MPR30198 180 IF (L4(NN1)-10)184,184,890 MPR30199 C MAX OF 10 OPERATIONS PER COMMAND MPR30200 184 IF(L2(NN1)) 870,185,186 MPR30201 185 L2(NN1) = 1 MPR30202 C CHECK ON LEGALITY OF L1,L5,L7,L2,L6,L8 MPR30203 186 IF (L8(NN1)) 870, 187, 188 MPR30204 187 L8(NN1) = 1 MPR30205 188 IF (L7(NN1)-2) 850,190,190 MPR30206 190 IF (L7(NN1)-5) 195,850,850 MPR30207 195 L = L7(NN1)-1 MPR30208 LN= L8(NN1)+L9(NN1) - 1 MPR30209 L8N = L8(NN1) MPR30210 GO TO (200,210,220),L MPR30211 200 IF(LN - 70) 205, 205, 870 MPR30212 205 DO 206 I = L8N,LN MPR30213 206 W(I) = DUMMY MPR30214 GO TO 230 MPR30215 210 IF (LN-150) 215,215,870 MPR30216 215 DO 216 I = L8N,LN MPR30217 216 X(I) = DUMMY MPR30218 GO TO 230 MPR30219 220 IF (LN-20) 225,225,870 MPR30220 225 DO 226 I = L8N,LN MPR30221 226 Y(I) = DUMMY MPR30222 230 J = 1 MPR30223 IF(L1(NN1)) 850, 258, 233 MPR30224 233 I1 = L1(NN1)-1 MPR30225 I2 = L2(NN1) MPR30226 1 12/01/64 PAGE 2 234 IF (I1) 258,258,235 MPR30227 235 IF (I1-4) 240,850,850 MPR30228 240 I3 = I2 +L4(NN1) -1 MPR30229 GO TO (245,250,255),I1 MPR30230 245 IF (I3-70) 246,246,870 MPR30231 246 DO 247 I = I2,I3 MPR30232 IF (W(I)) 880,880,247 MPR30233 247 CONTINUE MPR30234 GO TO 258 MPR30235 250 IF(I3-130) 251,251,870 MPR30236 251 DO 252 I = I2,I3 MPR30237 IF (X(I)) 880,880,252 MPR30238 252 CONTINUE MPR30239 GO TO 258 MPR30240 255 IF (I3-20) 256,256,870 MPR30241 256 DO 257 I = I2,I3 MPR30242 IF (Y(I)) 880,880,257 MPR30243 257 CONTINUE MPR30244 258 GO TO (260,270),J MPR30245 C MPR30246 260 J = 2 MPR30247 IF(L5(NN1)) 850,258,265 MPR30248 265 I1 = L5(NN1) -1 MPR30249 I2 = L6(NN1) MPR30250 GO TO 234 MPR30251 C MPR30252 270 NN1 = NN1 + 1 MPR30253 IF ( NN1 - 26) 100, 100, 1000 MPR30254 C ERROR RETURNS MPR30255 850 KK = 1 MPR30256 GO TO 1000 MPR30257 860 KK = 2 MPR30258 GO TO 1000 MPR30259 870 KK = 3 MPR30260 GO TO 1000 MPR30261 880 KK = 4 MPR30262 GO TO 1000 MPR30263 890 KK = 5 MPR30264 1000 WRITE OUTPUT TAPE 6, 910, RUN,NN1 ,KK MPR30265 PRINT 910, RUN,NN1,KK MPR30266 CALL EXIT MPR30267 C MPR30268 C NORMAL RETURN MPR30269 800 NN1 = NN1 - 1 MPR30270 RETURN MPR30271 900 FORMAT (2I6,F12.6,4I6,F12.6,2I6) MPR30272 910 FORMAT (107H ERROR IN -- (1) CLASS, (2) OP CODE, (3) INDEX, (4) A MPR30273 1OR B NOT PREVIOUSLY DEFINED, (5) INDEX OF OPERATIONS. /14H SET OF MPR30274 2DATA = F12.6, 12H, CARD NO. = I3, 10H, ERROR = I3, 1H./1H1) MPR30275 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 12/01/64 PAGE 3 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 518 01006 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14089 33411 NAMVAR 14081 33401 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 517 01005 B 507 00773 C 497 00761 CTEMP 432 00660 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | DUMMY 422 00646 I1 421 00645 I2 420 00644 I3 419 00643 I 418 00642 | J 417 00641 KK 416 00640 L8N 415 00637 LN 414 00636 L 413 00635 | 0 SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT STATEMENTS 0 EFN LOC EFN LOC EFN LOC EFN LOC EFN LOC | 8)S4 900 00625 8)SE 910 00620 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 1) 406 00626 2) 351 00537 3) 363 00553 6) 364 00554 C)G0 408 00630 C)G1 409 00631 C)G2 410 00632 C)G3 411 00633 C)G4 412 00634 D)20T 165 00245 | D)21N 265 00411 D)21R 295 00447 D)80L 140 00214 D)81N 262 00406 D)81T 304 00460 | D)820 314 00472 D)A0L 139 00213 D)A1T 303 00457 E)1 22 00026 E)2 55 00067 | E)H 126 00176 E)P 153 00231 E)10 185 00271 E)15 208 00320 E)16 212 00324 | E)1C 230 00346 E)1D 234 00352 E)1J 252 00374 E)1K 256 00400 E)1O 268 00414 E)21 342 00526 E)DR 294 00446 E)HN 264 00410 E)LT 302 00456 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | EXIT 5 00005 (FIL) 3 00003 (RTN) 1 00001 (SPH) 4 00004 (STH) 2 00002 | (TSH) 0 00000 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | 1 12/01/64 PAGE 4 | | | EXIT (FIL) (RTN) (SPH) (STH) (TSH) 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 100 10 00027 110 13 00070 120 14 00073 130 15 00075 140 16 00101 | 150 17 00105 160 19 00110 170 20 00114 180 21 00127 184 22 00133 | 185 23 00137 186 24 00141 187 25 00145 188 26 00147 190 27 00154 195 28 00160 200 32 00177 205 33 00203 206 34 00206 210 36 00215 | 215 37 00221 216 38 00224 220 40 00232 225 41 00236 226 42 00241 | 230 43 00246 233 45 00256 234 47 00265 235 48 00272 240 49 00276 | 245 51 00305 246 52 00311 247 54 00330 250 56 00333 251 57 00337 | 252 59 00356 255 61 00361 256 62 00365 257 64 00404 258 65 00412 260 66 00415 265 68 00425 270 71 00435 850 73 00450 860 75 00453 | 870 77 00461 880 79 00464 890 81 00467 1000 82 00473 800 87 00527 | 1 1 SUBROUTINE TO SET UP PERMANENT BCD CODES IN PROPER CELLS 12/01/64 PAGE 1 SUBROUTINE PERM(A,B,C,D,E,F,G,H,I,J,K,L,AA,BB,CC,DD,EE,FF,GG,HH, MPR30278 1 II,JJ,KK,LL) MPR30279 C MPR30280 AA = A MPR30281 BB = B MPR30282 CC = C MPR30283 DD = D MPR30284 EE = E MPR30285 FF = F MPR30286 GG = G MPR30287 HH = H MPR30288 II = I MPR30289 JJ = J MPR30290 KK = K MPR30291 LL = L MPR30292 RETURN MPR30293 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 SUBROUTINE TO SET UP PERMANENT BCD CODES IN PROPER CELLS 12/01/64 PAGE 2 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 91 00133 32561 77461 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 6) 85 00125 | 1 1 SUBROUTINE SWITCH 12/01/64 PAGE 1 SUBROUTINE SWITCH MPR30295 EQUIVALENCE (DATA,X) MPR30296 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30297 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30298 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30299 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30300 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30301 5, NAMVAR, NAMEOP MPR30302 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30303 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30304 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30305 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30306 4,NAMEOP(8),NAMVAR(4) MPR30307 C NN1 = NO. OF TRANSFORMATION COMMANDS MPR30308 DO 1000 N = 1, NN1 MPR30309 LL1 = L1(N) MPR30310 LL2 = L2(N) MPR30311 LL3 = L3(N) MPR30312 LL4 = L4(N) MPR30313 LL5 = L5(N) MPR30314 LL6 = L6(N) MPR30315 LL7 = L7(N)-1 MPR30316 LL8 = L8(N) MPR30317 LL9 = L9(N) MPR30318 C L9 = NO. OF RESULTS MPR30319 C L4 = NO. OF OPERATIONS MPR30320 M1 = LL3 MPR30321 M2 = 21 MPR30322 DO 110 I = 1,LL4 MPR30323 IF (LL1)5,5,10 MPR30324 5 A(I) = C1(N) MPR30325 GO TO 40 MPR30326 10 GO TO (15,20,25,30),LL1 MPR30327 15 A(I)=Z (LL2) + C1(N) MPR30328 GO TO 40 MPR30329 20 A(I) = W(LL2) + C1(N) MPR30330 GO TO 40 MPR30331 25 A(I) = X(LL2) + C1(N) MPR30332 GO TO 40 MPR30333 30 LL2 = NOVAR-1+L2(N) MPR30334 A(I) = X(LL2) + C1(N) MPR30335 GO TO 40 MPR30336 40 IF (LL5)45,45,50 MPR30337 45 B(I) = C2(N) MPR30338 GO TO 100 MPR30339 50 GO TO (55,60,65,70),LL5 MPR30340 55 B(I) = Z (LL6) + C2(N) MPR30341 GO TO 100 MPR30342 60 B(I) = W(LL6) + C2(N) MPR30343 GO TO 100 MPR30344 65 B(I) = X(LL6) + C2(N) MPR30345 GO TO 100 MPR30346 70 LL6 = NOVAR + L6(N) - 1 MPR30347 B(I) = X(LL6)+C2(N) MPR30348 100 LL2 = LL2 + 1 MPR30349 LL6 = LL6 + 1 MPR30350 1 SUBROUTINE SWITCH 12/01/64 PAGE 2 110 CONTINUE MPR30351 120 GO TO (210,220,230,240,250,260,270,280,290,300,310,310, MPR30352 1 310,310,310,310,310,310,390,400,500),M1 MPR30353 210 DO 215 K = 1,LL4 MPR30354 215 C(K) = A(K) + B(K) MPR30355 GO TO 450 MPR30356 C MPR30357 220 DO 225 K= 1,LL4 MPR30358 225 C(K) = A(K)-B(K) MPR30359 GO TO 450 MPR30360 C MPR30361 230 DO 235 K= 1,LL4 MPR30362 235 C(K) = A(K) * B(K) MPR30363 GO TO 450 MPR30364 C MPR30365 240 DO 245 K=1,LL4 MPR30366 245 C(K) = A(K) / B(K) MPR30367 GO TO 450 MPR30368 C MPR30369 250 DO 255 K=1,LL4 MPR30370 255 C(K) = LOGF(A(K)) MPR30371 GO TO 450 MPR30372 C MPR30373 260 DO 265 K = 1,LL4 MPR30374 265 C(K) = EXPF(A(K)) MPR30375 GO TO 450 MPR30376 C MPR30377 270 DO 275 K = 1,LL4 MPR30378 275 C(K) = ATANF(A(K)) MPR30379 GO TO 450 MPR30380 C MPR30381 280 DO 285 K = 1,LL4 MPR30382 285 C(K) = A(K) ** B(K) MPR30383 GO TO 450 MPR30384 C MPR30385 290 DO 295 K=1,LL4 MPR30386 CALL FCN1(A,B,C) MPR30387 295 CONTINUE MPR30388 GO TO 450 MPR30389 C MPR30390 300 IF(M2-21) 303,301,303 MPR30391 301 DO 302 K = 1,LL4 MPR30392 302 CTEMP(K) = A(K) MPR30393 GO TO 305 MPR30394 303 DO 304 K= 1,LL4 MPR30395 304 CTEMP(K) = C(K) MPR30396 305 KK = 1 MPR30397 DO 307 K1 = 1,LL4 MPR30398 C(KK) = CTEMP(K1) MPR30399 K1 = K1 MPR30400 KK=KK+1 MPR30401 DO 306 K=K1,LL4 MPR30402 C(KK) = CTEMP(K1) * CTEMP(K) MPR30403 306 KK = KK+1 MPR30404 307 CONTINUE MPR30405 M2 = 21 MPR30406 1 SUBROUTINE SWITCH 12/01/64 PAGE 3 GO TO 450 MPR30407 C MPR30408 310 M1 = M1-10 MPR30409 M2 = 10 MPR30410 GO TO 120 MPR30411 C MPR30412 390 DO 395 K= 1,LL4 MPR30413 CALL FCN2(A,B,C) MPR30414 395 CONTINUE MPR30415 GO TO 450 MPR30416 C MPR30417 400 DO 405 K=1,LL4 MPR30418 CALL FCN3 (A,B,C) MPR30419 405 CONTINUE MPR30420 GO TO 450 MPR30421 C MPR30422 450 M1 = M2 MPR30423 GO TO 120 MPR30424 C MPR30425 500 GO TO ( 510, 520, 530), LL7 MPR30426 510 DO 515 M = 1,LL9 MPR30427 W(LL8) = C(M) MPR30428 515 LL8 = LL8+1 MPR30429 GO TO 1000 MPR30430 520 DO 525 M = 1,LL9 MPR30431 X(LL8) = C(M) MPR30432 525 LL8= LL8+1 MPR30433 GO TO 1000 MPR30434 530 LL8 = LL8 + NOVAR-1 MPR30435 DO 535 M=1,LL9 MPR30436 X(LL8) = C(M) MPR30437 535 LL8 = LL8 + 1 MPR30438 GO TO 1000 MPR30439 C MPR30440 1000 CONTINUE MPR30441 C MPR30442 RETURN MPR30443 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 SUBROUTINE SWITCH 12/01/64 PAGE 4 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 565 01065 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 564 01064 B 554 01052 C 544 01040 CTEMP 479 00737 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | K1 469 00725 KK 468 00724 LL1 467 00723 LL2 466 00722 LL3 465 00721 | LL4 464 00720 LL5 463 00717 LL6 462 00716 LL7 461 00715 LL8 460 00714 | LL9 459 00713 M1 458 00712 M2 457 00711 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 2) 437 00665 6) 440 00670 C)G0 446 00676 C)G1 447 00677 C)G2 448 00700 | C)G3 449 00701 C)G4 450 00702 C)G5 451 00703 C)G6 452 00704 C)G8 453 00705 | C)G9 454 00706 C)GA 455 00707 C)GB 456 00710 D)204 69 00105 D)205 75 00113 D)206 81 00121 D)207 87 00127 D)209 102 00146 D)20B 109 00155 D)20C 115 00163 | D)20D 121 00171 D)20E 127 00177 D)402 62 00076 D)40I 153 00231 D)41M 309 00465 | D)424 384 00600 D)427 398 00616 D)42C 428 00654 D)602 61 00075 D)61M 308 00464 | D)80V 215 00327 D)822 377 00571 E)5 74 00112 E)6 80 00120 E)7 86 00126 | E)C 114 00162 E)D 120 00170 E)E 126 00176 E)M 184 00270 E)P 194 00302 E)S 204 00314 E)V 214 00326 E)12 227 00343 E)15 239 00357 E)18 251 00373 | E)1B 264 00410 E)1E 277 00425 E)1P 335 00517 E)1Q 339 00523 E)1R 348 00534 | E)1U 359 00547 E)22 376 00570 E)25 393 00611 E)28 407 00627 E)2B 426 00652 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT ATAN 2 00002 EXP 1 00001 EXP(3 3 00003 FCN1 4 00004 FCN2 5 00005 | FCN3 6 00006 LOG 0 00000 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | 1 SUBROUTINE SWITCH 12/01/64 PAGE 5 | | ATAN EXP EXP(3 FCN1 FCN2 FCN3 LOG | 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 5 20 00102 10 22 00106 15 23 00114 20 25 00122 25 27 00130 30 29 00134 40 32 00147 45 33 00152 50 35 00156 55 36 00164 | 60 38 00172 65 40 00200 70 42 00204 100 44 00215 110 46 00227 | 120 47 00232 210 48 00257 215 49 00262 220 51 00271 225 52 00274 | 230 54 00303 235 55 00306 240 57 00315 245 58 00320 250 60 00330 | 255 61 00333 260 63 00344 265 64 00347 270 66 00360 275 67 00363 280 69 00374 285 70 00377 290 72 00411 295 74 00422 300 76 00426 | 301 77 00433 302 78 00436 303 80 00443 304 81 00446 305 82 00452 | 306 89 00505 307 90 00513 310 93 00524 390 96 00535 395 98 00544 | 400 100 00550 405 102 00561 450 104 00564 500 106 00572 510 107 00575 | 515 109 00603 520 111 00613 525 113 00621 530 115 00631 535 118 00644 1000 120 00655 | 1 1 SUBROUTINE PRTTRF 12/01/64 PAGE 1 SUBROUTINE PRTTRF MPR30445 EQUIVALENCE (DATA,X) MPR30446 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30447 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30448 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30449 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30450 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30451 5, NAMVAR, NAMEOP MPR30452 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30453 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30454 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30455 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30456 4,NAMEOP(8),NAMVAR(4) MPR30457 WRITE OUTPUT TAPE 6, 10 MPR30458 DO 200 I = 1,NN1 MPR30459 LL1 = L1(I) MPR30460 LL2 = L2(I) MPR30461 LL3 = L3(I) MPR30462 LL4 = L4(I) MPR30463 LL5 = L5(I) MPR30464 LL6 = L6(I) MPR30465 LL7 = L7(I) MPR30466 LL8 = L8(I) MPR30467 NOPER = NAMEOP(LL3) MPR30468 NVAR1 = NAMVAR(LL1) MPR30469 NVAR5 = NAMVAR(LL5) MPR30470 NVAR7 = NAMVAR(LL7) MPR30471 GO TO (20, 20, 20, 20, 40, 40, 40, 20, 60, 80, 100, 100, 100, 100,MPR30472 1 120, 120, 120, 100, 60, 60), LL3 MPR30473 C MPR30474 C MPR30475 20 DO 25 J = 1, LL4 MPR30476 WRITE OUTPUT TAPE 6, 11, NVAR1, LL2, C1(I), NOPER, NVAR5, MPR30477 1LL6, MPR30478 1 C2(I), NVAR7, LL8 MPR30479 LL2 = LL2 + 1 MPR30480 LL6 = LL6 + 1 MPR30481 LL8 = LL8 + 1 MPR30482 25 CONTINUE MPR30483 GO TO 200 MPR30484 C MPR30485 C MPR30486 40 DO 45 J = 1, LL4 MPR30487 WRITE OUTPUT TAPE 6, 12, NOPER, NVAR1, LL2, C1(I), NVAR7, LL8 MPR30488 LL2 = LL2 + 1 MPR30489 LL8 = LL8 + 1 MPR30490 45 CONTINUE MPR30491 GO TO 200 MPR30492 C MPR30493 C MPR30494 60 IF (LL3 - 19) 62, 64, 66 MPR30495 62 NOFCN = 1 MPR30496 GO TO 70 MPR30497 64 NOFCN = 2 MPR30498 GO TO 70 MPR30499 66 NOFCN = 3 MPR30500 1 SUBROUTINE PRTTRF 12/01/64 PAGE 2 70 DO 75 J = 1,LL4 MPR30501 WRITE OUTPUT TAPE 6, 16, NOFCN, NVAR7, LL8 MPR30502 LL8 = LL8 + 1 MPR30503 75 CONTINUE MPR30504 GO TO 200 MPR30505 C MPR30506 C MPR30507 80 NOPER = NAMEOP(3) MPR30508 DO 95 J= 1, LL4 MPR30509 WRITE OUTPUT TAPE 6, 13, NVAR1, LL2, C1(I), NVAR7, LL8 MPR30510 J = J MPR30511 LL8 = LL8 + 1 MPR30512 J2 = LL2 MPR30513 DO 85 J1 = J,LL4 MPR30514 WRITE OUTPUT TAPE 6, 11, NVAR1, LL2, C1(I), NOPER, NVAR1, J2, MPR30515 1C1(I), NVAR7, LL8 MPR30516 J2 = J2 +1 MPR30517 LL8 = LL8 + 1 MPR30518 85 CONTINUE MPR30519 LL2 = LL2 + 1 MPR30520 95 CONTINUE MPR30521 GO TO 200 MPR30522 C MPR30523 C MPR30524 100 LL3 = LL3 - 10 MPR30525 NOPER = NAMEOP(LL3) MPR30526 DO 110 J = 1,LL4 MPR30527 WRITE OUTPUT TAPE 6, 11, NVAR1, LL2, C1(I), NOPER, NVAR5, LL6, MPR30528 2C2(I), NVAR7, LL8 MPR30529 J = J MPR30530 LL8= LL8 + 1 MPR30531 J2 = LL2 MPR30532 J6 = LL6 MPR30533 DO 105 J1 = J, LL4 MPR30534 WRITE OUTPUT TAPE 6, 14, NVAR1, LL2, C1(I), NOPER, NVAR5, LL6, MPR30535 3C2(I), NVAR1, J2, MPR30536 1 C1(I), NOPER, NVAR5, J6, C2(I), NVAR7, LL8 MPR30537 J2 = J2 + 1 MPR30538 J6 = J6 + 1 MPR30539 LL8 = LL8 + 1 MPR30540 105 CONTINUE MPR30541 LL2 = LL2 + 1 MPR30542 LL6 = LL6 + 1 MPR30543 110 CONTINUE MPR30544 GO TO 200 MPR30545 C MPR30546 C MPR30547 120 LL3 = LL3 - 10 MPR30548 NOPER = NAMEOP(LL3) MPR30549 DO 130 J = 1, LL4 MPR30550 WRITE OUTPUT TAPE 6, 12, NOPER, NVAR1, LL2, C1(I), NVAR7, LL8 MPR30551 J = J MPR30552 LL8= LL8 + 1 MPR30553 J2 = LL2 MPR30554 DO 125 J1 = J, LL4 MPR30555 WRITE OUTPUT TAPE 6, 15, NOPER, NVAR1, LL2, C1(I), NOPER, NVAR1, MPR30556 1 SUBROUTINE PRTTRF 12/01/64 PAGE 3 4J2, C1(I), NVAR7 MPR30557 1, LL8 MPR30558 J2 = J2 + 1 MPR30559 LL8 = LL8 + 1 MPR30560 125 CONTINUE MPR30561 LL2 = LL2 + 1 MPR30562 130 CONTINUE MPR30563 GO TO 200 MPR30564 C MPR30565 C MPR30566 200 CONTINUE MPR30567 C MPR30568 RETURN MPR30569 C MPR30570 10 FORMAT (72H1 VARIABLMPR30571 1E TRANSFORMATIONS / / ) MPR30572 11 FORMAT (5H ( A2,1H( I3, 3H) + F12.6, 2H) A2, 3H ( A2, 1H( MPR30573 1 I3, 3H) + F12.6, 4H) = A2, 1H( I3, 1H) ) MPR30574 12 FORMAT (4H A6, 1H( A2, 1H( I3, 3H) + F12.6, 4H) = A2, 1H( MPR30575 1 I3, 1H) ) MPR30576 13 FORMAT (5H ( A2, 1H( I3, 3H) + F12.6, 4H) = A2, 1H( I3, MPR30577 1 1H) ) MPR30578 14 FORMAT (6H (( A2, 1H( I3, 3H) + F12.6, 1H) A2, 1H( A2, 1H( MPR30579 1 I3, 3H) + F12.6, 7H)) * (( A2, 1H( I3, 3H) + F12.6, 1H) A2, MPR30580 2 1H( A2, 1H( I3, 3H) + F12.6, 4H)) = A2, 1H( I3, 1H) ) MPR30581 15 FORMAT (4H A6, 1H( A2, 1H( I3, 3H) + F12.6, 4H) * A6, 1H( MPR30582 1 A2, 1H( I3, 3H) + F12.6, 4H) = A2, 1H( I3, 1H) ) MPR30583 16 FORMAT (7H FCN I1, 13H CALCULATES A2, 1H( I3, 1H) ) MPR30584 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 SUBROUTINE PRTTRF 12/01/64 PAGE 4 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 707 01303 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 706 01302 B 696 01270 C 686 01256 CTEMP 621 01155 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | J2 611 01143 J6 610 01142 J 609 01141 LL1 608 01140 LL2 607 01137 | LL3 606 01136 LL4 605 01135 LL5 604 01134 LL6 603 01133 LL7 602 01132 | LL8 601 01131 NOFCN 600 01130 NOPER 599 01127 NVAR1 598 01126 NVAR5 597 01125 NVAR7 596 01124 | 0 SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT STATEMENTS | 0 EFN LOC EFN LOC EFN LOC EFN LOC EFN LOC | 8)A 10 01114 8)B 11 01076 8)C 12 01054 8)D 13 01040 8)E 14 01024 | 8)F 15 00765 8)G 16 00741 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 2) 461 00715 6) 467 00723 C)G0 589 01115 C)G1 590 01116 C)G2 591 01117 | C)G3 592 01120 C)G4 593 01121 C)203 594 01122 C)205 595 01123 D)20U 453 00705 | D)403 76 00114 D)406 119 00167 D)40G 195 00303 D)40L 275 00423 D)40Q 383 00577 D)60G 194 00302 D)60L 274 00422 E)2 72 00110 E)5 115 00163 E)8 149 00225 | E)F 185 00271 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | (FIL) 1 00001 (STH) 0 00000 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | 1 SUBROUTINE PRTTRF 12/01/64 PAGE 5 | | | (FIL) (STH) 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 20 28 00111 25 34 00160 40 36 00164 45 41 00222 60 43 00226 | 62 44 00232 64 46 00235 66 48 00240 70 49 00242 75 53 00266 | 80 55 00272 85 67 00374 95 69 00401 100 71 00405 105 86 00547 110 89 00557 120 91 00563 125 104 00674 130 106 00701 200 108 00706 | 1 1 SUBROUTINE SUM 12/01/64 PAGE 1 SUBROUTINE SUM MPR30586 C CALCUALTION OF SUMS WITH WEIGHTS = 1.0 MPR30587 C MPR30588 EQUIVALENCE (DATA,X) MPR30589 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30590 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30591 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30592 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30593 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30594 5, NAMVAR, NAMEOP MPR30595 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30596 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30597 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30598 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30599 4,NAMEOP(8),NAMVAR(4) MPR30600 C MPR30601 161 WRITE TAPE 3, (X(L), L = 1, NOVAR), RUN MPR30602 180 DO 190 I = 1, NOVAR MPR30603 200 VECTOR(I, NOVAR + 1) = VECTOR(I, NOVAR + 1) + DATA(I) MPR30604 210 DO 220 J = I, NOVAR MPR30605 220 VECTOR (I,J) = VECTOR(I,J) + DATA(I) * DATA(J) MPR30606 190 CONTINUE MPR30607 170 VECTOR(NVP1, NVP1) = VECTOR(NVP1, NVP1) + 1.0 MPR30608 RETURN MPR30609 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 SUBROUTINE SUM 12/01/64 PAGE 2 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 209 00321 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 208 00320 B 198 00306 C 188 00274 CTEMP 123 00173 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | I 113 00161 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT 1) 104 00150 2) 94 00136 3) 97 00141 6) 98 00142 A)102 81 00121 | C)G1 108 00154 C)G2 109 00155 C)100 110 00156 C)101 111 00157 C)102 112 00160 | D)404 58 00072 D)604 57 00071 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT (STB) 0 00000 (WLR) 1 00001 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | | (STB) (WLR) | 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 161 6 00015 180 11 00031 200 12 00050 210 13 00053 220 14 00073 | 190 15 00102 170 16 00112 | 1 1 CALCULATION SUMS WHEN VARIABLE WEIGHTS 12/01/64 PAGE 1 SUBROUTINE WGTSUM MPR30612 C MPR30613 EQUIVALENCE (DATA,X) MPR30614 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30615 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30616 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30617 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30618 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30619 5, NAMVAR, NAMEOP MPR30620 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30621 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30622 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30623 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30624 4,NAMEOP(8),NAMVAR(4) MPR30625 C MPR30626 521 WRITE TAPE 3, (X(L), L = 1, NOVAR), RUN MPR30627 530 DO 540 I = 1, NOVAR MPR30628 550 VECTOR (I,NOVAR + 1) = VECTOR (I, NOVAR + 1) + DATA (I) * WHT MPR30629 560 DO 540 J = I, NOVAR MPR30630 540 VECTOR (I, J) = VECTOR (I, J) + DATA (I) * DATA (J) * WHT MPR30631 510 VECTOR (NVP1, NVP1) = VECTOR (NVP1, NVP1) +WHTMPR30632 RETURN MPR30633 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 CALCULATION SUMS WHEN VARIABLE WEIGHTS 12/01/64 PAGE 2 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 211 00323 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 210 00322 B 200 00310 C 190 00276 CTEMP 125 00175 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | I 115 00163 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT 1) 106 00152 2) 97 00141 6) 100 00144 A)102 84 00124 C)G1 110 00156 | C)G2 111 00157 C)100 112 00160 C)101 113 00161 C)102 114 00162 D)804 59 00073 | D)A04 58 00072 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT (STB) 0 00000 (WLR) 1 00001 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | | (STB) (WLR) | 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 521 6 00015 530 11 00031 550 12 00050 560 13 00054 540 14 00074 | 510 15 00115 | 1 1 CALCULATES THE NORMAL EQUATIONS 12/01/64 PAGE 1 C CALCULATION OF RESIDUAL SUMS OF SQUARES AND CROSS PRODUCTS MPR30636 SUBROUTINE RESID MPR30637 C MPR30638 EQUIVALENCE (DATA,X) MPR30639 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30640 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30641 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30642 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30643 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30644 5, NAMVAR, NAMEOP MPR30645 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30646 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30647 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30648 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30649 4,NAMEOP(8),NAMVAR(4) MPR30650 C MPR30651 650 IF(IFCNST) 900,651,735 MPR30652 651 IF(VECTOR(NOVPL,NOVPL)) 652,652,655 MPR30653 652 PRINT 654 MPR30654 GO TO 910 MPR30655 655 DO 660 I = 1, NOVAR MPR30656 670 DO 660 J = I, NOVAR MPR30657 660 VECTOR (I,J) = VECTOR (I,J) - (VECTOR(I,NOVPL) * VECTOR (J,NOVPL) MPR30658 - / VECTOR (NOVPL, NOVPL)) MPR30659 680 DO 690 I = 1, NOVAR MPR30660 690 AVE(I) = VECTOR(I,NOVPL) / VECTOR(NOVPL,NOVPL) MPR30661 700 IF (IFAVE) 900, 710, 735 MPR30662 710 WRITE OUTPUT TAPE 6, 50 MPR30663 720 WRITE OUTPUT TAPE 6, 20, (I,AVE(I), I=1,NOVMI) MPR30664 730 WRITE OUTPUT TAPE 6, 25, AVE(NOVAR) MPR30665 735 IF (IFRESD) 900, 740, 780 MPR30666 740 WRITE OUTPUT TAPE 6, 55 MPR30667 750 WRITE OUTPUT TAPE 6, 35, ((I,J,VECTOR(I,J),J=I,NOVMI),I=1,NOVMI) MPR30668 760 WRITE OUTPUT TAPE 6, 40, (I,VECTOR(I,NOVAR),I=1,NOVMI) MPR30669 770 WRITE OUTPUT TAPE 6, 45, VECTOR(NOVAR,NOVAR) MPR30670 780 RETURN MPR30671 C MPR30672 C ERROR RETURNS MPR30673 900 PRINT 905 MPR30674 910 REWIND 3 MPR30675 CALL EXIT MPR30676 C MPR30677 20 FORMAT (1H 11H SUM X( I2,3H) = F12.4,8H SUM X( I2,3H) =F12.4,MPR30678 18H SUM X( I2,3H) =F12.4,8H SUM X(I2,3H) =F12.4 ) MPR30679 25 FORMAT (17H SUM Y =F12.4) MPR30680 35 FORMAT (1H 7H X(I2,7H) VS X(I2,3H) = F15.6, MPR30681 1 6H X(I2,7H) VS X(I2,3H) = F15.6, MPR30682 2 6H X(I2,7H) VS X(I2,3H) = F15.6 ) MPR30683 40 FORMAT (1H 7H X(I2,12H) VS Y =F15.6, MPR30684 1 6H X(I2,12H) VS Y =F15.6, MPR30685 2 6H X(I2,12H) VS Y =F15.6 ) MPR30686 45 FORMAT (1H 21H Y VS Y =F15.6) MPR30687 50 FORMAT (1H063H AVERAGE VALUE OFMPR30688 - VARIABLES// ) MPR30689 55 FORMAT(1H077H RESIDUAL SUMS OF SQUAMPR30690 1 CALCULATES THE NORMAL EQUATIONS 12/01/64 PAGE 2 -RES AND CROSS PRODUCTS//) MPR30691 654 FORMAT (31H ZERO NUMBER OF DATA. SO LONG.) MPR30692 905 FORMAT (42H ERROR IN CONTROL CARD, PROBLEM TERMINATED) MPR30693 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 CALCULATES THE NORMAL EQUATIONS 12/01/64 PAGE 3 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 536 01030 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 535 01027 B 525 01015 C 515 01003 CTEMP 450 00702 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | I 440 00670 J 439 00667 | 0 SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT STATEMENTS | 0 EFN LOC EFN LOC EFN LOC EFN LOC EFN LOC 8)K 20 00652 8)P 25 00627 8)13 35 00622 8)18 40 00565 8)1D 45 00531 | 8)1I 50 00523 8)1N 55 00506 8)KE 654 00467 8)S9 905 00460 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 1) 427 00653 2) 287 00437 6) 291 00443 A)100 261 00405 A)105 274 00422 C)G0 431 00657 C)G1 432 00660 C)100 433 00661 C)101 434 00662 C)102 435 00663 | C)105 436 00664 C)201 437 00665 C)203 438 00666 D)405 72 00110 D)808 109 00155 | D)80D 148 00224 D)80L 233 00351 D)80M 239 00357 D)A0D 147 00223 D)A0M 238 00356 | D)C08 108 00154 D)C0L 232 00350 E)1 23 00027 E)E 153 00231 E)1M 241 00361 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | EXIT 4 00004 (FIL) 1 00001 (RWT) 3 00003 (SPH) 0 00000 (STH) 2 00002 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | | EXIT (FIL) (RWT) (SPH) (STH) 1 CALCULATES THE NORMAL EQUATIONS 12/01/64 PAGE 4 | | 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 650 15 00023 651 16 00030 652 17 00033 655 19 00041 670 20 00061 660 21 00111 680 22 00136 690 23 00156 700 24 00164 710 25 00170 | 720 26 00174 730 31 00214 735 33 00225 740 34 00232 750 35 00240 | 760 43 00307 770 48 00341 780 50 00352 900 52 00362 910 53 00367 | 1 1 PRINT TRANSFORMED VARIABLES (X) 12/01/64 PAGE 1 SUBROUTINE PRINTX MPR30696 EQUIVALENCE (DATA,X) MPR30697 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30698 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30699 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30700 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30701 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30702 5, NAMVAR, NAMEOP, NODATA MPR30703 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30704 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30705 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30706 3L9(25), C1(26),C2(26) MPR30707 4, NAMEOP(8), NAMVAR(4), IK(9) MPR30708 C MPR30709 WRITE OUTPUT TAPE 6, 300 MPR30710 DO 100 I = 1, NOVAR, 9 MPR30711 REWIND 3 MPR30712 IK = I MPR30713 DO 10 J = 1, 8 MPR30714 J = J MPR30715 10 IK(J+1) = IK + J MPR30716 IL1 = NOVAR - (IK +8) MPR30717 IF (IL1) 20, 30, 30 MPR30718 20 IL1 = 9 + IL1 MPR30719 IL2 = NOVAR MPR30720 GO TO 40 MPR30721 30 IL1 = 9 MPR30722 IL2 = IK +8 MPR30723 40 WRITE OUTPUT TAPE 6, 310, (IK(L), L = 1, IL1) MPR30724 DO 50 N = 1, NODATA MPR30725 READ TAPE 3, (X(M), M = 1, NOVAR), RUN MPR30726 WRITE OUTPUT TAPE 6, 320, RUN, (X(M), M= I, IL2) MPR30727 50 CONTINUE MPR30728 100 CONTINUE MPR30729 RETURN MPR30730 300 FORMAT(75H1 TRANSFMPR30731 1ORMED VARABLES (X) / / ) MPR30732 310 FORMAT( / 12H RUN NO. I10, 8I12 ) MPR30733 320 FORMAT ( F8.0, 4H 9F12.6) MPR30734 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 PRINT TRANSFORMED VARIABLES (X) 12/01/64 PAGE 2 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 153 00231 14076 33374 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NODATA 14077 33375 NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 | NOSTEP 31539 75463 NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 | VECTOR 31537 75461 WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 | Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT IK 152 00230 | 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | IL1 143 00217 IL2 142 00216 I 141 00215 J 140 00214 | 0 SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT STATEMENTS 0 EFN LOC EFN LOC EFN LOC EFN LOC EFN LOC | 8)9C 300 00213 8)9M 310 00175 8)A0 320 00167 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 2) 105 00151 6) 110 00156 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | (FIL) 1 00001 (RLR) 4 00004 (RWT) 2 00002 (STH) 0 00000 (TSB) 3 00003 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | (FIL) (RLR) (RWT) (STH) (TSB) | 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 10 15 00040 20 18 00056 30 21 00064 40 23 00071 50 40 00140 | 100 41 00142 1 1 SUBROUTINE CALC 12/01/64 PAGE 1 SUBROUTINE CALC MPR30736 EQUIVALENCE (DATA,X) MPR30737 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30738 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30739 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30740 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30741 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30742 5, NAMVAR, NAMEOP MPR30743 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30744 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30745 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30746 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30747 4,NAMEOP(8),NAMVAR(4) MPR30748 781 ASSIGN 1320 TO NUMBER MPR30749 1000 NOSTEP = NOSTEP + 1 MPR30750 1001 IF (VECTOR( NOVAR,NOVAR)) 1002,1002,1010 MPR30751 1002 NSTPM1 = NOSTEP - 1 MPR30752 PRINT 1004, NSTPM1 MPR30753 WRITE OUTPUT TAPE 6, 1004, NSTPM1 MPR30754 GO TO 1381 MPR30755 1010 SIGY = SIGMA(NOVAR) * SQRTF (VECTOR(NOVAR,NOVAR)/ DEFR) MPR30756 1011 IF(1.0 - VECTOR(NOVAR, NOVAR)) 1013,1013,1012 MPR30757 1012 R = SQRTF(1.0 - VECTOR(NOVAR, NOVAR) ) MPR30758 GO TO 1015 MPR30759 1013 R = 0.0 MPR30760 1015 DEFR=DEFR-1.0 MPR30761 1016 IF(DEFR)1017,1020,1020 MPR30762 1017 PRINT 1019, NOSTEP MPR30763 1018 WRITE OUTPUT TAPE 6, 1019, NOSTEP MPR30764 GO TO 1381 MPR30765 1020 VMIN = 0.0 MPR30766 1030 VMAX = 0.0 MPR30767 1035 NOIN = 0 MPR30768 1040 DO 1050 I = 1,NOVMI MPR30769 1041 IF (VECTOR (I,I)) 1042,1050,1060 MPR30770 1042 PRINT 1044, I, NOSTEP MPR30771 2043 WRITE OUTPUT TAPE 6, 1044, I, NOSTEP MPR30772 1045 GO TO 920 MPR30773 1060 IF(VECTOR(I,I) - TOL) 1050,1080,1080 MPR30774 1080 VAR = VECTOR(I,NOVAR) * VECTOR(NOVAR,I) / VECTOR(I,I) MPR30775 1090 IF(VAR)1100,1050,1110 MPR30776 1100 NOIN = NOIN + 1 MPR30777 1120 INDEX(NOIN) = I MPR30778 1130 COEN(NOIN) = VECTOR(I,NOVAR) * SIGMA(NOVAR) / SIGMA (I) MPR30779 1140 SIGMCO(NOIN) = (SIGY / SIGMA(I)) * SQRTF(VECTOR(I,I)) MPR30780 1150 IF (VMIN) 1160,1170,904 MPR30781 904 WRITE OUTPUT TAPE 6, 906 MPR30782 PRINT 906 MPR30783 GO TO 910 MPR30784 1170 VMIN = VAR MPR30785 1180 NOMIN = I MPR30786 1190 GO TO 1050 MPR30787 1160 IF(VAR - VMIN)1050,1050,1170 MPR30788 1110 IF (VAR - VMAX)1050,1050,1210 MPR30789 1210 VMAX = VAR MPR30790 1220 NOMAX = I MPR30791 1 SUBROUTINE CALC 12/01/64 PAGE 2 1050 CONTINUE MPR30792 1230 IF (NOIN) 903,1240,1245 MPR30793 903 WRITE OUTPUT TAPE 6, 907 MPR30794 PRINT 907 MPR30795 GO TO 910 MPR30796 1240 WRITE OUTPUT TAPE 6, 65, SIGY MPR30797 1260 GO TO 1350 MPR30798 1245 IF (IFCNST) 900,1250,1246 MPR30799 1246 CNST = 0.0 MPR30800 1247 GO TO 1300 MPR30801 1250 CNST = AVE(NOVAR) MPR30802 1270 DO 1280 I = 1,NOIN MPR30803 1290 J = INDEX(I) MPR30804 1280 CNST = CNST - (COEN(I) * AVE(J)) MPR30805 1300 IF(IFSTEP) 900,1310,1320 MPR30806 1310 IF (NOENT) 1311,1311,1313 MPR30807 1311 WRITE OUTPUT TAPE 6, 91,NOSTEP,K MPR30808 1312 GO TO 1314 MPR30809 1313 WRITE OUTPUT TAPE 6, 92,NOSTEP,K MPR30810 1314 WRITE OUTPUT TAPE 6, 70, FLEVEL, SIGY, R, CNST, MPR30811 1 (INDEX(J),COEN(J),SIGMCO(J), J = 1, NOIN ) MPR30812 1315 GO TO (1316,1580),NIX MPR30813 1316 GO TO NUMBER, (1320,1580) MPR30814 1320 FLEVEL = VMIN *(DEFR+1.0)/VECTOR(NOVAR,NOVAR) MPR30815 1330 IF(EFOUT + FLEVEL) 1350, 1350, 1340 MPR30816 1340 K = NOMIN MPR30817 1345 NOENT = 0 MPR30818 GO TO 1391 MPR30819 1350 DENOM = VECTOR(NOVAR,NOVAR) - VMAX MPR30820 IF(DENOM) 1351,1351,1352 MPR30821 1351 NIX = 2 MPR30822 GO TO 1370 MPR30823 1352 FLEVEL = VMAX * DEFR / DENOM MPR30824 1360 IF (EFIN - FLEVEL) 1370,1361,1380 MPR30825 1361 IF (EFIN) 1380,1380,1370 MPR30826 1370 K = NOMAX MPR30827 1390 NOENT = K MPR30828 1391 IF(K) 1392,1392,1400 MPR30829 1392 WRITE OUTPUT TAPE 6, 1395, NOSTEP MPR30830 1393 PRINT 1395, NOSTEP MPR30831 1394 GO TO 910 MPR30832 1400 DO 1410 I = 1,NOVAR MPR30833 1420 IF (I-K) 1430,1410,1430 MPR30834 1430 DO 1440 J = 1, NOVAR MPR30835 1450 IF (J-K) 1460,1440,1460 MPR30836 1460 VECTOR(I,J) = VECTOR(I,J) - (VECTOR(I,K) * VECTOR (K,J) / VECTOR MPR30837 -(K,K)) MPR30838 1440 CONTINUE MPR30839 1410 CONTINUE MPR30840 1470 DO 1480 I = 1, NOVAR MPR30841 1490 IF (I-K) 1500,1480,1500 MPR30842 1500 VECTOR (I,K) = - VECTOR (I,K) / VECTOR (K,K) MPR30843 1480 CONTINUE MPR30844 1510 DO 1520 J = 1, NOVAR MPR30845 1530 IF (J-K) 1540,1520,1540 MPR30846 1540 VECTOR(K,J) = VECTOR (K,J) / VECTOR (K,K) MPR30847 1 SUBROUTINE CALC 12/01/64 PAGE 3 1520 CONTINUE MPR30848 1550 VECTOR(K,K) = 1.0 / VECTOR(K,K) MPR30849 1560 GO TO (1000,1561),NIX MPR30850 1561 WRITE OUTPUT TAPE 6, 1004, NOSTEP MPR30851 R = 1.00 MPR30852 SIGY = 0.0 MPR30853 IFSTEP = 0 MPR30854 GO TO 1015 MPR30855 1380 PRINT 75, NOSTEP MPR30856 1381 IF (IFSTEP) 900, 1580,1570 MPR30857 1570 ASSIGN 1580 TO NUMBER MPR30858 1571 GO TO 1310 MPR30859 1580 RETURN MPR30860 C MPR30861 C ERROR RETURNS MPR30862 900 PRINT 905 MPR30863 910 REWIND 3 MPR30864 CALL EXIT MPR30865 1800 WRITE OUTPUT TAPE 6, 93, INVAR MPR30866 GO TO 910 MPR30867 C MPR30868 920 IF(IFSTEP)900,910,921 MPR30869 921 NIX=3 MPR30870 GO TO 1310 MPR30871 65 FORMAT (25H0 STANDARD ERROR OF Y = F12.6 ) MPR30872 70 FORMAT (11H F LEVEL F12.4 /25H STANDARD ERROR OF Y F12.4/MPR30873 135H MULTIPLE CORRELATION COEFFICIENT F12.5/ 12HMPR30874 2 CONSTANT F13.5/56H VARIABLE COEFFICIENT STD ERRMPR30875 3OR OF COEF // (16H X-I3,F15.5, F15.5 )) MPR30876 75 FORMAT (10H COMPLETED I5,20H STEPS OF REGRESSION) MPR30877 91 FORMAT (9H0STEP NO.I5 /19H VARIABLE REMOVED I8) MPR30878 92 FORMAT (9H0STEP NO.I5 /20H VARIABLE ENTERING I8) MPR30879 93 FORMAT (45H1INVAR SHOULD NOT BE GREATER THAN 150. IT IS I5 ) MPR30880 906 FORMAT (25H ERROR, VMIN PLUS. SOLONG) MPR30881 907 FORMAT (26H ERROR,NOIN MINUS, SOLONG ) MPR30882 1004 FORMAT (1H037HY SQUARE NON-POSITIVE,TERMINATE STEP I 5) MPR30883 1019 FORMAT (1H029H NO MORE DEGREES FREEDOM STEP I 5 ) MPR30884 1044 FORMAT (1H010H SQUARE X-I5,17H NEGATIVE. SOLONG I5,6H STEPS) MPR30885 1395 FORMAT (12H K=0. STEP I6, 7H SOLONG) MPR30886 905 FORMAT (42H ERROR IN CONTROL CARD, PROBLEM TERMINATED) MPR30887 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 SUBROUTINE CALC 12/01/64 PAGE 4 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 959 01677 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 958 01676 B 948 01664 C 938 01652 CTEMP 873 01551 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | DENOM 863 01537 INVAR 862 01536 I 861 01535 J 860 01534 NSTPM1 859 01533 | NUMBER 858 01532 R 857 01531 SIGY 856 01530 VMAX 855 01527 VMIN 854 01526 | 0 SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT STATEMENTS 0 EFN LOC EFN LOC EFN LOC EFN LOC EFN LOC | 8)21 65 01500 8)26 70 01471 8)2B 75 01417 8)2R 91 01407 8)2S 92 01377 | 8)2T 93 01367 8)S9 905 01301 8)SA 906 01355 8)SB 907 01347 8)VC 1004 01341 | 8)VR 1019 01330 8)10K 1044 01320 8)1BJ 1395 01307 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 1) 833 01501 2) 684 01254 3) 690 01262 5) 682 01252 6) 692 01264 | A)100 656 01220 A)107 669 01235 C)G0 837 01505 C)G1 838 01506 C)G2 839 01507 | C)G3 840 01510 C)G4 841 01511 C)100 842 01512 C)101 843 01513 C)102 844 01514 | C)103 845 01515 C)104 846 01516 C)105 847 01517 C)106 848 01520 C)107 849 01521 C)200 850 01522 C)201 851 01523 C)202 852 01524 C)203 853 01525 D)201 28 00034 | D)20P 284 00434 D)20T 310 00466 D)216 391 00607 D)21G 484 00744 D)21I 499 00763 | D)220 602 01132 D)223 622 01156 D)224 631 01167 D)40J 246 00366 D)414 368 00560 | D)422 615 01147 D)601 27 00033 D)60T 309 00465 D)614 367 00557 D)616 390 00606 | D)620 601 01131 D)622 614 01146 D)623 621 01155 D)624 630 01166 D)806 86 00126 D)80S 303 00457 D)81M 530 01022 D)A06 85 00125 D)A0S 302 00456 D)A1G 483 00743 | D)C0J 245 00365 D)E01 26 00032 D)E0T 308 00464 D)E16 389 00605 D)E1G 482 00742 | D)E20 600 01130 D)E23 620 01154 D)E24 629 01165 E)3 58 00072 E)A 144 00220 | E)D 183 00267 E)U 314 00472 E)13 365 00555 E)1E 468 00724 E)1I 497 00761 | E)21 607 01137 E)4J 248 00370 E)5O 540 01034 E)5S 561 01061 E)62 612 01144 E)8J 242 00362 E)GJ 242 00362 E)KJ 248 00370 | 1 SUBROUTINE CALC 12/01/64 PAGE 5 | | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT EXIT 5 00005 SQRT 3 00003 (FIL) 1 00001 (RWT) 4 00004 (SPH) 0 00000 | (STH) 2 00002 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | | EXIT SQRT (FIL) (RWT) (SPH) (STH) 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 781 19 00027 1000 20 00035 1001 21 00040 1002 22 00043 1010 28 00074 | 1011 29 00107 1012 30 00114 1013 32 00123 1015 33 00127 1016 34 00132 | 1017 35 00135 1018 37 00146 1020 40 00161 1030 41 00163 1035 42 00165 1040 43 00171 1041 44 00215 1042 45 00222 2043 47 00235 1045 49 00251 | 1060 50 00252 1080 51 00257 1090 52 00264 1100 53 00271 1120 54 00276 | 1130 55 00300 1140 56 00306 1150 57 00320 904 58 00324 1170 61 00336 | 1180 62 00340 1190 63 00342 1160 64 00343 1110 65 00350 1210 66 00355 | 1220 67 00357 1050 68 00372 1230 69 00403 903 70 00406 1240 73 00416 1260 75 00424 1245 76 00425 1246 77 00431 1247 78 00433 1250 79 00435 | 1270 80 00437 1290 81 00442 1280 82 00446 1300 83 00460 1310 84 00467 | 1311 85 00473 1312 87 00507 1313 88 00510 1314 90 00524 1315 96 00553 | 1316 97 00556 1320 98 00561 1330 99 00567 1340 100 00574 1345 101 00602 | 1350 103 00610 1351 105 00616 1352 107 00623 1360 108 00627 1361 109 00634 1370 110 00640 1390 111 00646 1391 112 00650 1392 113 00653 1393 115 00664 | 1394 117 00674 1400 118 00675 1420 119 00727 1430 120 00732 1450 121 00745 | 1460 122 00750 1440 123 00764 1410 124 00773 1470 125 01002 1490 126 01023 | 1500 127 01026 1480 128 01036 1510 129 01042 1530 130 01052 1540 131 01055 | 1520 132 01063 1550 133 01067 1560 134 01073 1561 135 01076 1380 141 01117 1381 143 01133 1570 144 01141 1571 145 01143 1580 146 01150 900 148 01157 | 910 149 01170 1800 151 01175 920 154 01210 921 155 01214 | 1 1 SUBROUTINE NORM 12/01/64 PAGE 1 SUBROUTINE NORM MPR30889 C MPR30890 EQUIVALENCE (DATA,X) MPR30891 COMMON DATA, AVE, SIGMA, COEN, SIGMCO, INDEX, Z, W, X, Y, RUN,MPR30892 1 WHT, TOL, EFIN, EFOUT, IFSTEP, IFAVE, IFRESD, IFCOEN, IFCNST, MPR30893 2NOIN, VAR, K, FLEVEL, NOENT, NOMIN, NOMAX, NOVAR, NIX, NOVMI, MPR30894 3NOVPL, CNST, NOSTEP, DEFR, VECTOR, NVP1 MPR30895 4,L1,L2,L3,L4,L5,L6,L7,L8,L9,C1,C2,NN1 MPR30896 5, NAMVAR, NAMEOP MPR30897 DIMENSION DATA(130), AVE(130), SIGMA(130), COEN(130), SIGMCO(130)MPR30898 1, INDEX(130), Z(130), W(70), X(130), Y(20), VECTOR(131,131), MPR30899 2L1(26), L2(26), L3(26), L4(26), L5(26), L6(26), L7(26),L8(26), MPR30900 3L9(25),A(10),B(10), C(65), C1(26),C2(26),CTEMP(10) MPR30901 4,NAMEOP(8),NAMVAR(4) MPR30902 C MPR30903 780 NOSTEP = -1 MPR30904 782 DEFR = VECTOR(NOVPL,NOVPL) - 1.0 MPR30905 790 DO 800 I = 1,NOVAR MPR30906 791 IF(VECTOR(I,I)) 792,794,810 MPR30907 792 PRINT 793, I MPR30908 GO TO 910 MPR30909 794 WRITE OUTPUT TAPE 6, 795, I MPR30910 796 SIGMA(I) = 1.0 MPR30911 797 GO TO 800 MPR30912 810 SIGMA(I) = SQRTF (VECTOR (I,I)) MPR30913 800 VECTOR(I,I) = 1.0 MPR30914 820 DO 830 I = 1,NOVMI MPR30915 840 IP1 = I + 1 MPR30916 841 DO 830 J = IP1, NOVAR MPR30917 850 VECTOR(I,J) = VECTOR(I,J) /( SIGMA(I)* SIGMA(J)) MPR30918 830 VECTOR(J,I) = VECTOR(I,J) MPR30919 860 IF (IFCOEN) 900, 870, 1000 MPR30920 870 WRITE OUTPUT TAPE 6, 60 MPR30921 874 NOVM2 = NOVMI - 1 MPR30922 875 DO 885 I = 1, NOVM2 MPR30923 880 IP1 = I + 1 MPR30924 885 WRITE OUTPUT TAPE 6, 35,(I,J,VECTOR(I,J),J=IP1,NOVMI ) MPR30925 890 WRITE OUTPUT TAPE 6, 40, (I,VECTOR(I,NOVAR),I=1,NOVMI) MPR30926 1000 RETURN MPR30927 C MPR30928 C ERROR RETURNS MPR30929 900 PRINT 905 MPR30930 910 REWIND 3 MPR30931 CALL EXIT MPR30932 C MPR30933 35 FORMAT (1H 7H X(I2,7H) VS X(I2,3H) = F15.6, MPR30934 1 6H X(I2,7H) VS X(I2,3H) = F15.6, MPR30935 2 6H X(I2,7H) VS X(I2,3H) = F15.6 ) MPR30936 40 FORMAT (1H 7H X(I2,12H) VS Y =F15.6, MPR30937 1 6H X(I2,12H) VS Y =F15.6, MPR30938 2 6H X(I2,12H) VS Y =F15.6 ) MPR30939 60 FORMAT(1H069H SIMPLE CORRELATIMPR30940 -ON COEFFICIENTS//) MPR30941 793 FORMAT (31H ERROR RESIDUAL SQUARE VARIABLE I4,31H IS NEGATIVE,PROBMPR30942 1LEM TERMINATED ) MPR30943 795 FORMAT ( 10H VARIABLE I5,13H IS CONSTANT ) MPR30944 1 SUBROUTINE NORM 12/01/64 PAGE 2 905 FORMAT (42H ERROR IN CONTROL CARD, PROBLEM TERMINATED) MPR30945 END(1,0,0,0,0,0,1,0,0,1,0,0,0,0,0) 1 SUBROUTINE NORM 12/01/64 PAGE 3 | | 0 STORAGE NOT USED BY PROGRAM | 0 DEC OCT DEC OCT | 454 00706 14077 33375 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | AVE 32431 77257 C1 14142 33476 C2 14116 33444 CNST 31540 75464 COEN 32171 76653 | DATA 32561 77461 DEFR 31538 75462 EFIN 31558 75506 EFOUT 31557 75505 FLEVEL 31548 75474 | IFAVE 31555 75503 IFCNST 31552 75500 IFCOEN 31553 75501 IFRESD 31554 75502 IFSTEP 31556 75504 INDEX 31911 76247 K 31549 75475 L1 14375 34047 L2 14349 34015 L3 14323 33763 | L4 14297 33731 L5 14271 33677 L6 14245 33645 L7 14219 33613 L8 14193 33561 | L9 14167 33527 NAMEOP 14085 33405 NAMVAR 14089 33411 NIX 31543 75467 NN1 14090 33412 | NOENT 31547 75473 NOIN 31551 75477 NOMAX 31545 75471 NOMIN 31546 75472 NOSTEP 31539 75463 | NOVAR 31544 75470 NOVMI 31542 75466 NOVPL 31541 75465 NVP1 14376 34050 RUN 31561 75511 SIGMA 32301 77055 SIGMCO 32041 76451 TOL 31559 75507 VAR 31550 75476 VECTOR 31537 75461 | WHT 31560 75510 W 31651 75643 X 32561 77461 Y 31581 75535 Z 31781 76045 | 0 STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | A 453 00705 B 443 00673 C 433 00661 CTEMP 368 00560 0 STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENCE STATEMENT | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | IP1 358 00546 I 357 00545 J 356 00544 NOVM2 355 00543 | 0 SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT STATEMENTS | 0 EFN LOC EFN LOC EFN LOC EFN LOC EFN LOC 8)13 35 00530 8)18 40 00473 8)1S 60 00437 8)OP 793 00421 8)OR 795 00404 | 8)S9 905 00375 | 0 LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT | 1) 345 00531 2) 235 00353 3) 239 00357 6) 240 00360 A)100 222 00336 C)G0 349 00535 C)100 350 00536 C)101 351 00537 C)102 352 00540 C)103 353 00541 | C)201 354 00542 D)208 99 00143 D)20K 207 00317 D)408 101 00145 D)60K 206 00316 | D)E0K 205 00315 E)3 41 00051 | 0 LOCATIONS OF NAMES IN TRANSFER VECTOR | 0 DEC OCT DEC OCT DEC OCT DEC OCT DEC OCT EXIT 5 00005 SQRT 3 00003 (FIL) 1 00001 (RWT) 4 00004 (SPH) 0 00000 | (STH) 2 00002 | 0 ENTRY POINTS TO SUBROUTINES NOT OUTPUT FROM LIBRARY | | EXIT SQRT (FIL) (RWT) (SPH) (STH) 1 SUBROUTINE NORM 12/01/64 PAGE 4 | | 0 EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS | 0 EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC EFN IFN LOC | 780 12 00021 782 13 00023 790 14 00026 791 15 00034 792 16 00037 794 19 00053 796 21 00065 797 22 00067 810 23 00070 800 24 00075 | 820 25 00103 840 26 00111 841 27 00114 850 28 00146 830 29 00154 | 860 30 00171 870 31 00175 874 32 00201 875 33 00204 880 34 00212 | 885 35 00215 890 40 00254 1000 45 00306 900 47 00312 910 48 00320 | 1 1 G2 ER MPR3 09-04-64 3145 SY 0037 FAP SY + DATA 12/01/64 PAGE 1 G2 ER MPR3 09-04-64 3145 SY 0037 FAP SY + DATA * FAP MPR30947 1 DUMMY ENTRY LOCATIONS FOR UNNECESSARY SUBROUTINES 12/01/64 PAGE 1 | | 00000 ENTRY FCN1 MPR30949 | 00000 ENTRY FCN2 MPR30950 | 00000 ENTRY FCN3 MPR30951 00001 ENTRY OUTPUT MPR30952 | 00000 0020 00 4 00004 FCN1 TRA 4,4 MPR30954 | 00001 0020 00 4 00001 OUTPUT TRA 1,4 MPR30955 | 00000 FCN2 SYN FCN1 MPR30956 | 00000 FCN3 SYN FCN1 MPR30957 END MPR30958 | 1 DUMMY ENTRY LOCATIONS FOR UNNECESSARY SUBROUTINES 12/01/64 PAGE 1 | POST PROCESSOR ASSEMBLY DATA | | 2 IS THE FIRST LOCATION NOT USED BY THIS PROGRAM | 0REFERENCES TO DEFINED SYMBOLS | 0 FCN1 0, 2 | 0 FCN2 0, 2 | 0 FCN3 0, 2 1 OUTPUT 0 | | 0 NO ERROR IN ABOVE ASSEMBLY. | 1 * DATA 0LIBRARY ENTRY POINTS, (FPT) (TSHM) (RTN) (RWT) (STHM) (FIL) (SPH) (TSB) (RLR) EXIT LOG EXP ATAN EXP(3 (STB) (WLR) SQRT 0LOGICAL MACHINE TOTAL TOTAL NOISE RECORDS TOTAL REDUNDANCIES POSITIONING TAPE TAPE WRITES READS WRITING READING WRITING READING ERRORS 1 A1 0 31 0 0 0 0 0 2 A4 1360 1386 0 0 0 0 0 3 B4 224 245 0 0 0 0 0 4 A3 893 935 0 0 0 0 0 5 A2 0 980 0 0 0 0 0 6 B1 1214 11 0 0 0 0 0 7 B2 241 238 0 0 0 0 0 0 EXECUTION 1 VARIABLE TRANSFORMATIONS ( Z( 1) + -0. ) - (00( -0) + 96.250000) = W( 1) ( Z( 2) + -0. ) - (00( -0) + 9.000000) = W( 2) ( W( 1) + -0. ) + (00( -0) + -0. ) = X( 1) ( W( 2) + -0. ) + (00( -0) + 0. ) = X( 2) ( W( 2) + -0. ) * ( W( 2) + -0. ) = X( 3) ( W( 1) + -0. ) * ( W( 2) + -0. ) = X( 4) ( Z( 3) + -0. ) + (00( -0) + 0. ) = Y( 1) 1STEPWISE REGRESSION PROBLEM NO 1 NO OF DATA = 16 NO OF VARIABLES = 5 WEIGHTED DEGREES OF FREEDOM = 16.00 F LEVEL TO ENTER VARIABLE = 0.001 F LEVEL TO REMOVE VARIABLE = 0.000 SUM OF VARIABLES SUM X( 1) = -0. SUM X( 2) = -0. SUM X( 3) = 80.0000 SUM X( 4) = -0. SUM Y = 1483.3199 0 RAW SUM OF SQUARES AND CROSS PRODUCTS X( 1) VS X( 1) = 107.000000 X( 1) VS X( 2) = 0. X( 1) VS X( 3) = 0. X( 1) VS X( 4) = 0. X( 2) VS X( 2) = 80.000000 X( 2) VS X( 3) = 0. X( 2) VS X( 4) = 0. X( 3) VS X( 3) = 656.000000 X( 3) VS X( 4) = 0. X( 4) VS X( 4) = 535.000000 X( X( 1) VS Y = 55.410030 X( 2) VS Y = -4.600002 X( 3) VS Y = 7320.599670 X( 4) VS Y = 1.150040 X( Y VS Y = 137579.843750 0 AVERAGE VALUE OF VARIABLES SUM X( 1) = -0. SUM X( 2) = -0. SUM X( 3) = 5.0000 SUM X( 4) = -0. SUM Y = 92.7075 0 RESIDUAL SUMS OF SQUARES AND CROSS PRODUCTS X( 1) VS X( 1) = 107.000000 X( 1) VS X( 2) = -0. X( 1) VS X( 3) = 0. X( 1) VS X( 4) = -0. X( 2) VS X( 2) = 80.000000 X( 2) VS X( 3) = 0. X( 2) VS X( 4) = -0. X( 3) VS X( 3) = 256.000000 X( 3) VS X( 4) = 0. X( 4) VS X( 4) = 535.000000 X( X( 1) VS Y = 55.410030 X( 2) VS Y = -4.600002 X( 3) VS Y = -95.999939 X( 4) VS Y = 1.150040 X( Y VS Y = 64.968750 0 SIMPLE CORRELATION COEFFICIENTS X( 1) VS X( 2) = -0. X( 1) VS X( 3) = 0. X( 1) VS X( 4) = -0. X( 2) VS X( 3) = 0. X( 2) VS X( 4) = -0. X( X( 3) VS X( 4) = 0. X( X( 1) VS Y = 0.664575 X( 2) VS Y = -0.063806 X( 3) VS Y = -0.744387 X( 4) VS Y = 0.006169 X( 0 STANDARD ERROR OF Y = 2.081165 0STEP NO. 1 VARIABLE ENTERING 3 F LEVEL 17.3980 STANDARD ERROR OF Y 1.4385 MULTIPLE CORRELATION COEFFICIENT 0.74439 CONSTANT 94.58249 VARIABLE COEFFICIENT STD ERROR OF COEF X- 3 -0.37500 0.08990 0STEP NO. 2 VARIABLE ENTERING 1 F LEVEL 1358.0820 STANDARD ERROR OF Y 0.1454 MULTIPLE CORRELATION COEFFICIENT 0.99788 CONSTANT 94.58249 VARIABLE COEFFICIENT STD ERROR OF COEF X- 1 0.51785 0.01405 X- 3 -0.37500 0.00908 0STEP NO. 3 VARIABLE ENTERING 2 F LEVEL 312.1196 STANDARD ERROR OF Y 0.0291 MULTIPLE CORRELATION COEFFICIENT 0.99992 CONSTANT 94.58249 VARIABLE COEFFICIENT STD ERROR OF COEF X- 1 0.51785 0.00281 X- 2 -0.05750 0.00325 X- 3 -0.37500 0.00182 0STEP NO. 4 VARIABLE ENTERING 4 F LEVEL 3.5330 STANDARD ERROR OF Y 0.0265 MULTIPLE CORRELATION COEFFICIENT 0.99994 CONSTANT 94.58249 VARIABLE COEFFICIENT STD ERROR OF COEF X- 1 0.51785 0.00256 X- 2 -0.05750 0.00296 X- 3 -0.37500 0.00165 X- 4 0.00215 0.00114 0 DIAGONAL ELEMENTS VAR.NO. VALUE 1 1.000000 2 1.000000 3 1.000000 4 1.000000 0 PREDICTED VS ACTUAL RESULTS RUN NO. ACTUAL PREDICTED DEVIATION 16. 94.65000 94.65177 -0.00177 16. 94.54000 94.53999 0.00001 16. 91.43000 91.42822 0.00178 16. 93.62000 93.62037 -0.00037 16. 93.50000 93.49999 0.00001 16. 90.38000 90.37962 0.00038 16. 92.59000 92.58896 0.00104 16. 92.46000 92.45999 0.00001 6. 89.33000 89.33102 -0.00102 6. 91.76000 91.76355 -0.00355 6. 90.74000 90.74074 -0.00074 6. 89.72000 89.71794 0.00206 6. 93.30000 93.29775 0.00225 6. 96.20000 96.19887 0.00113 6. 96.10000 96.09999 0.00001 6. 93.00000 93.00112 -0.00112 1 TRANSFORMED VARABLES (X) RUN NO. 1 2 3 4 5 16. 0.750000 -1.000000 1.000000 -0.750000 94.650000 16. 0.750000 1.000000 1.000000 0.750000 94.540000 16. 0.750000 3.000000 9.000000 2.250000 91.429999 16. -1.250000 -1.000000 1.000000 1.250000 93.620000 16. -1.250000 1.000000 1.000000 -1.250000 93.500000 16. -1.250000 3.000000 9.000000 -3.750000 90.379999 16. -3.250000 -1.000000 1.000000 3.250000 92.589999 16. -3.250000 1.000000 1.000000 -3.250000 92.459999 6. -3.250000 3.000000 9.000000 -9.750000 89.330000 6. 0.750000 -3.000000 9.000000 -2.250000 91.759999 6. -1.250000 -3.000000 9.000000 3.750000 90.740000 6. -3.250000 -3.000000 9.000000 9.750000 89.719999 6. 3.750000 -3.000000 9.000000 -11.250000 93.299999 6. 3.750000 -1.000000 1.000000 -3.750000 96.200000 6. 3.750000 1.000000 1.000000 3.750000 96.099999 6. 3.750000 3.000000 9.000000 11.250000 93.000000 1 2042 LINES OUTPUT THIS JOB. 1 0 FORTRAN MONITOR RETURNING TO IBSYS $STOP | 1 PERIPHERAL UNIT POSITIONS AT END OF JOBS | 1 SYSPP1 IS B2 REC. 00000, FILE 00002 SYSOU1 IS B1 REC. 01359, FILE 00000 SYSIN1 IS A2 REC. 00002, FILE 00001| 1 END OF JOBS | 1