C - SIMEQR-F2-11/13/71-PIKE                    PG 1 OF 3

01.10 G 10.10
C - ZEROING LOOPS, SUBSCRIPT POINTERS
01.40 S T=FP(420+I,0)
01.50 S T=FP(400+I,0)
01.60 F K=1,N;D 1.70;S T=FP(IK,0)
01.70 S IK=I+20*(K-1)
01.80 S IJ=I+20*(J-1)
01.90 S KJ=K+20*(J-1)

C - SUBTRACTION FOR ROW K+1
02.10 D 1.80;D 1.90
02.30 S T=FP(IJ,FP(IJ)-RA*FP(KJ))

C - ELIMINATION BELOW DIAGONAL, COLUMN K
03.20 D 1.70;S KK=K+20*(K-1)
03.30 S RA=FP(IK)/FP(KK)
03.40 F J=KP,N;D 2
03.50 S X=FP(400+I,FP(400+I)-RA*FP(400+K))

C - ROW INTERCHANGE, RIGHT SIDE OF DIAGONAL
04.20 D 1.90;S LJ=L+20*(J-1)
04.30 S T=FP(KJ)
04.40 S X=FP(KJ,FP(LJ))
04.50 S X=FP(LJ,T)

C - SEARCH ROUTINE DOWN COLUMN K
05.20 D 1.70;S LK=L+20*(K-1)
05.30 I (FABS(FP(IK)) - FABS(FP(LK)))5.50,5.50,5.40
05.40 S L=I
05.50 R

C - CONDENSATION, PIVOT ROW K
06.10 S KP=K+1;S L=K
06.15 F I=KP,N;D 5
06.20 I (L-K)6.25,6.45,6.25
06.25 F J=K,N;D 4
06.30 S T=FP(400+K)
06.35 S X=FP(400+K,FP(400+L))
06.40 S X=FP(400+L,T)
06.45 F I=KP,N;D 3

C - GAUSS ELIMINATION ROUTINE
07.10 S NM=N-1;I (NM-1)7.30;
07.20 F K=1,NM;D 6
C - BACK SUBSTITUTION
07.30 S NN=N+20*(N-1)
07.40 S X=FP(420+N,FP(400+N)/FP(NN))
07.50 S I=NM;I (NM-1)7.90;
07.60 S IP=I+1;S SM=0
07.70 F J=IP,N;D 1.80;S SM=SM+FP(IJ)*FP(420+J)
07.75 S II=I+20*(I-1)
07.80 S X=FP(420+I,(FP(400+I)-SM)/FP(II))
07.85 S I=I-1;I (I-1)7.90,7.60,7.60
07.90 R










C - SIMEQR-F2-11/13/71                         PG 2 OF 3

C - MAINLINE, I/O HANDLERS
10.10 T !"SIMEQR-F2"!!
10.12 A "TAPE INPUT? TYPE Y OR N: "X,!
10.14 I (X-0Y)10.16,10.18;
10.16 S X=0;T !"ENTER NO. OF EQS: ";G 10.20
10.18 S X=-1;*
10.20 A N;S NQ=N*N
10.21 I (N)12.30,12.30,10.22
10.22 I (N-20)10.24,10.24,12.30
10.24 F I=1,N;D 1.40;D 1.50;D 1.60
10.26 I (X)10.28;T !!"ENTER MATRIX A:"!
10.28 S I=0
10.30 S I=I+1;D 11
10.32 I (20-RO)10.34;I (I-NQ)10.30;
10.34 I (X)10.36;T !!"ENTER VECTOR B:"!
10.36 S I=0
10.38 S I=I+1;D 12
10.40 I (20-RO)10.41;I (I-N)10.38;
10.41 I (X)10.42,10.43,10.43
10.42 *
10.43 T !!!!"INPUT DATA====="
10.44 D 13
10.46 T %2,!!!!"VECTOR B WITH"N," CONSTANTS"!!
10.48 S I=0
10.50 S I=I+1
10.52 T %2,"     ROW"I,%,"  "FP(400+I),!
10.54 I (I-N)10.50;
10.56 D 7
10.58 T !!!!!"TABULATED SOLUTIONS====="!
10.60 F I=1,N;D 12.20
10.70 T !!!!!;QUIT

C - MATRIX INPUT ROUTINE
11.10 I (X)11.12;T !"ROW: "
11.12 A RO;I (20-RO)11.24;
11.16 I (X)11.18;T "  COL: "
11.18 A CL;I (X)11.20;T "   A? "
11.20 A T;S RC=RO+20*(CL-1)
11.22 S T=FP(RC,T)
11.24 R

C - VECTOR INPUT ROUTINE
12.10 I (X)12.12;T !"ROW: "
12.12 A RO;I (20-RO)12.18;
12.14 I (X)12.16;T "   B? "
12.16 A T;S T=FP(400+RO,T)
12.18 R
12.20 T %2,!"     X("I,")  = "%,FP(420+I),!
12.30 T !!!"IMPROPER NO. OF EQS.  QUIT"
12.40 G 10.70














C - SIMEQR-F2-11/13/71                         PG 3 OF 3

C - PUT ROUTINE FOR SQUARE MATRIX
13.10 S J=1;S ND=4;S LD=20
13.11 S LS=1
13.12 T %2,!!!"SQUARE MATRIX A  - "N," ROWS  AND"N," COLS"!
13.13 S JN=J+ND-1
13.14 I (JN-N)13.16,13.16;S JN=N
13.16 T %2,!"         COLUMN   "
13.18 S JC=J;T JC
13.20 S JC=JC+1;I (JN-JC)13.26;
13.22 T "            "JC;G 13.20
13.26 T !
13.28 S LT=LS+LD-1
13.30 S L=LS-1
13.32 S L=L+1
13.34 T %2,!"     ROW"L
13.36 S K=0
13.38 S K=K+1;S KK=K;S JT=J+K-1
13.40 S LJ=L+20*(JT-1)
13.42 S T(K)=FP(LJ)
13.44 I (JT-N)13.45,13.46,13.46
13.45 I (K-ND)13.38;
13.46 F JW=1,KK;T %,"  "T(JW)
13.48 I (N-L)13.56,13.56;
13.50 I (L-LT)13.32;
13.52 S LS=LS+LD
13.54 G 13.12
13.56 I (JT-N)13.58,13.62,13.62
13.58 S J=JT+1
13.60 G 13.11
13.62 R
*
|