C ----------------------------------------------------
C
C SIMULATION OF FORTRAN PROGRAM
C RUNNING UNDER 1130 DISK MONITOR RELEASE 2 VERSION 12
C ON WINDOWS 10
C
C SOLUTION OF SET OF SIMULTANEOUS LINEAR EQUATIONS
C USING EITHER
C JORDAN METHOD
C OR GAUSS-SIEDEL METHOD
C
C ROGER K. SIMPSON
C PROGRAM WRITTEN NOV 15, 2010
C
C PROGRAM RESURRECTED AND RUN AGAIN OCT 27, 2016
C
C ----------------------------------------------------
C *** TYPE STATEMENTS
REAL CONS,TOLER,TRM
REAL MULT,OLDPV,SUM,COEFF,SOLN
INTEGER NEQU,NROW,NCOL,IROW,COL,IN,OUT,S,T,U,V
INTEGER CWZ,Q
INTEGER TERMF,ORDR,METHD,K
INTEGER KROW,KCOL,SOLCL,ITER
REAL A(20,21)
REAL B(20,21)
REAL SOLTN(20)
REAL SOLNX(20)
INTEGER DDROW(20)
C *** LOGICAL FILE NUMBERS
C 1442=2 1132=3
IN=2
OUT=3
C **********************************
C READ IN DATA
C **********************************
METHD = 1
GO TO (15,18), METHD
15 WRITE (OUT,4)
GO TO 20
18 WRITE (OUT,5)
C -- GET ORDER OF SIMULTANEOUS LINEAR EQUATION SYSTEM
C NEQU IS THE NUMBER OF ROWS IN THE
C INPUT MATRIX.
C NEQU MUST BE IN RANGE OF 16 TO 20
20 READ (IN, 1) NEQU
C -- SETUP NUMBER OF ROWS AND COLUMNS IN MATRIX
NROW = NEQU
NCOL = NROW + 1
C -- POPULATE THE MATRIX
DO 30 IROW = 1, NROW
READ (IN,2) (A(IROW,COL), COL = 1,6)
READ (IN,2) (A(IROW,COL), COL = 7,12)
READ (IN,2) (A(IROW,COL), COL = 13,18)
READ (IN,2) (A(IROW,COL), COL = 19,20)
30 CONTINUE
C ***********************************
C COPY ARRAY A TO ARRAY B
C ***********************************
DO 70 IROW = 1,NROW
DO 75 COL = 1,NCOL
75 B(IROW,COL) = A(IROW,COL)
70 CONTINUE
C ***********************************
C EITHER RUN GAUSS-JORDAN ELIMINATION
C METHOD
C OR GAUSS-JORDAN ITERATIVE METHOD
C METHD = 1 IS JORDAN
C METHD = 2 IS GAUSS-SIEDEL
C ***********************************
GO TO (100,810), METHD
C ***********************************
C PERFORM JORDAN METHOD
C ***********************************
100 DO 800 S=1,NROW
C --- (A) ADJUST 'CURRENT ROW' SO PIVOT VALUE IS 1
DO 200 T=S,COL
IF (T-S) 150,125,150
C SET PIVOT VALUE TO 1
125 OLDPV = A(S,T)
A(S,T) = 1.0
GO TO 200
C NORMALIZE ALL COEFFICIENTS IN CURRENT ROW AS
C REQUIRED SINCE THE PIVOT VALUE WAS SET TO ONE
150 A(S,T) = A(S,T) / OLDPV
200 CONTINUE
C --- (B) ADJUST ALL COEFFICIENTS IN ALL ROWS BESIDES
C THE CURRENT ROW
C 1. THE COEFFICIENT OF THE PIVOT COLUMN IS NOW 0
C 2. MODIFY ALL COEFFICIENTS IN OTHER ROWS AS
C REQUIRED SINCE THE COEFFICIENT IN THE
C PIVOT COLUMN IS NOW ZERO
DO 750 U=1, NROW
IF (U-S) 240, 720, 240
C THIS IS A ROW WE HAVE TO WORK ON
240 CWZ = S
DO 700 V=CWZ,NCOL
IF (V-CWZ) 260,250,260
250 MULT= A(U,V)
A(U,V) = 0.0
GO TO 700
260 A(U,V) = A(U,V) - (MULT * A(S,V))
700 CONTINUE
720 CONTINUE
750 CONTINUE
800 CONTINUE
C ---- MOVE SOLUTION TO 'SOLTN' ARRAY
SOLCL = NEQU+1
DO 805 Q=1,NEQU
805 SOLTN(Q) = A(Q, SOLCL)
GO TO 890
C ***********************************
C GAUSS - SEIDEL ITERATIVE METHOD
C ***********************************
810 ORDR = NEQU
TERMF = 0
ITER = 0
TOLER = 1.0001
C === A. DETERMINE IF MATRIX IS DIAGONALLY DOMINANT
DO 600 KROW=1,ORDR
C --- AT THE OUTSET I ASSUME ROW IS NOT DIAGONALLY
C DOMINANT
DDROW(KROW) = 0
DIAG = A(KROW,KROW)
SUM = 0.0
DO 610 KCOL=1,ORDR
IF (KCOL-KROW) 605,610,605
605 TRM = ABS(A(KROW,KCOL))
SUM = SUM + TRM
610 CONTINUE
C ---- DETERMINE IF THIS ROW MEETS CRITERION FOR
C DIAGONAL DOMINANCE
IF (DIAG-SUM) 600,600,620
C -- THIS ROW IS DIAGONALLY DOMINANT
620 DDROW(KROW) = 1
600 CONTINUE
C ---- IN ORDER FOR THE MATRIX TO BE DIAGONALLY DOMINANT
C THE CRITERION FOR EACH ROW MUST BE MET
C --- INITIALLY I ASSUME MATRIX IS DIAGONALLY DOMINANT
DIAGC = 1
DO 630 KROW=1,ORDR
ROWC = DDROW(KROW)
IF (ROWC) 625,625,630
C --- MATRIX CAN'T BE DIAGONALLY DOMINANT BECAUSE
C THIS ROW DOES NOT MEET CRITERION
625 DIAGC = 0
630 CONTINUE
IF (DIAGC) 634,634,632
632 WRITE (OUT,1005)
GO TO 812
C --- BECAUSE MATRIX IS NOT DIAGONALLY DOMINANT
C I CAN'T PERFORM THE GAUSS-SIEDEL METHOD
634 WRITE (OUT,1006)
GO TO 950
C === B. INITIALIZE SOLUTION MATRIX TO ZERO
812 DO 813 KROW=1,ORDR
813 SOLTN(KROW) = 0.0
C === C. PERFORM ITERATION
C -- BEGIN DO LOOP 1 (OUTER LOOP)
815 DO 850 KROW=1,ORDR
SUM = 0.0
C -- BEGIN DO LOOP 2 (INNER LOOP)
DO 830 KCOL=1,ORDR
IF (KROW-KCOL) 820,830,820
820 COEFF = -A(KROW,KCOL)
SUM = SUM + (COEFF * SOLTN(KCOL))
830 CONTINUE
C -- END DO LOOP 2
K = ORDR+1
CONS = A(KROW, K)
SOLN = (SUM + CONS) / A(KROW, KROW)
SOLNX(KROW) = SOLN
C WRITE(OUT,1002) SOLN
850 CONTINUE
C -- END DO LOOP 1
ITER=ITER+1
C --- DETERMINE IF SOLUTION HAS CONVERGED
TERMF = 1
C -- BEGIN DO LOOP 3
DO 880 KROW=1,ORDR
RATIO = SOLTN(KROW) / SOLNX(KROW)
C WRITE (OUT,1001) RATIO
C --- IF NUMBER 1
IF (RATIO-0.0000001) 852,852,860
852 TERMF = 0
GO TO 880
C --- ELSE NUMBER 1
C -- IF NUMBER 2
860 IF(RATIO-1.0) 865,868,868
C -- IF RATIO .LT. 1.0 THEN
C SET RATIO TO RECIPROCAL OF
C RATIO
865 RATIO = 1.0 / RATIO
C -- ENDIF NUMBER 2
C -- IF NUMBER 3
868 IF(RATIO-TOLER) 880,880,869
C -- IF RATIO IS .GT. TOLERANCE THEN DONT
C TERMINATE
869 TERMF=0
C -- ENDIF NUMBER 3
C -- ENDIF NUMBER 1
880 CONTINUE
C -- END DO LOOP 3
C
C --- PREPARE FOR NEXT ITERATION BY
C COPYING NEXT SOLUTION 'SOLNX' TO
C SOLUTION 'SOLTN'
C --- BEGIN DO LOOP 4
DO 885 KROW=1,ORDR
885 SOLTN(KROW) = SOLNX(KROW)
C --- END DO LOOP 4
IF (TERMF) 815,815,887
C --- PRINT OUT NUMBER OF ITERATIONS
887 WRITE (OUT,1003) ITER
C ***********************************
C DISPLAY SOLUTION
C ***********************************
890 WRITE (OUT,6)
DO 900 Q=1,NEQU
900 WRITE (OUT,8) Q, SOLTN(Q)
C ***********************************
C CHECK THE SOLUTION
C ***********************************
WRITE (OUT,9)
DO 910 KROW=1,NEQU
SUM = 0.0
DO 920 KCOL=1,NEQU
920 SUM = SUM + (SOLTN(KCOL) * B(KROW,KCOL))
910 WRITE (OUT,10) KROW, SUM
950 STOP 16
C ***
C ***
C *** FORMAT STATEMENTS
C ***
C --- INPUT
1 FORMAT(I6)
2 FORMAT(6F11.0)
C --- OUTPUT
4 FORMAT (2X,'JORDAN METHOD -- RK SIMPSON -- OCT 2016')
5 FORMAT (2X,'GAUSS SIEDEL METHOD -- RK SIMPSON -- OCT 2016')
6 FORMAT (/2X,'RESULT -- ')
8 FORMAT (2X,'X( ',I2,' )= ',F10.4)
9 FORMAT (/2X,'CHECK --- ')
10 FORMAT (2X,'C( ',I2,' )= ',F10.4)
1003 FORMAT (2x,'NUMBER OF ITERATIONS = ',I4)
1005 FORMAT (2X,'MATRIX IS DIAGNALLY DOMINANT')
1006 FORMAT (2X,'MATRIX IS NOT DIAGNALLY DOMINANT')
END
// DUP
PAGE 1
// JOB
LOG DRIVE CART SPEC CART AVAIL PHY DRIVE
0000 1234 1234 0000
V2 M12 ACTUAL 16K CONFIG 16K
// DUP
*DELETE SIMEQ
CART ID 1234 DB ADDR 216A DB CNT 0046
// FOR
*IOCS (CARD,1132 PRINTER)
*ONE WORD INTEGERS
*EXTENDED PRECISION
*LIST ALL
C ----------------------------------------------------
C
C SIMULATION OF 1130 DISK MONITOR RELEASE 2 VERSION 12
C ON WINDOWS 10
C
C SOLUTION OF SET OF SIMULTANEOUS LINEAR EQUATIONS
C USING EITHER
C JORDAN METHOD
C OR GAUSS-SIEDEL METHOD
C
C ROGER K. SIMPSON
C PROGRAM WRITTEN NOV 15, 2010
C
C PROGRAM RESURRECTED AND RUN AGAIN OCT 27, 2016
C
C
C ----------------------------------------------------
C *** TYPE STATEMENTS
REAL CONS,TOLER,TRM
REAL MULT,OLDPV,SUM,COEFF,SOLN
INTEGER NEQU,NROW,NCOL,IROW,COL,IN,OUT,S,T,U,V
INTEGER CWZ,Q
INTEGER TERMF,ORDR,METHD,K
INTEGER KROW,KCOL,SOLCL,ITER
REAL A(20,21)
REAL B(20,21)
REAL SOLTN(20)
REAL SOLNX(20)
INTEGER DDROW(20)
C *** LOGICAL FILE NUMBERS
C 1442=2 1132=3
IN=2
OUT=3
C **********************************
C READ IN DATA
C **********************************
METHD = 1
GO TO (15,18), METHD
15 WRITE (OUT,4)
GO TO 20
18 WRITE (OUT,5)
C -- GET ORDER OF SIMULTANEOUS LINEAR EQUATION SYSTEM
PAGE 2
C NEQU IS THE NUMBER OF ROWS IN THE
C INPUT MATRIX.
C NEQU MUST BE IN RANGE OF 16 TO 20
20 READ (IN, 1) NEQU
C -- SETUP NUMBER OF ROWS AND COLUMNS IN MATRIX
NROW = NEQU
NCOL = NROW + 1
C -- POPULATE THE MATRIX
DO 30 IROW = 1, NROW
READ (IN,2) (A(IROW,COL), COL = 1,6)
READ (IN,2) (A(IROW,COL), COL = 7,12)
READ (IN,2) (A(IROW,COL), COL = 13,18)
READ (IN,2) (A(IROW,COL), COL = 19,20)
30 CONTINUE
C ***********************************
C COPY ARRAY A TO ARRAY B
C ***********************************
DO 70 IROW = 1,NROW
DO 75 COL = 1,NCOL
75 B(IROW,COL) = A(IROW,COL)
70 CONTINUE
C ***********************************
C EITHER RUN GAUSS-JORDAN ELIMINATION
C METHOD
C OR GAUSS-JORDAN ITERATIVE METHOD
C METHD = 1 IS JORDAN
C METHD = 2 IS GAUSS-SIEDEL
C ***********************************
GO TO (100,810), METHD
C ***********************************
C PERFORM JORDAN METHOD
C ***********************************
100 DO 800 S=1,NROW
C --- (A) ADJUST 'CURRENT ROW' SO PIVOT VALUE IS 1
DO 200 T=S,COL
IF (T-S) 150,125,150
C SET PIVOT VALUE TO 1
125 OLDPV = A(S,T)
A(S,T) = 1.0
GO TO 200
C NORMALIZE ALL COEFFICIENTS IN CURRENT ROW AS
C REQUIRED SINCE THE PIVOT VALUE WAS SET TO ONE
150 A(S,T) = A(S,T) / OLDPV
200 CONTINUE
C --- (B) ADJUST ALL COEFFICIENTS IN ALL ROWS BESIDES
C THE CURRENT ROW
C 1. THE COEFFICIENT OF THE PIVOT COLUMN IS NOW 0
C 2. MODIFY ALL COEFFICIENTS IN OTHER ROWS AS
C REQUIRED SINCE THE COEFFICIENT IN THE
C PIVOT COLUMN IS NOW ZERO
DO 750 U=1, NROW
IF (U-S) 240, 720, 240
C THIS IS A ROW WE HAVE TO WORK ON
240 CWZ = S
DO 700 V=CWZ,NCOL
IF (V-CWZ) 260,250,260
250 MULT= A(U,V)
A(U,V) = 0.0
GO TO 700
PAGE 3
260 A(U,V) = A(U,V) - (MULT * A(S,V))
700 CONTINUE
720 CONTINUE
750 CONTINUE
800 CONTINUE
C ---- MOVE SOLUTION TO 'SOLTN' ARRAY
SOLCL = NEQU+1
DO 805 Q=1,NEQU
805 SOLTN(Q) = A(Q, SOLCL)
GO TO 890
C ***********************************
C GAUSS - SEIDEL ITERATIVE METHOD
C ***********************************
810 ORDR = NEQU
TERMF = 0
ITER = 0
TOLER = 1.0001
C === A. DETERMINE IF MATRIX IS DIAGONALLY DOMINANT
DO 600 KROW=1,ORDR
C --- AT THE OUTSET I ASSUME ROW IS NOT DIAGONALLY
C DOMINANT
DDROW(KROW) = 0
DIAG = A(KROW,KROW)
SUM = 0.0
DO 610 KCOL=1,ORDR
IF (KCOL-KROW) 605,610,605
605 TRM = ABS(A(KROW,KCOL))
SUM = SUM + TRM
610 CONTINUE
C ---- DETERMINE IF THIS ROW MEETS CRITERION FOR
C DIAGONAL DOMINANCE
IF (DIAG-SUM) 600,600,620
C -- THIS ROW IS DIAGONALLY DOMINANT
620 DDROW(KROW) = 1
600 CONTINUE
C ---- IN ORDER FOR THE MATRIX TO BE DIAGONALLY DOMINANT
C THE CRITERION FOR EACH ROW MUST BE MET
C --- INITIALLY I ASSUME MATRIX IS DIAGONALLY DOMINANT
DIAGC = 1
DO 630 KROW=1,ORDR
ROWC = DDROW(KROW)
IF (ROWC) 625,625,630
C --- MATRIX CAN'T BE DIAGONALLY DOMINANT BECAUSE
C THIS ROW DOES NOT MEET CRITERION
625 DIAGC = 0
630 CONTINUE
IF (DIAGC) 634,634,632
632 WRITE (OUT,1005)
GO TO 812
C --- BECAUSE MATRIX IS NOT DIAGONALLY DOMINANT
C I CAN'T PERFORM THE GAUSS-SIEDEL METHOD
634 WRITE (OUT,1006)
GO TO 950
C === B. INITIALIZE SOLUTION MATRIX TO ZERO
812 DO 813 KROW=1,ORDR
813 SOLTN(KROW) = 0.0
C === C. PERFORM ITERATION
C -- BEGIN DO LOOP 1 (OUTER LOOP)
815 DO 850 KROW=1,ORDR
PAGE 4
SUM = 0.0
C -- BEGIN DO LOOP 2 (INNER LOOP)
DO 830 KCOL=1,ORDR
IF (KROW-KCOL) 820,830,820
820 COEFF = -A(KROW,KCOL)
SUM = SUM + (COEFF * SOLTN(KCOL))
830 CONTINUE
C -- END DO LOOP 2
K = ORDR+1
CONS = A(KROW, K)
SOLN = (SUM + CONS) / A(KROW, KROW)
SOLNX(KROW) = SOLN
C WRITE(OUT,1002) SOLN
850 CONTINUE
C -- END DO LOOP 1
ITER=ITER+1
C --- DETERMINE IF SOLUTION HAS CONVERGED
TERMF = 1
C -- BEGIN DO LOOP 3
DO 880 KROW=1,ORDR
RATIO = SOLTN(KROW) / SOLNX(KROW)
C WRITE (OUT,1001) RATIO
C --- IF NUMBER 1
IF (RATIO-0.0000001) 852,852,860
852 TERMF = 0
GO TO 880
C --- ELSE NUMBER 1
C -- IF NUMBER 2
860 IF(RATIO-1.0) 865,868,868
C -- IF RATIO .LT. 1.0 THEN
C SET RATIO TO RECIPROCAL OF
C RATIO
865 RATIO = 1.0 / RATIO
C -- ENDIF NUMBER 2
C -- IF NUMBER 3
868 IF(RATIO-TOLER) 880,880,869
C -- IF RATIO IS .GT. TOLERANCE THEN DONT
C TERMINATE
869 TERMF=0
C -- ENDIF NUMBER 3
C -- ENDIF NUMBER 1
880 CONTINUE
C -- END DO LOOP 3
C
C
C --- PREPARE FOR NEXT ITERATION BY
C COPYING NEXT SOLUTION 'SOLNX' TO
C SOLUTION 'SOLTN'
C --- BEGIN DO LOOP 4
DO 885 KROW=1,ORDR
885 SOLTN(KROW) = SOLNX(KROW)
C --- END DO LOOP 4
IF (TERMF) 815,815,887
C --- PRINT OUT NUMBER OF ITERATIONS
887 WRITE (OUT,1003) ITER
C ***********************************
C DISPLAY SOLUTION
C ***********************************
890 WRITE (OUT,6)
PAGE 5
DO 900 Q=1,NEQU
900 WRITE (OUT,8) Q, SOLTN(Q)
C ***********************************
C CHECK THE SOLUTION
C ***********************************
WRITE (OUT,9)
DO 910 KROW=1,NEQU
SUM = 0.0
DO 920 KCOL=1,NEQU
920 SUM = SUM + (SOLTN(KCOL) * B(KROW,KCOL))
910 WRITE (OUT,10) KROW, SUM
950 STOP 16
C ***
C ***
C *** FORMAT STATEMENTS
C ***
C --- INPUT
1 FORMAT(I6)
2 FORMAT(6F11.0)
C --- OUTPUT
4 FORMAT (2X,'JORDAN METHOD -- RK SIMPSON -- OCT 2016')
5 FORMAT (2X,'GAUSS SIEDEL METHOD -- RK SIMPSON -- OCT 2016')
6 FORMAT (/2X,'RESULT -- ')
8 FORMAT (2X,'X( ',I2,' )= ',F10.4)
9 FORMAT (/2X,'CHECK --- ')
10 FORMAT (2X,'C( ',I2,' )= ',F10.4)
1003 FORMAT (2X,'NUMBER OF ITERATIONS = ',I4)
1005 FORMAT (2X,'MATRIX IS DIAGNALLY DOMINANT')
1006 FORMAT (2X,'MATRIX IS NOT DIAGNALLY DOMINANT')
END
VARIABLE ALLOCATIONS
CONS(R )=0000 TOLER(R )=0003 TRM(R )=0006 MULT(R )=0009 OLDPV(R )=000C SUM(R )=000F
COEFF(R )=0012 SOLN(R )=0015 A(R )=0501-0018 B(R )=09ED-0504 SOLTN(R )=0A29-09F0 SOLNX(R )=0A65-0A2C
DIAG(R )=0A68 DIAGC(R )=0A6B ROWC(R )=0A6E RATIO(R )=0A71 NEQU(I )=0A74 NROW(I )=0A75
NCOL(I )=0A76 IROW(I )=0A77 COL(I )=0A78 IN(I )=0A79 OUT(I )=0A7A S(I )=0A7B
T(I )=0A7C U(I )=0A7D V(I )=0A7E CWZ(I )=0A7F Q(I )=0A80 TERMF(I )=0A81
ORDR(I )=0A82 METHD(I )=0A83 K(I )=0A84 KROW(I )=0A85 KCOL(I )=0A86 SOLCL(I )=0A87
ITER(I )=0A88 DDROW(I )=0A9C-0A89
STATEMENT ALLOCATIONS
1 =0AB9 2 =0ABB 4 =0ABE 5 =0AD5 6 =0AEF 8 =0AF8 9 =0B02 10 =0B0B 1003 =0B15 1005 =0B25
1006 =0B36 15 =0B72 18 =0B78 20 =0B7C 30 =0BEF 75 =0C00 70 =0C14 100 =0C23 125 =0C31 150 =0C44
200 =0C51 240 =0C64 250 =0C72 260 =0C85 700 =0C9D 720 =0CA6 750 =0CA6 800 =0CAF 805 =0CC2 810 =0CDF
605 =0D15 610 =0D2C 620 =0D3C 600 =0D45 625 =0D66 630 =0D6B 632 =0D79 634 =0D7F 812 =0D85 813 =0D89
815 =0D9B 820 =0DAD 830 =0DC6 850 =0DF8 852 =0E21 860 =0E27 865 =0E2E 868 =0E34 869 =0E3B 880 =0E3F
885 =0E4C 887 =0E62 890 =0E68 900 =0E70 920 =0E96 910 =0EB7 950 =0EC8
FEATURES SUPPORTED
ONE WORD INTEGERS
EXTENDED PRECISION
IOCS
CALLED SUBPROGRAMS
EABS EADD ESUB EMPYX EDIV EDIVX ELD ELDX ESTO ESTOX ESBRX FLOAT CARDZ PRNTZ SRED
SWRT SCOMP SFIO SIOFX SIOF SIOI SUBSC STOP SNR
REAL CONSTANTS
.100000000E 01=0AA0 .000000000E 00=0AA3 .100010000E 01=0AA6 .100000000E-06=0AA9
PAGE 6
INTEGER CONSTANTS
2=0AAC 3=0AAD 1=0AAE 6=0AAF 7=0AB0 12=0AB1 13=0AB2 18=0AB3 19=0AB4 20=0AB5
0=0AB6 16=0AB7 22=0AB8
CORE REQUIREMENTS FOR
COMMON 0 VARIABLES 2720 PROGRAM 1066
END OF COMPILATION
// DUP
*STORE WS UA SIMEQ
CART ID 1234 DB ADDR 21B2 DB CNT 0046
// XEQ SIMEQ L
R 41 24FC (HEX) WDS UNUSED BY CORE LOAD
CALL TRANSFER VECTOR
EABS 17EE
LIBF TRANSFER VECTOR
EBCTB 1A5F
HOLTB 1A23
GETAD 19E0
NORM 19B6
EGETP 19A0
EADDX 180B
XMD 195E
FARC 193C
XDD 18F2
HOLEZ 18BC
PAUSE 18A6
IFIX 187A
STOP 186E
SIOF 11A8
EDIVX 1756
SNR 1866
FLOAT 185C
ESUB 17FA
EADD 1805
ESBRX 17D6
EMPYX 17A8
EDIV 175A
ESTOX 170A
ELDX 1720
SIOFX 11D0
SUBSC 1738
SIOI 11D4
SRED 10CD
SCOMP 11AC
SWRT 10C8
ESTO 170E
ELD 1724
PRNTZ 1622
CARDZ 1572
SFIO 11E9
SYSTEM SUBROUTINES
ILS04 00C4
ILS02 00B3
ILS01 1A64
PAGE 7
ILS00 1A7D
0D47 (HEX) IS THE EXECUTION ADDR
JORDAN METHOD -- RK SIMPSON -- OCT 2016
RESULT --
X( 1 )= -0.0187
X( 2 )= 0.0008
X( 3 )= 0.0078
X( 4 )= -0.0781
X( 5 )= -0.0122
X( 6 )= -0.0508
X( 7 )= 0.0318
X( 8 )= -0.0330
X( 9 )= -0.0605
X( 10 )= -0.0296
X( 11 )= 0.0030
X( 12 )= -0.0656
X( 13 )= -0.0118
X( 14 )= 0.0662
X( 15 )= -0.0263
X( 16 )= -0.0403
X( 17 )= 0.0724
X( 18 )= -0.0182
X( 19 )= 0.0281
CHECK ---
C( 1 )= -27.1962
C( 2 )= 6.7745
C( 3 )= 29.3174
C( 4 )= -80.4052
C( 5 )= -21.9417
C( 6 )= -57.9263
C( 7 )= 51.5459
C( 8 )= -35.4564
C( 9 )= -77.0753
C( 10 )= -32.0387
C( 11 )= 5.9642
C( 12 )= -76.0911
C( 13 )= -44.9409
C( 14 )= 94.5933
C( 15 )= -34.7711
C( 16 )= -67.1532
C( 17 )= 72.1178
C( 18 )= -17.8736
C( 19 )= 51.9958
Updated - October 28, 2016