| 1 $JOB SPFAC $EXECUTE IBJOB IBJOB VERSION 5 HAS CONTROL. | $IBJOB GO,MAP | $IBFTC CF1 1 12/01/64 PAGE 1 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE SPFAC(G,H,NP,IEND,EPS,IFLAG) SPFAC001 | C SPFAC002 C THIS IS A SUBROUTINE TO PERFORM THE FACTORIZATION SPFAC003 | C G(S) = H(-S)*H(S), WHERE G(S) = SUM(I=0 TO NP) G(I)*S**I IS A SPFAC004 | C REAL POLYNOMIAL WITH THE PROPERTIES SPFAC005 | C G(S) = G(-S) SPFAC006 | C G(JP) .GT. 0 FOR ALL REAL P, J = SQRT(-1) SPFAC007 C AND H(S) IS A POLYNOMIAL WITH ALL ITS ROOTS IN THE LEFT HAND SPFAC008 | C PLANE. THIS PROBLEM ARISES IN THE SOLUTION OF WIENER-HOPF SPFAC009 | C EQUATIONS AND IN CIRCUIT SYNTHESIS. THE ROOTS OF THE POLYNOMIAL, SPFAC010 | C BEING REDUNDANT, ARE NOT CALCULATED. THE PROCEDURE FOLLOWS...... SPFAC011 | C (1) THE LEFT-HAND S PLANE IS MAPPED INTO THE UNIT CIRCLE UNDER THESPFAC012 C TRANSFORMATION S = (Z-1)/(Z+1). THUS SPFAC013 | C G(Z-1/Z+1) = SUM(K=-M TO M) C(K)*Z**-K /((1+Z)*(1+1/Z))**M, M=NP/2SPFAC014 | C WHERE THE C(K) ARE COMPUTED BY MULTIPLYING OUT THE TRANSFORMATION.SPFAC015 | C (2) THE SET OF NON-LINEAR ALGEBRAIC EQUATIONS SPFAC016 | C SPFAC017 C X(I) = C(I) + Y(0,I) I = 0,1,...,M SPFAC018 | C SPFAC019 | C Y(I-1,J-1) = Y(I,J) - X(I)*X(J)/X(0) I=1,2,...,M, J=I,I+1,...,M SPFAC020 | C SPFAC021 | C Y(I,J) = Y(J,I) SPFAC022 C SPFAC023 | C Y(M,I) = 0 FOR ALL I SPFAC024 | C SPFAC025 | C IS SOLVED RECURSIVELY WITH THE INITIAL ITERATES SPFAC026 | C SPFAC027 C Y(I,J) = 0, I,J=0,1,...,M-1 SPFAC028 | C SPFAC029 | C X(I) = C(I), I=0,1,...,M SPFAC030 | C SPFAC031 | C UNTIL THE SUBWAY (ABSOLUTE VALUE) NORM OF THE DIFFERENCE BETWEEN SPFAC032 C SUCCESSIVE ITERATES OF THE MATRIX Y IS LESS THAN EPS OR UNTIL SPFAC033 | C IEND ITERATIONS HAVE BEEN CALCULATED. SPFAC034 | C (3) FROM THE SOLUTION OF (2) THE QUANTITIES SPFAC035 | C SPFAC036 | C D(K) = X(K)/SQRT(X(0)) K=0,1,...,M SPFAC037 C SPFAC038 | C ARE CALCULATED. DEFINE SPFAC039 | C D(Z) = SUM(K=0 TO M) D(K)*Z**-K. THEN IT CAN BE SHOWN THAT SPFAC040 | C G(Z-1/Z+1) = D(1/Z)/(1+Z)**M * D(Z)/(1+1/Z)**M SPFAC041 | C (4) THE SOLUTION H(S) IS OBTAINED BY INVERTING THE TRANSFORMATION.SPFAC042 C H(S) = D(1+S/1-S) * ((1+S)/2)**M SPFAC043 | C IS THE REQUIRED STABLE FACTOR. SPFAC044 | C SPFAC045 | C SPFAC046 | C USE OF THE SUBROUTINE SPFAC047 C SPFAC048 | C SPFAC049 | C THE ARRAY G CONTAINS THE COEFFICIENTS OF THE POLYNOMIAL G(S) SPFAC050 | C ARRANGED IN ASCENDING POWERS OF S. H IS RETURNED WITH THE STABLE SPFAC051 | C FACTOR, ALSO IN ASCENDING POWERS OF S. NP IS THE HIGHEST POWER OFSPFAC052 C S IN G(S). SPFAC053 | C IEND IS THE MAXIMUM NUMBER OF ITERATIONS TO BE PERFORMED. IF THE SPFAC054 | C ITERATION CONVERGES PRIOR TO IEND, IEND IS SET THE THE NUMBER OF SPFAC055 | 1 12/01/64 PAGE 2 | CF1 - EFN SOURCE STATEMENT - IFN(S) - 0C STEPS REQUIRED. IF THE ITERATION FAILS TO CONVERGE IN IEND STEPS,SPFAC056 | C IEND IS SET TO ZERO. IN EITHER CASE, IFLAG=0 AND H CONTAINS SPFAC057 | C THE LATEST APPROXIMATE TO THE FACTOR. IF AT ANY STEP X(M+1) SPFAC058 | C BECOMES NEGATIVE, IFLAG IS SET TO 1 AND IEND IS SET TO THE SPFAC059 | C NEGATIVE OF THE CURRENT ITERATION NUMBER. IN THE LATTER CASE, H SPFAC060 C IS UNDEFINED. SPFAC061 | C EPS IS THE NORM CONVERGENCE CRITERION. THE USER MAY SET HIS OWN, SPFAC062 | C OR MAY LET THE PROGRAM SELECT EPS BY ENTERING THE SUBROUTINE WITH SPFAC063 | C EPS=0. THE PROGRAM WILL CALCULATE EPS = C(0)*10.**-(INX-1), WHERE SPFAC064 | C INX IS THE NUMBER OF DECIMAL DIGITS IN THE FLOATING POINT MANTISSASPFAC065 C THIS HELPS INSURE THAT EPS IS LARGER THAN THE TRUNCATION NOISE SPFAC066 | C LEVEL ENCOUNTERED IN COMPUTING THE ERROR NORM. SPFAC067 | C INX IS CALCULATED BY THE PROGRAM. SPFAC068 | C THE ARRAYS IN THE SUBROUTINE MUST HAVE AT LEAST THE FOLLOWINGSPFAC069 | C DIMENSIONS SPFAC070 C DIMENSION G(NP+1), H(M+1) SPFAC071 | C DIMENSION WORK(NP+1),C(M+1),JH(NP+1),JTEM(NP+1),JWORK(NP+1) SPFAC072 | C DIMENSION Y(M,M), X(M+1), D(M+1) SPFAC073 | C SINCE INTEGER ARITHMETIC IS USED IN THE TRANSFORMATION, SPFAC074 | C NP MUST NOT EXCEED THE BINARY WORD LENGTH OF THE MACHINE ON WHICH SPFAC075 C THE PROGRAM IS RUN. E.G. ON THE 7094, NP MUST NOT EXCEED 36. SPFAC076 | C THIS VERSION IS SET FOR A MAXIMUM VALUE OF NP=16. SPFAC077 | C SPFAC078 | C ******************************************************************SPFAC079 | C SPFAC080 DIMENSION G(17),H(9) SPFAC081 | DIMENSION WORK(17),C(9),JH(17),JTEM(17),JWORK(17) SPFAC082 | DIMENSION Y(8,8),X(9),D(9) SPFAC083 | EQUIVALENCE (WORK(1),X(1)), (C(1),D(1)) SPFAC084 | C SPFAC085 C SPFAC086 | C SPFAC087 | NPP1=NP+1 SPFAC088 | M=NP/2 SPFAC089 | MP1=M+1 SPFAC090 C SPFAC091 | C SPFAC092 | C PERFORM BILINEAR TRANSFORMATION S=(Z-1)/(Z+1) AND EXPRESS SPFAC093 | C RESULT AS A POLYNOMIAL IN POWERS OF (1/Z). SPFAC094 | C SPFAC095 DO 101 K=1,NPP1 SPFAC096 | JTEM(K)=0 SPFAC097 | 101 WORK(K) = 0. SPFAC098 | JTEM(1)=1 SPFAC099 | C SPFAC100 C GENERATE TRANSFORMATION POLYNOMIALS SPFAC101 | C (1+Z)**(I-J)*(-1+Z)**J SPFAC102 | C SPFAC103 | C OUTER LOOP WITH J SPFAC104 | C SPFAC105 DO 102 J=1,NPP1 SPFAC106 | DO 103 JJ=1,NPP1 SPFAC107 | 103 JH(JJ)=JTEM(JJ) SPFAC108 | K=NPP1-J+1 SPFAC109 | JM1=J-1 SPFAC110 C SPFAC111 | 1 12/01/64 PAGE 3 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0C MULTIPLY POLYNOMIAL BY COEFFICIENT OF G AND ADD TO WORK. SPFAC112 | C SPFAC113 DO 104 I=1,K SPFAC114 | IPJM1=I+JM1 SPFAC115 | IF(K-I) 106,105,106 SPFAC116 | 105 SAVE=G(J) SPFAC117 | DO 107 KK=1,NPP1 SPFAC118 FLJ=JH(KK) SPFAC119 | 107 WORK(KK) = FLJ*SAVE+WORK(KK) SPFAC120 | GO TO 104 SPFAC121 | C SPFAC122 | C INNER LOOP WITH I SPFAC123 C SPFAC124 | 106 DO 108 II=1,IPJM1 SPFAC125 | 108 JWORK(II)=JH(II) SPFAC126 | DO 109 II=1,IPJM1 SPFAC127 | 109 JH(II+1)=JWORK(II)+JH(II+1) SPFAC128 104 CONTINUE SPFAC129 | C SPFAC130 | C END OF INNER LOOP SPFAC131 | C SPFAC132 | IF(NPP1-J) 110,102,110 SPFAC133 110 CONTINUE SPFAC134 | DO 111 II=1,NPP1 SPFAC135 | 111 JWORK(II)=JTEM(II) SPFAC136 | DO 112 II=1,J SPFAC137 | 112 JTEM(II+1)=JWORK(II)-JTEM(II+1) SPFAC138 JTEM(1)=-JTEM(1) SPFAC139 | 102 CONTINUE SPFAC140 | C SPFAC141 | C END OF OUTER LOOP SPFAC142 | C SPFAC143 C WORK(K) CONTAINS A POLYNOMIAL IN Z, WITH WORK(I)=WORK(NP+2-I) SPFAC144 | C DIVIDE BY Z**NP AND RETAIN ONLY POWERS OF (1/Z) SPFAC145 | C SPFAC146 | JJ=MP1/2 SPFAC147 | IF(JJ) 114,114,113 SPFAC148 113 DO 115 II=1,JJ SPFAC149 | SAVE = WORK(II) SPFAC150 | K=MP1-(II-1) SPFAC151 | WORK(II) = WORK(K) SPFAC152 | 115 WORK(K) = SAVE SPFAC153 114 CONTINUE SPFAC154 | DO 116 II=MP1,NP SPFAC155 | 116 WORK(II+1) = 0. SPFAC156 | C SPFAC157 | C SPFAC158 C SPFAC159 | DO 2 I=1,MP1 SPFAC160 | 2 C(I) = WORK(I) SPFAC161 | C SPFAC162 | IF (EPS) 6,6,5 SPFAC163 C SPFAC164 | C CALCULATE ERROR CRITERION TO BE ABOVE TRUNCATION NOISE LEVEL SPFAC165 | C SPFAC166 | C DETERMINE EFFECTIVE MANTISSA WORD LENGTH SPFAC167 | 1 12/01/64 PAGE 4 CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0C SPFAC168 | 6 SAVE=0.5 SPFAC169 | DO 7 INX=1,256 SPFAC170 | SAVE=0.5*SAVE SPFAC171 TEMP=0.5+SAVE SPFAC172 | IF(TEMP-.5) 8,8,7 SPFAC173 | 7 CONTINUE SPFAC174 | INX=256 SPFAC175 | 8 SAVE=INX SPFAC176 SAVE=0.30103*SAVE SPFAC177 | INX=SAVE SPFAC178 | C SPFAC179 | C CALCULATE EPS SPFAC180 | C SPFAC181 EPS=C(1)*10.**(-INX+1) SPFAC182 | C SPFAC183 | C SOLVE THE SET OF NON-LINEAR ALGEBRAIC EQUATIONS BY ITERATION SPFAC184 | C SPFAC185 | C INITIALIZE THE ITERATION SPFAC186 C SPFAC187 113 | 5 DO 201 I=1,MP1 SPFAC188 | 201 X(I) = C(I) SPFAC189 | DO 202 I=1,M SPFAC190 | DO 202 J=I,M SPFAC191 202 Y(I,J)=0. SPFAC192 | IFLAG=0 SPFAC193 | C SPFAC194 | C IF X(0) IS NOT POSITIVE, ABORT THE ROUTINE SPFAC195 | C SPFAC196 IF(X(1)) 203,203,204 SPFAC197 | 203 JEND=-1 SPFAC198 | IFLAG=1 SPFAC199 | GO TO 205 SPFAC200 | 204 XINV=1.0/X(1) SPFAC201 JEND=1 SPFAC202 | C SPFAC203 | C GO INTO LOOP TO CALCULATE Y(I,J) SPFAC204 | C SPFAC205 | GO TO 206 SPFAC206 C SPFAC207 | C SPFAC208 | C ITERATION LOOP SPFAC209 | C SPFAC210 | C SPFAC211 207 JEND=JEND+1 SPFAC212 | IF(JEND-IEND) 209,209,208 SPFAC213 | C TERMINATE THE LOOP DUE TO POOR CONVERGENCE SPFAC214 | C SPFAC215 | 208 JEND=0 SPFAC216 GO TO 205 SPFAC217 | C CALCULATE X(I) SPFAC218 | 209 DO 212 I=1,M SPFAC219 | 212 X(I) = C(I) + Y(1,I) SPFAC220 | IF(X(1)) 210,210,211 SPFAC221 210 JEND = -JEND SPFAC222 | IFLAG=1 SPFAC223 | 1 12/01/64 PAGE 5 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 GO TO 205 SPFAC224 C SPFAC225 | 211 XINV = 1.0/X(1) SPFAC226 | C SPFAC227 | C CALCULATE Y(I,J) AND TEST FOR CONVERGENCE USING SUBWAY NORM OF SPFAC228 | C Y(I,J) - PREVIOUS Y(I,J). SPFAC229 C ONLY THE UPPER RIGHT HALF OF Y IS CALCULATED SPFAC230 | C SPFAC231 | 206 ENRM=-EPS SPFAC232 | DO 213 I=1,M SPFAC233 | SAVE=X(I+1)*XINV SPFAC234 DO 214 J=I,M SPFAC235 | SUM=-SAVE*X(J+1) SPFAC236 | IF(J-M) 215,216,216 SPFAC237 | 215 IF(I-M) 217,216,216 SPFAC238 | 217 SUM=SUM+Y(I+1,J+1) SPFAC239 216 TEMP = SUM-Y(I,J) SPFAC240 | Y(I,J) = SUM SPFAC241 | C COMPUTE ERROR NORM SPFAC242 | IF(TEMP) 218,219,220 SPFAC243 | 218 TEMP=-TEMP SPFAC244 220 IF(J-I) 221,222,221 SPFAC245 | 222 ENRM=ENRM+TEMP SPFAC246 | GO TO 219 SPFAC247 | 221 ENRM=ENRM+2.0*TEMP SPFAC248 | 219 CONTINUE SPFAC249 214 CONTINUE SPFAC250 | 213 CONTINUE SPFAC251 | C SPFAC252 | C CHECK FOR CONVERGENCE SPFAC253 | C SPFAC254 IF(ENRM) 205,205,207 SPFAC255 | C SPFAC256 | C SPFAC257 | 205 IEND=JEND SPFAC258 | C SPFAC259 C SPFAC260 | C SPFAC261 | IF(IFLAG) 1,3,1 SPFAC262 | 3 TEMP = 1.0/SQRT(X(1)) SPFAC263 194 | DO 4 I=1,MP1 SPFAC264 4 D(I) = X(I)*TEMP SPFAC265 | C SPFAC266 | C SPFAC267 | C SPFAC268 | C PERFORM THE INVERSE BILINEAR TRANSFORM SPFAC269 C Z=(1+S)/(1-S) SPFAC270 | C SPFAC271 | DO 301 I=1,MP1 SPFAC272 | JTEM(I)=0 SPFAC273 | 301 H(I)=0. SPFAC274 JTEM(1)=1 SPFAC275 | C SPFAC276 | C GENERATE TRANSFORMATION POLYNOMIALS SPFAC277 | C (1-S)**J*(1+S)**(I-J) SPFAC278 | C SPFAC279 1 12/01/64 PAGE 6 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0C OUTER LOOP WITH J SPFAC280 | C SPFAC281 | DO 302 J=1,MP1 SPFAC282 DO 303 K=1,MP1 SPFAC283 | 303 JH(K)=JTEM(K) SPFAC284 | JM1=J-1 SPFAC285 | K=MP1-JM1 SPFAC286 | C SPFAC287 C MULTIPLY POLYNOMIALS BY COEFFICIENTS OF D AND ADD TO H SPFAC288 | C SPFAC289 | DO 304 I=1,K SPFAC290 | IPJM1=I+JM1 SPFAC291 | IF(K-I) 306,305,306 SPFAC292 305 SAVE = D(J) SPFAC293 | DO 307 KK=1,IPJM1 SPFAC294 | FLJ=JH(KK) SPFAC295 | 307 H(KK)=FLJ*SAVE + H(KK) SPFAC296 | GO TO 304 SPFAC297 C SPFAC298 | C INNER LOOP WITH I SPFAC299 | C SPFAC300 | 306 DO 308 II=1,IPJM1 SPFAC301 | 308 JWORK(II)=JH(II) SPFAC302 DO 309 II=1,IPJM1 SPFAC303 | 309 JH(II+1)=JWORK(II)+JH(II+1) SPFAC304 | 304 CONTINUE SPFAC305 | C SPFAC306 | C END OF INNER LOOP SPFAC307 C SPFAC308 | IF(MP1-J) 310,302,310 SPFAC309 | 310 DO 311 II=1,MP1 SPFAC310 | 311 JWORK(II)=JTEM(II) SPFAC311 | DO 312 II=1,J SPFAC312 312 JTEM(II+1)=-JWORK(II)+JTEM(II+1) SPFAC313 | 302 CONTINUE SPFAC314 | C SPFAC315 | C END OF OUTER LOOP SPFAC316 | C SPFAC317 C SPFAC318 | C DIVIDE BY 2**M SPFAC319 | C SPFAC320 | TEMP=0.5**M SPFAC321 270 | DO 317 I=1,MP1 SPFAC322 317 H(I) = H(I)*TEMP SPFAC323 | C SPFAC324 | C SPFAC325 | C SPFAC326 | C SOLUTION IS STORED IN H IN ASCENDING POWERS OF S. SPFAC327 C SPFAC328 | 1 RETURN SPFAC329 | END SPFAC330 | 1 12/01/64 PAGE 7 | 0| -| $IBFTC MAIN | 1 12/01/64 PAGE 8 | MAIN - EFN SOURCE STATEMENT - IFN(S) - 0| C TEST CASE GENERATOR FOR SUBROUTINE TRY 001 | C SPFAC(G,H,NP,IEND,EPS,IFLAG) TRY 002 | C TRY 003 | DIMENSION G(17),H(9),W(9),JPOW(17) TRY 004 C LOGICAL OUTPUT UNIT (PRINTER) DESIGNATION TRY 005 | LGOUT=6 TRY 006 | NPMAX=16 TRY 007 | NMXP1=NPMAX+1 TRY 008 | DO 9 I=1,NMXP1 TRY 009 9 JPOW(I) = I-1 TRY 010 | ESAVE=0. TRY 011 | ISAVE=1000 TRY 012 | C TRY 013 | C GENERATE A PSUEDO-RANDOM TEST CASE FOR EACH VALUE OF NP. TRY 014 C TRY 015 | ISECT=1 TRY 016 | RV = 3209. TRY 017 | NP=2 TRY 018 | 19 NPP1=NP+1 TRY 019 MP1=NP/2+1 TRY 020 | C TRY 021 | C G(S) = W(S)*W(-S), WHERE W(S) IS AN ARBITRARY POLYNOMIAL WITH TRY 022 | C VALUES ON (-5.0,5.0) TRY 023 | C TRY 024 DO 2 I=1,MP1 TRY 025 | C TRY 026 | C GENERATE RANDOM NUMBER ON (11,9999) TRY 027 | C TRY 028 | RV=RV*RV/100. TRY 029 IV=RV/10000. TRY 030 | FV=IV TRY 031 | IV=RV-FV*10000.+1.0 TRY 032 | IF(IV-10) 21,21,22 TRY 033 | 21 IV=IV+7164 TRY 034 22 RV=IV TRY 035 | C TRY 036 | 2 W(I) = RV*.001-5.0 TRY 037 | C TRY 038 | C FORM G(S) TRY 039 C TRY 040 | DO 5 I=1,NPP1 TRY 041 | 5 G(I) = 0. TRY 042 | DO 3 I=1,MP1 TRY 043 | SIG=-1.0 TRY 044 DO 4 J=1,MP1 TRY 045 | SIG=-SIG TRY 046 | K=I+J-1 TRY 047 | 4 G(K) = G(K) + W(I)*W(J)*SIG TRY 048 | 3 CONTINUE TRY 049 C TRY 050 | C FORCE ODD POWERS OF S TO ZERO TO ACCOUNT FOR POSSIBLE TRY 051 | C ROUND-OFF ERRORS TRY 052 | C TRY 053 | DO 6 I=2,NP,2 TRY 054 6 G(I) = 0. TRY 055 | 1 12/01/64 PAGE 9 | MAIN - EFN SOURCE STATEMENT - IFN(S) - | 0C TRY 056 | 201 WRITE(LGOUT,7) TRY 057 57 7 FORMAT(///8H G(S)) TRY 058 | WRITE(LGOUT,8) (G(I),JPOW(I),I=1,NPP1) TRY 059 58 | 8 FORMAT(3X,E14.7,3HS**,I2,E14.7,3HS**,I2,E14.7,3HS**,I2,E14.7,3HS**TRY 060 | 1,I2) TRY 061 | C TRY 062 EPS=ESAVE TRY 063 | IEND = ISAVE TRY 064 | C ********************************** TRY 065 | CALL SPFAC(G,H,NP,IEND,EPS,IFLAG) TRY 066 | C ********************************** TRY 067 68 IF(IFLAG) 10,11,10 TRY 068 | 10 IEND=-IEND TRY 069 | WRITE(LGOUT,12) IEND TRY 070 72 | 12 FORMAT(30H0 ITERATION DIVERGED AT STEP, I10) TRY 071 | GO TO 1 TRY 072 11 IF(IEND) 10,13,14 TRY 073 | 13 WRITE(LGOUT,15) TRY 074 76 | 15 FORMAT(32H0 ITERATION FAILED TO CONVERGE ) TRY 075 | IEND=ISAVE TRY 076 | GO TO 16 TRY 077 14 WRITE(LGOUT,17) TRY 078 78 | 17 FORMAT(23H0 ITERATION CONVERGED ) TRY 079 | 16 WRITE(LGOUT,18) IEND,EPS TRY 080 79 | 18 FORMAT(16H RESULT AFTER,I10,33H ITERATIONS. ERROR CRITERION WATRY 081 | 1S,E10.3 /8H H(S)) TRY 082 WRITE(LGOUT,8) (H(I),JPOW(I),I=1,MP1) TRY 083 80 | 1 GO TO (20,104),ISECT TRY 084 | C TRY 085 | 20 NP=NP+2 TRY 086 | IF(NP-NPMAX) 19,19,101 TRY 087 C TRY 088 | C TRY 089 | C FACTORIZATION OF BUTTERWORTH FUNCTIONS TRY 090 | C 1 + (-S**2)**M, M=1,2,...,NPMAX/2 TRY 091 | C TRY 092 101 NP=2 TRY 093 | ESAVE=0.00001 TRY 094 | ISECT=2 TRY 095 | SIG=1.0 TRY 096 | 103 NPP1=NP+1 TRY 097 MP1=NP/2+1 TRY 098 | SIG=-SIG TRY 099 | DO 102 I=1,NPP1 TRY 100 | 102 G(I) = 0. TRY 101 | G(1)=1.0 TRY 102 G(NPP1)=SIG TRY 103 | GO TO 201 TRY 104 | 104 NP=NP+2 TRY 105 | IF(NP-NPMAX) 103,103,105 TRY 106 | C TRY 107 C END OF TEST CASE PROGRAM TRY 108 | 105 CALL EXIT TRY 109 | C TRY 110 | C PRINT-OUT WHEN PROGRAM RAN ON 7094 WAS | 1 12/01/64 PAGE 10 MAIN - EFN SOURCE STATEMENT - IFN(S) - | 0C | C | C | C G(S) C 0.4092529E 01S** 0 0.0000000E-38S** 1-0.1314788E 02S** 2 | C | C ITERATION CONVERGED | C RESULT AFTER 7 ITERATIONS. ERROR CRITERION WAS 0.345E-05 | C H(S) C 0.2023000E 01S** 0 0.3626000E 01S** 1 | C | C | C | C G(S) C 0.8482410E 00S** 0 0.0000000E-38S** 1-0.6383223E 01S** 2 0.0000000E-38S** 3 | C 0.5890329E 01S** 4 | C | C ITERATION CONVERGED | C RESULT AFTER 10 ITERATIONS. ERROR CRITERION WAS 0.532E-05 C H(S) | C 0.9210000E 00S** 0 0.3294504E 01S** 1 0.2427000E 01S** 2 | C | C | C C G(S) | C 0.1153282E 02S** 0 0.0000000E-38S** 1-0.2236772E 02S** 2 0.0000000E-38S** 3 | C 0.1053597E 02S** 4 0.0000000E-38S** 5-0.1876901E-01S** 6 | C | C ITERATION CONVERGED C RESULT AFTER 70 ITERATIONS. ERROR CRITERION WAS 0.363E-04 | C H(S) | C 0.3396000E 01S** 0 0.6803302E 01S** 1 0.3521378E 01S** 2 0.1370020E 00S** 3 | C | C C | C G(S) | C 0.2214144E 01S** 0 0.0000000E-38S** 1-0.4938340E 01S** 2 0.0000000E-38S** 3 | C 0.4808792E 02S** 4 0.0000000E-38S** 5 0.2759710E 02S** 6 0.0000000E-38S** 7 | C 0.1896603E 02S** 8 C | C ITERATION CONVERGED | C RESULT AFTER 15 ITERATIONS. ERROR CRITERION WAS 0.154E-03 | C H(S) | C 0.1488000E 01S** 0 0.6425740E 01S** 1 0.1221499E 02S** 2 0.8876680E 01S** 3 C 0.4355000E 01S** 4 | C | C | C | C G(S) C 0.2592100E-01S** 0 0.0000000E-38S** 1-0.2011566E 01S** 2 0.0000000E-38S** 3 | C 0.9321917E 01S** 4 0.0000000E-38S** 5-0.1331592E 02S** 6 0.0000000E-38S** 7 | C 0.2451339E 02S** 8 0.0000000E-38S** 9-0.9672099E-01S**10 | C | C ITERATION CONVERGED C RESULT AFTER 58 ITERATIONS. ERROR CRITERION WAS 0.105E-03 | C H(S) | 1 12/01/64 PAGE 11 | MAIN - EFN SOURCE STATEMENT - IFN(S) - | 0C 0.1610000E 00S** 0 0.2045796E 01S** 1 0.6750671E 01S** 2 0.9292692E 01S** 3 C 0.5503948E 01S** 4 0.3110001E 00S** 5 | C | C | C | C G(S) C 0.8596624E 01S** 0 0.0000000E-38S** 1-0.1415172E 02S** 2 0.0000000E-38S** 3 | C -0.1856197E 02S** 4 0.0000000E-38S** 5-0.3464703E 02S** 6 0.0000000E-38S** 7 | C -0.5289680E 00S** 8 0.0000000E-38S** 9 0.6780199E 00S**10 0.0000000E-38S**11 | C 0.5207524E 01S**12 | C C ITERATION CONVERGED | C RESULT AFTER 14 ITERATIONS. ERROR CRITERION WAS 0.140E-02 | C H(S) | C 0.2932000E 01S** 0 0.1237152E 02S** 1 0.2368738E 02S** 2 0.2844816E 02S** 3 | C 0.2118690E 02S** 4 0.9798928E 01S** 5 0.2282000E 01S** 6 C | C | C | C G(S) | C 0.2231618E 02S** 0 0.0000000E-38S** 1-0.2558518E 02S** 2 0.0000000E-38S** 3 C 0.6405988E 02S** 4 0.0000000E-38S** 5-0.5459582E 02S** 6 0.0000000E-38S** 7 | C 0.6937521E 02S** 8 0.0000000E-38S** 9-0.3831577E 02S**10 0.0000000E-38S**11 | C 0.2768096E 02S**12 0.0000000E-38S**13-0.1029768E 02S**14 | C | C ITERATION CONVERGED C RESULT AFTER 14 ITERATIONS. ERROR CRITERION WAS 0.138E-01 | C H(S) | C 0.4724000E 01S** 0 0.2410740E 02S** 1 0.5880415E 02S** 2 0.8658881E 02S** 3 | C 0.8266239E 02S** 4 0.5095314E 02S** 5 0.1883343E 02S** 6 0.3209000E 01S** 7 | C C | C | C G(S) | C 0.1261129E 01S** 0 0.0000000E-38S** 1-0.1293495E 02S** 2 0.0000000E-38S** 3 | C 0.6943691E 02S** 4 0.0000000E-38S** 5-0.3177253E 02S** 6 0.0000000E-38S** 7 C 0.4119043E 02S** 8 0.0000000E-38S** 9-0.2200648E 02S**10 0.0000000E-38S**11 | C 0.3566076E 01S**12 0.0000000E-38S**13-0.8167420E 01S**14 0.0000000E-38S**15 | C 0.1896129E 01S**16 | C | C ITERATION CONVERGED C RESULT AFTER 16 ITERATIONS. ERROR CRITERION WAS 0.809E-02 | C H(S) | C 0.1123000E 01S** 0 0.9224860E 01S** 1 0.3212961E 02S** 2 0.6109579E 02S** 3 | C 0.7316351E 02S** 4 0.5789598E 02S** 5 0.3012367E 02S** 6 0.9546098E 01S** 7 | C 0.1377000E 01S** 8 C | C | C | C G(S) | C 0.1000000E 01S** 0 0.0000000E-38S** 1-0.1000000E 01S** 2 C | C ITERATION CONVERGED | C RESULT AFTER 1 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) | C 0.1000000E 01S** 0 0.1000000E 01S** 1 1 12/01/64 PAGE 12 | MAIN - EFN SOURCE STATEMENT - IFN(S) - | 0C | C | C C G(S) | C 0.1000000E 01S** 0 0.0000000E-38S** 1 0.0000000E-38S** 2 0.0000000E-38S** 3 | C 0.1000000E 01S** 4 | C | C ITERATION CONVERGED C RESULT AFTER 7 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) | C 0.1000000E 01S** 0 0.1414214E 01S** 1 0.1000000E 01S** 2 | C | C C | C G(S) | C 0.1000000E 01S** 0 0.0000000E-38S** 1 0.0000000E-38S** 2 0.0000000E-38S** 3 | C 0.0000000E-38S** 4 0.0000000E-38S** 5-0.1000000E 01S** 6 | C C ITERATION CONVERGED | C RESULT AFTER 13 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) | C 0.1000000E 01S** 0 0.2000000E 01S** 1 0.2000000E 01S** 2 0.1000000E 01S** 3 | C C | C | C G(S) | C 0.1000000E 01S** 0 0.0000000E-38S** 1 0.0000000E-38S** 2 0.0000000E-38S** 3 | C 0.0000000E-38S** 4 0.0000000E-38S** 5 0.0000000E-38S** 6 0.0000000E-38S** 7 C 0.1000000E 01S** 8 | C | C ITERATION CONVERGED | C RESULT AFTER 19 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) C 0.1000000E 01S** 0 0.2613126E 01S** 1 0.3414214E 01S** 2 0.2613126E 01S** 3 | C 0.1000000E 01S** 4 | C | C | C C G(S) | C 0.1000000E 01S** 0 0.0000000E-38S** 1 0.0000000E-38S** 2 0.0000000E-38S** 3 | C 0.0000000E-38S** 4 0.0000000E-38S** 5 0.0000000E-38S** 6 0.0000000E-38S** 7 | C 0.0000000E-38S** 8 0.0000000E-38S** 9-0.1000000E 01S**10 | C C ITERATION CONVERGED | C RESULT AFTER 27 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) | C 0.1000000E 01S** 0 0.3236068E 01S** 1 0.5236068E 01S** 2 0.5236068E 01S** 3 | C 0.3236068E 01S** 4 0.1000000E 01S** 5 C | C | C | C G(S) | C 0.1000000E 01S** 0 0.0000000E-38S** 1 0.0000000E-38S** 2 0.0000000E-38S** 3 C 0.0000000E-38S** 4 0.0000000E-38S** 5 0.0000000E-38S** 6 0.0000000E-38S** 7 | C 0.0000000E-38S** 8 0.0000000E-38S** 9 0.0000000E-38S**10 0.0000000E-38S**11 | C 0.1000000E 01S**12 | 1 12/01/64 PAGE 13 | MAIN - EFN SOURCE STATEMENT - IFN(S) - 0C | C ITERATION CONVERGED | C RESULT AFTER 33 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) | C 0.1000000E 01S** 0 0.3863703E 01S** 1 0.7464102E 01S** 2 0.9141620E 01S** 3 C 0.7464102E 01S** 4 0.3863703E 01S** 5 0.1000000E 01S** 6 | C | C | C | C G(S) C 0.1000000E 01S** 0 0.0000000E-38S** 1 0.0000000E-38S** 2 0.0000000E-38S** 3 | C 0.0000000E-38S** 4 0.0000000E-38S** 5 0.0000000E-38S** 6 0.0000000E-38S** 7 | C 0.0000000E-38S** 8 0.0000000E-38S** 9 0.0000000E-38S**10 0.0000000E-38S**11 | C 0.0000000E-38S**12 0.0000000E-38S**13-0.1000000E 01S**14 | C C ITERATION CONVERGED | C RESULT AFTER 43 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) | C 0.1000000E 01S** 0 0.4493959E 01S** 1 0.1009783E 02S** 2 0.1459179E 02S** 3 | C 0.1459179E 02S** 4 0.1009783E 02S** 5 0.4493959E 01S** 6 0.1000000E 01S** 7 C | C | C | C G(S) | C 0.1000000E 01S** 0 0.0000000E-38S** 1 0.0000000E-38S** 2 0.0000000E-38S** 3 C 0.0000000E-38S** 4 0.0000000E-38S** 5 0.0000000E-38S** 6 0.0000000E-38S** 7 | C 0.0000000E-38S** 8 0.0000000E-38S** 9 0.0000000E-38S**10 0.0000000E-38S**11 | C 0.0000000E-38S**12 0.0000000E-38S**13 0.0000000E-38S**14 0.0000000E-38S**15 | C 0.1000000E 01S**16 | C C ITERATION CONVERGED | C RESULT AFTER 49 ITERATIONS. ERROR CRITERION WAS 0.100E-04 | C H(S) | C 0.1000000E 01S** 0 0.5125831E 01S** 1 0.1313707E 02S** 2 0.2184615E 02S** 3 | C 0.2568836E 02S** 4 0.2184615E 02S** 5 0.1313707E 02S** 6 0.5125831E 01S** 7 C 0.1000000E 01S** 8 111 | END | 0ERROR MESSAGE NUMBER 1 | 1 12/01/64 PAGE 14 | CF1 STORAGE MAP 0 SUBROUTINE SPFAC | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | WORK 00001 R C 00022 R JH 00033 I | JTEM 00054 I JWORK 00075 I Y 00116 R X 00001 R D 00022 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | NPP1 00216 I M 00217 I MP1 00220 I | K 00221 I J 00222 I JJ 00223 I JM1 00224 I I 00225 I IPJM1 00226 I | SAVE 00227 R FLJ 00230 R II 00231 I | INX 00232 I TEMP 00233 R JEND 00234 I | XINV 00235 R ENRM 00236 R SUM 00237 R | - ENTRY POINTS 0 SPFAC SECTION 2 | - SUBROUTINES CALLED | 0 .XP2. SECTION 3 SQRT 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 | 101 9A 00276 102 69A 00452 103 18A 00325 | 104 52A 00424 106 38A 00403 105 28A 00361 107 33A 00374 108 41A 00407 109 47A 00417 | 110 55A 00432 111 58A 00433 112 64A 00443 | 114 86A 00517 113 74A 00464 115 82A 00510 | 116 90A 00523 2 97A 00532 6 103A 00541 | 5 114A 00621 7 108A 00562 8 111A 00566 201 118A 00625 202 128A 00653 203 135A 00666 | 204 138A 00673 205 190A 01053 206 160A 00740 | 207 140A 00701 209 145A 00712 208 143A 00710 | 212 149A 00717 210 156A 00730 211 159A 00735 | 213 186A 01046 214 184A 01042 215 169A 01006 216 173A 01015 217 171A 01012 218 177A 01025 | 219 183A 01042 220 178A 01027 221 182A 01036 | 222 180A 01032 1 279A 01277 3 193A 01057 | 4 198A 01073 301 207A 01105 302 267A 01256 | 303 216A 01131 304 250A 01232 306 236A 01211 305 226A 01164 307 231A 01202 308 239A 01215 | 1 12/01/64 PAGE 15 | CF1 STORAGE MAP | 0309 245A 01225 310 253A 01240 311 256A 01241 | 312 262A 01251 317 274A 01272 THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 01354.| 1 12/01/64 PAGE 16 | MAIN STORAGE MAP | 0 MAIN PROGRAM | - DIMENSIONED PROGRAM VARIABLES 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | G 00001 R H 00022 R W 00033 R | JPOW 00044 I | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE LGOUT 00065 I NPMAX 00066 I NMXP1 00067 I | I 00070 I ESAVE 00071 R ISAVE 00072 I | ISECT 00073 I RV 00074 R NP 00075 I | NPP1 00076 I MP1 00077 I IV 00100 I | FV 00101 R SIG 00102 R J 00103 I K 00104 I EPS 00105 R IEND 00106 I | IFLAG 00107 I | - ENTRY POINTS | 0 ...... SECTION 2 | - SUBROUTINES CALLED 0 .FVIO. SECTION 3 .FWRD. SECTION 4 SPFAC SECTION 5 | EXIT SECTION 6 .FXEM. SECTION 7 .FFIL. 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 | 9 8A 00232 19 14A 00250 2 25A 00341 | 21 23A 00326 22 24A 00331 5 32A 00353 3 47A 00425 4 41A 00411 6 53A 00433 | 201 57A 00436 7 FORMAT 00134 8 FORMAT 00137 | 10 71A 00515 11 74A 00535 12 FORMAT 00152 | 1 88A 00636 13 76A 00541 14 78A 00557 | 15 FORMAT 00161 16 79A 00572 17 FORMAT 00170 18 FORMAT 00175 20 89A 00645 104 107A 00715 | 101 92A 00654 103 95A 00664 102 101A 00704 | 105 110A 00724 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00755.| 1 12/01/64 PAGE 17 COMPILATION MAIN | 0 DIAGNOSTIC MESSAGES | -1 SOURCE ERROR 20 LEVEL 1 - WARNING ONLY| THE PROGRAM SHOULD END WITH A TRANSFER. A RETURN HAS BEEN GENERATED. 1 12/01/64 PAGE 18 | | 0$ENTRY MAIN 1 IBLDR 12/01/64 PAGE 19 | * MEMORY MAP * | 0 | 0SYSTEM 00000 THRU 02717| FILE BLOCK ORIGIN 02720 FILES 1. UNIT01 | 2. UNIT02 | 3. UNIT03 | 4. UNIT04 | 5. UNIT05 6. UNIT06 | 7. UNIT07 | 8. UNIT08 | FILE LIST ORIGIN 03060| PRE-EXECUTION INITIALIZATION 03100 CALL ON OBJECT PROGRAM 03137| OBJECT PROGRAM 03144 THRU 20476| 0 DECK ORIGIN CONTROL SECTIONS (/NAME/=NON 0 LENGTH, (LOC)=DELETED, 'LOC'=MOVED, *=NOT REFERENCED)| | 1. CF1 03144 SPFAC 04446 2. MAIN 04520 ...... 05457 * | 3. .LXCON 05475 .LXSTR 05475 * .LXSTP 05500 .LXOUT 05546 * .LXERR 05555 .LXCAL 05560 * | .LXRTN 05560 IBEXIT 05560 * .DBCLS 05742 * .LXARG 06111 * .L0 06134 * | .CLSE 06142 .LFBL 06143 * .LUNB 06144 .DFOUT 06145 | 4. .IODEF 06151 .DEFIN 06151 .ATTAC 06155 * .CLOSE 06157 .OPEN 06161 .READ 06163 .WRITE 06165 .BSR 06175 * .READR 06205 .RELES 06207 * .LAREA 06220 | .LFBLK 06236 .LTSX 06241 * .RLHLD 06251 * .AREA1 06254 .LUNBL 06262 | .ENTRY 06266 .GOA 06323 .GO 06327 .DERR 06343 .NOPXI 06344 | .COMXI 06346 .EX34 06370 | 5. .LXSL 06375 .LXSLB 06375 * .LXSEL 06376 .LXSL1 06377 .LXTST 06402 * .LXOVL 06445 * .LXMOD 06507 * .LXIND 06533 * .LXDIS 06536 * .LXFLG 06537 * .LTCH 06540 | 6. .FPTRP 06547 .FFPT. 06547 * .FPOUT 06676 * .FPARG 06704 * /.COUNT/ 06706 * OVFLOW 06752 * | 7. .ERAS. 06757 E.1 06757 E.2 06760 E.3 06761 E.4 06762 | 8. .XCC. 06763 CC.1 06763 CC.2 06764 CC.3 06765 CC.4 06766 | 9. XIT 06767 EXIT 06767 .EXIT. 06767 * 10. FXEM 06770 .FXEM. 06770 TRACE 06776 * /.OPTW./ 07333 * | 11. FCOM 07347 .FCOM. 07347 /FORCOM/ 07452 * | 12. FRDD 07464 .FRDD. 07464 * .FWRD. 07467 | 13. FIOS 07601 .FIOS. 07601 .FSEL. 07770 .FRTB. 10005 .FRTD. 10022 /.FBF1./ 10254 * | 14. FIOH 10260 .FIOH. 10260 .FCNV. 11037 .FFIL. 11134 .FRTN. 11134 * 15. FCNV 11253 .FCON. 11253 | 16. FVIO 14341 .FVIO. 14341 ..FVIO 14341 * | 17. FIOT 14344 .FIOT. 14344 .FIOE. 14435 | 18. UN01 14463 .UN01. 14463 | 19. UN02 14464 .UN02. 14464 20. UN03 14465 .UN03. 14465 | 21. UN04 14466 .UN04. 14466 | 22. UN05 14467 .UN05. 14467 | 23. UN06 14470 .UN06. 14470 ..UN06 14470 * | 24. UN07 14471 .UN07. 14471 25. UN08 14472 .UN08. 14472 | 26. FSQR 14473 SQRT 14473 | 27. FXP2 14546 .XP2. 14546 | 28. .IOCS 14664 .L(0) 14664 .MONSW 14704 .TEOR 14753 .DEFI. 15033 .JOINX 15077 * | .CLOS. 15116 .ATTC. 15131 .SH1 15343 * .SH9 15405 * .OPEN. 15426 .OP4 15454 * .OP7 15505 * .OP9.2 15521 * .RLSE. 15573 .RER2. 15573 | 1 IBLDR 12/01/64 PAGE 20 | * MEMORY MAP * | 0 .READ. 15574 .RER1. 15617 .WRIT. 15621 .MNT1A 16011 * .EOFEX 16072 * | .FEEIT 16142 .GTIOX 16163 .RW7 16301 * .RE7 16720 * .ENDTR 17361 .SEL59 17363 * .BSR. 20002 .EOTOF 20127 .ETOF3 20135 * .SWITC 20164 | .TCHEX 20471 .BASIO 20474 * | 29. .IOCSM 20477 | 0I/O BUFFERS 20477 THRU 77754| 0UNUSED CORE 77755 THRU 77777 0 G(S) 0.4092529E 01S** 0 0. S** 1-0.1314788E 02S** 2 0 ITERATION CONVERGED RESULT AFTER 7 ITERATIONS. ERROR CRITERION WAS 0.345E-05 H(S) 0.2023000E 01S** 0 0.3626000E 01S** 1 G(S) 0.8482410E 00S** 0 0. S** 1-0.6383223E 01S** 2 0. S** 3 0.5890329E 01S** 4 0 ITERATION CONVERGED RESULT AFTER 10 ITERATIONS. ERROR CRITERION WAS 0.532E-05 H(S) 0.9210000E 00S** 0 0.3294504E 01S** 1 0.2427000E 01S** 2 G(S) 0.1153282E 02S** 0 0. S** 1-0.2236772E 02S** 2 0. S** 3 0.1053597E 02S** 4 0. S** 5-0.1876901E-01S** 6 0 ITERATION CONVERGED RESULT AFTER 70 ITERATIONS. ERROR CRITERION WAS 0.363E-04 H(S) 0.3396000E 01S** 0 0.6803302E 01S** 1 0.3521378E 01S** 2 0.1370020E 00S** 3 G(S) 0.2214144E 01S** 0 0. S** 1-0.4938340E 01S** 2 0. S** 3 0.4808792E 02S** 4 0. S** 5 0.2759710E 02S** 6 0. S** 7 0.1896603E 02S** 8 0 ITERATION CONVERGED RESULT AFTER 15 ITERATIONS. ERROR CRITERION WAS 0.154E-03 H(S) 0.1488000E 01S** 0 0.6425740E 01S** 1 0.1221499E 02S** 2 0.8876680E 01S** 3 0.4355000E 01S** 4 G(S) 0.2592100E-01S** 0 0. S** 1-0.2011566E 01S** 2 0. S** 3 0.9321917E 01S** 4 0. S** 5-0.1331592E 02S** 6 0. S** 7 0.2451339E 02S** 8 0. S** 9-0.9672099E-01S**10 0 ITERATION CONVERGED RESULT AFTER 58 ITERATIONS. ERROR CRITERION WAS 0.105E-03 H(S) 0.1610000E 00S** 0 0.2045796E 01S** 1 0.6750671E 01S** 2 0.9292692E 01S** 3 0.5503948E 01S** 4 0.3110001E 00S** 5 G(S) 0.8596624E 01S** 0 0. S** 1-0.1415172E 02S** 2 0. S** 3 -0.1856197E 02S** 4 0. S** 5-0.3464703E 02S** 6 0. S** 7 -0.5289680E 00S** 8 0. S** 9 0.6780199E 00S**10 0. S**11 0.5207524E 01S**12 0 ITERATION CONVERGED RESULT AFTER 14 ITERATIONS. ERROR CRITERION WAS 0.140E-02 H(S) 0.2932000E 01S** 0 0.1237152E 02S** 1 0.2368738E 02S** 2 0.2844816E 02S** 3 0.2118690E 02S** 4 0.9798928E 01S** 5 0.2282000E 01S** 6 G(S) 0.2231618E 02S** 0 0. S** 1-0.2558518E 02S** 2 0. S** 3 0.6405988E 02S** 4 0. S** 5-0.5459582E 02S** 6 0. S** 7 0.6937521E 02S** 8 0. S** 9-0.3831577E 02S**10 0. S**11 0.2768096E 02S**12 0. S**13-0.1029768E 02S**14 0 ITERATION CONVERGED RESULT AFTER 14 ITERATIONS. ERROR CRITERION WAS 0.138E-01 H(S) 0.4724000E 01S** 0 0.2410740E 02S** 1 0.5880415E 02S** 2 0.8658881E 02S** 3 0.8266239E 02S** 4 0.5095314E 02S** 5 0.1883343E 02S** 6 0.3209000E 01S** 7 G(S) 0.1261129E 01S** 0 0. S** 1-0.1293495E 02S** 2 0. S** 3 0.6943691E 02S** 4 0. S** 5-0.3177253E 02S** 6 0. S** 7 0.4119043E 02S** 8 0. S** 9-0.2200648E 02S**10 0. S**11 0.3566076E 01S**12 0. S**13-0.8167420E 01S**14 0. S**15 0.1896129E 01S**16 0 ITERATION CONVERGED RESULT AFTER 16 ITERATIONS. ERROR CRITERION WAS 0.809E-02 H(S) 0.1123000E 01S** 0 0.9224860E 01S** 1 0.3212961E 02S** 2 0.6109579E 02S** 3 0.7316351E 02S** 4 0.5789598E 02S** 5 0.3012367E 02S** 6 0.9546098E 01S** 7 0.1377000E 01S** 8 G(S) 0.1000000E 01S** 0 0. S** 1-0.1000000E 01S** 2 0 ITERATION CONVERGED RESULT AFTER 1 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.1000000E 01S** 1 G(S) 0.1000000E 01S** 0 0. S** 1 0. S** 2 0. S** 3 0.1000000E 01S** 4 0 ITERATION CONVERGED RESULT AFTER 7 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.1414214E 01S** 1 0.1000000E 01S** 2 G(S) 0.1000000E 01S** 0 0. S** 1 0. S** 2 0. S** 3 0. S** 4 0. S** 5-0.1000000E 01S** 6 0 ITERATION CONVERGED RESULT AFTER 13 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.2000000E 01S** 1 0.2000000E 01S** 2 0.1000000E 01S** 3 G(S) 0.1000000E 01S** 0 0. S** 1 0. S** 2 0. S** 3 0. S** 4 0. S** 5 0. S** 6 0. S** 7 0.1000000E 01S** 8 0 ITERATION CONVERGED RESULT AFTER 19 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.2613126E 01S** 1 0.3414214E 01S** 2 0.2613126E 01S** 3 0.1000000E 01S** 4 G(S) 0.1000000E 01S** 0 0. S** 1 0. S** 2 0. S** 3 0. S** 4 0. S** 5 0. S** 6 0. S** 7 0. S** 8 0. S** 9-0.1000000E 01S**10 0 ITERATION CONVERGED RESULT AFTER 27 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.3236068E 01S** 1 0.5236068E 01S** 2 0.5236068E 01S** 3 0.3236068E 01S** 4 0.1000000E 01S** 5 G(S) 0.1000000E 01S** 0 0. S** 1 0. S** 2 0. S** 3 0. S** 4 0. S** 5 0. S** 6 0. S** 7 0. S** 8 0. S** 9 0. S**10 0. S**11 0.1000000E 01S**12 0 ITERATION CONVERGED RESULT AFTER 33 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.3863703E 01S** 1 0.7464102E 01S** 2 0.9141620E 01S** 3 0.7464102E 01S** 4 0.3863703E 01S** 5 0.1000000E 01S** 6 G(S) 0.1000000E 01S** 0 0. S** 1 0. S** 2 0. S** 3 0. S** 4 0. S** 5 0. S** 6 0. S** 7 0. S** 8 0. S** 9 0. S**10 0. S**11 0. S**12 0. S**13-0.1000000E 01S**14 0 ITERATION CONVERGED RESULT AFTER 43 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.4493959E 01S** 1 0.1009783E 02S** 2 0.1459179E 02S** 3 0.1459179E 02S** 4 0.1009783E 02S** 5 0.4493959E 01S** 6 0.1000000E 01S** 7 G(S) 0.1000000E 01S** 0 0. S** 1 0. S** 2 0. S** 3 0. S** 4 0. S** 5 0. S** 6 0. S** 7 0. S** 8 0. S** 9 0. S**10 0. S**11 0. S**12 0. S**13 0. S**14 0. S**15 0.1000000E 01S**16 0 ITERATION CONVERGED RESULT AFTER 49 ITERATIONS. ERROR CRITERION WAS 0.100E-04 H(S) 0.1000000E 01S** 0 0.5125831E 01S** 1 0.1313707E 02S** 2 0.2184615E 02S** 3 0.2568836E 02S** 4 0.2184615E 02S** 5 0.1313707E 02S** 6 0.5125831E 01S** 7 0.1000000E 01S** 8 1 1028 LINES OUTPUT. $IBSYS $STOP | 1 PERIPHERAL UNIT POSITIONS AT END OF JOBS | 1 SYSPP1 IS B2 REC. 00069, FILE 00000 SYSOU1 IS B1 REC. 00364, FILE 00000 SYSIN1 IS A2 REC. 00002, FILE 00001| 1 END OF JOBS | 1