IMD 1.16: 31/08/2007 11:39:40 fortran programs  LGDSDD un Q@SIMUL FORTMAT PRN TMAT REL TMAT PRNTMAT BAK TMAT PRN!EQUIL0 FOR TMAT COM` !"#$CALCN0 FOR TMAT FOR%&EQUILX PRN')*+,-./EQUILX REL2(1:EEQUILX PRN02345678EQUILX PRN9;<=>?@AEQUILX PRN0BCDEQKON2 $$0EQKON2 BAK qEQKON2 FOR rEQUILX COMuKLMOPQRSEQUILX FOR+FGHEQUILX BAK+IJNEQKON1 PRN6Z]^_BDERV BAKTBDERV FOR USIDERV BAKVSIDERV FOR WEQKON BAKXEQKON FORYEQKON1 REL \EQKON1 COM2; B[A U8b*@7XT:iSqT"b`#u(Ff7iSqt@ZF`Pouҧ€@ kH\sXd#3[@U92~-* i:Z#9 HFf7@W1hZn[HR)~dE[@R)T:%6Z2mkm *=OIǕ-*+m`"`#mk HGIDZU!&kV?X02 ***** 022E' ADD HL,DE ***** 022F' LD A,(HL) ***** 0230' INC HL ***** 0231' LD H,(HL) ***** 0232' LD L,A ***** 0233' LD (IROWK),HL ***** 0236' LD HL,(T:000000) ***** 0239' LD DE,JCOL-0002 ***** 023C' ADD HL,DE ***** 023D' LD A,(HL) ***** 023E' INC HL ***** 023F' LD H,(HL) ***** 0240' LD L,A ***** 0241' LD (JCOLK),HL ***** 0244' LD HL,PIVOT ***** 0247' CALL $L3 ***** 024A' LD HL,DETER ***** 024D' CALL $MU ***** 0250' LD HL,DETER ***** 0253' CALL $T3 36 14 A(IROWK,J) = A(IROWK,J)/PI 02AA' CALL $M9 ***** 02AD' EX DE,HL ***** 02AE' LD HL,(T:000000) ***** 02B1' ADD HL,DE ***** 02B2' LD DE,A-0060 ***** 02B5' ADD HL,DE ***** 02B6' LD (T:010000),HL ***** 02B9' LD HL,[ 00 00 00 81] ***** 02BC' CALL $L1 ***** 02BF' CALL $CK ***** 02C2' LD HL,PIVOT ***** 02C5' CALL $DU ***** 02C8' LD HL,(T:010000) ***** 02CB' CALL $T3 39 AIJCK = A(I,JCOLK) ***** 02CE' LD HL,0001 ***** 02D1' LD (I),HL 40 IF ( I.EQ.IROWK) GO TO 18 ***** 02D4' LD HL,(I) ***** 02D7' ADD HL,HL ****L ***** 013D' LD A,D ***** 013E' SBC H ***** 013F' JP P,00C1' ***** 0142' LD HL,(I) ***** 0145' ADD HL,HL ***** 0146' ADD HL,HL ***** 0147' ADD HL,HL ***** 0148' LD (T:000000),HL ***** 014B' LD HL,(J) ***** 014E' EX DE,HL ***** 014F' LD HL,0058 ***** 0152' CALL $M9 ***** 0155' EX DE,HL ***** 0156' LD HL,(T:000000) ***** 0159' ADD HL,DE ***** 015A' LD DE,A-0060 ***** 015D' ADD HL,DE ***** 015E' LD (T:010000),HL ***** 0161' LD HL,PIVOT ***** 0164' CALL DABS ***** 0167' LD HL,T:000003 ******** 01C0' LD (HL),D ***** 01C1' LD HL,(T:010000) ***** 01C4' LD DE,JCOL-0002 ***** 01C7' ADD HL,DE ***** 01C8' LD (T:030000),HL ***** 01CB' LD HL,(J) ***** 01CE' EX DE,HL ***** 01CF' LD HL,(T:030000) ***** 01D2' LD (HL),E ***** 01D3' INC HL ***** 01D4' LD (HL),D 29 IF (DABS(PIVOT) .GT. EPS) GO TO 13 ***** 01D5' LD HL,(J) ***** 01D8' INC HL ***** 01D9' EX DE,HL ***** 01DA' LD HL,(N) ***** 01DD' LD A,(HL) ***** 01DE' INC HL ***** 01DF' LD H,(HL) ***** 01E0' LD L,A ***** 01E1' EX DfVfPf2>2` v kHR"d ڦkGY 3kXV DZ#"-*V9Btk Xcg)Yb ~TyZpVX#3g8 *ͮ "ͮ#ֱ>2` v! &Ͳʹ,V DZ#'"polf fͤ BE[ R)Y:%6ʀZ2fԀ!ftͪBd#uH='UdU"H\ PBX3khU!Ak tfр&mPmkm@-`*HE[(!, * k4h*V=MZF>2` v*'"ۀpomVE!*kVE"-t* `ͥαV }ePZRAZF`Pouҧʾ@ kHR"d z kGDZU"HG Y$U<UnV?Y*zN<@0Xb ~TyW`afPvLb-Lp #@ *$);# Q Nct3P 2DI^dQH<d16HdqfHsdjvHgNdP6Hdk1&Hcwdʱ6HSd3U`xSIMULHG Y$U<UnV?Y*zN<@0Xb ~TyW`afPvLb-Lp #@ *$) 0281' CALL $DU ***** 0284' LD HL,(T:010000) ***** 0287' CALL $T3 ***** 028A' LD HL,(J) ***** 028D' INC HL ***** 028E' EX DE,HL ***** 028F' LD HL,(MAX) ***** 0292' EX DE,HL ***** 0293' LD A,E ***** 0294' SUB L ***** 0295' LD A,D ***** 0296' SBC H ***** 0297' JP P,0259' 38 DO 18 I = 1,N ***** 029A' LD HL,(IROWK) ***** 029D' ADD HL,HL ***** 029E' ADD HL,HL ***** 029F' ADD HL,HL ***** 02A0' LD (T:000000),HL ***** 02A3' LD HL,(JCOLK) ***** 02A6' EX DE,HL ***** 02A7' LD HL,0058 ********* 011B' LD (T:000002),A 23 8 CONTINUE ***** 011E' OR A ***** 011F' JP NZ,11L 24 9 IF (DABS(A(I,J)).LE.DABS(PIVOT)) GO TO 11 ***** 0122' LD HL,(JSCAN) ***** 0125' INC HL ***** 0126' EX DE,HL ***** 0127' LD HL,(KM1) ***** 012A' EX DE,HL ***** 012B' LD A,E ***** 012C' SUB L ***** 012D' LD A,D ***** 012E' SBC H ***** 012F' JP P,00C7' ***** 0132' LD HL,(ISCAN) ***** 0135' INC HL ***** 0136' EX DE,HL ***** 0137' LD HL,(KM1) ***** 013A' EX DE,HL ***** 013B' LD A,E ***** 013C' SUB ** 0191' LD HL,0058 ***** 0194' CALL $M9 ***** 0197' EX DE,HL ***** 0198' LD HL,(T:000000) ***** 019B' ADD HL,DE ***** 019C' LD DE,A-0060 ***** 019F' ADD HL,DE ***** 01A0' CALL $L3 ***** 01A3' LD HL,PIVOT ***** 01A6' CALL $T3 ***** 01A9' LD HL,(K) ***** 01AC' ADD HL,HL ***** 01AD' LD (T:010000),HL ***** 01B0' LD DE,IROW-0002 ***** 01B3' ADD HL,DE ***** 01B4' LD (T:020000),HL ***** 01B7' LD HL,(I) ***** 01BA' EX DE,HL ***** 01BB' LD HL,(T:020000) ***** 01BE' LD (HL),E ***** 01BF' INC HL  DEC A ***** 020A' SUB 7F ***** 020C' SBC A ***** 020D' LD (T:000002),A 30 SIMUL = 0. ***** 0210' OR A ***** 0211' JP NZ,13L 31 RETURN ***** 0214' LD HL,[ 00 00 00 00] ***** 0217' CALL $L1 ***** 021A' CALL $CK ***** 021D' LD HL,SIMUL ***** 0220' CALL $T3 ***** 0223' RET 32 13 IROWK = IROW(K) 33 JCOLK = JCOL(K) 34 DETER = DETER*PIVOT 35 DO 14 J= 1, MAX ***** 0224' LD HL,(K) ***** 0227' ADD HL,HL ***** 0228' LD (T:000000),HL ***** 022B' LD DE,IROW-00]fcy*o=' :Z[VtdE31EZR`?X\g!VuFrm@jL{ZF`Pouҧ *TyZ`AVu@@h)u*o=' :Z1[ f`Fmm|HFf7@W1h=c S-Vaͮͯ(k "`)kk fcxh G?O31E[HR)pdE[@R)Y:%6:Z2mkmkU@31{J+Z !ZU"H\ BX3op+U!AkVP66 *=OIǕk U@31{J+N !ZBY DZ#VOT ***** 0256' LD HL,0001 ***** 0259' LD (J),HL 37 A(IROWK,JCOLK) = 1./PIVOT ***** 025C' LD HL,(IROWK) ***** 025F' ADD HL,HL ***** 0260' ADD HL,HL ***** 0261' ADD HL,HL ***** 0262' LD (T:000000),HL ***** 0265' LD HL,(J) ***** 0268' EX DE,HL ***** 0269' LD HL,0058 ***** 026C' CALL $M9 ***** 026F' EX DE,HL ***** 0270' LD HL,(T:000000) ***** 0273' ADD HL,DE ***** 0274' LD DE,A-0060 ***** 0277' ADD HL,DE ***** 0278' LD (T:010000),HL ***** 027B' CALL $L3 ***** 027E' LD HL,PIVOT *****  ***** 0531' LD H,(HL) ***** 0532' LD L,A ***** 0533' EX DE,HL ***** 0534' LD A,E ***** 0535' SUB L ***** 0536' LD A,D ***** 0537' SBC H ***** 0538' JP P,04B0' ***** 053B' LD HL,(I) ***** 053E' INC HL ***** 053F' EX DE,HL ***** 0540' LD HL,(NM1) ***** 0543' EX DE,HL ***** 0544' LD A,E ***** 0545' SUB L ***** 0546' LD A,D ***** 0547' SBC H ***** 0548' JP P,04A6' ***** 054B' LD HL,(INTCH) ***** 054E' EX DE,HL ***** 054F' LD HL,0002 ***** 0552' CALL $D9 ***** 0555' ADD HL,HL ***** 0556'6 DO 27 I=1,N ***** 05A4' LD HL,0001 ***** 05A7' LD (J),HL 67 IROWI = IROW(I) ***** 05AA' LD HL,0001 ***** 05AD' LD (I),HL 68 JCOLI = JCOL(I) 69 27 Y(JCOLI) = A(IROWI,J) ***** 05B0' LD HL,(I) ***** 05B3' ADD HL,HL ***** 05B4' LD (T:000000),HL ***** 05B7' LD DE,IROW-0002 ***** 05BA' ADD HL,DE ***** 05BB' LD A,(HL) ***** 05BC' INC HL ***** 05BD' LD H,(HL) ***** 05BE' LD L,A ***** 05BF' LD (IROWI),HL ***** 05C2' LD HL,(T:000000) ***** 05C5' LD DE,JCOL-0002 ***** 05C13' JP P,05AD' 71 28 A(I,J) = Y(I) ***** 0616' LD HL,0001 ***** 0619' LD (I),HL 72 DO 30 I=1,N ***** 061C' LD HL,(I) ***** 061F' ADD HL,HL ***** 0620' ADD HL,HL ***** 0621' ADD HL,HL ***** 0622' LD (T:000000),HL ***** 0625' LD HL,(J) ***** 0628' EX DE,HL ***** 0629' LD HL,0058 ***** 062C' CALL $M9 ***** 062F' EX DE,HL ***** 0630' LD HL,(T:000000) ***** 0633' ADD HL,DE ***** 0634' LD DE,A-0060 ***** 0637' ADD HL,DE ***** 0638' LD (T:010000),HL ***** 063B' LD HL,(T:000000) ******** 0689' ADD HL,DE ***** 068A' LD A,(HL) ***** 068B' INC HL ***** 068C' LD H,(HL) ***** 068D' LD L,A ***** 068E' LD (IROWJ),HL ***** 0691' LD HL,(T:000000) ***** 0694' LD DE,JCOL-0002 ***** 0697' ADD HL,DE ***** 0698' LD A,(HL) ***** 0699' INC HL ***** 069A' LD H,(HL) ***** 069B' LD L,A ***** 069C' LD (JCOLJ),HL 77 DO 30 J=1,N ***** 069F' LD HL,(IROWJ) ***** 06A2' ADD HL,HL ***** 06A3' ADD HL,HL ***** 06A4' ADD HL,HL ***** 06A5' LD DE,Y-0008 ***** 06A8' ADD HL,DE ***** 06A9' LALL $M9 ***** 06FE' EX DE,HL ***** 06FF' LD HL,(T:000000) ***** 0702' ADD HL,DE ***** 0703' LD DE,A-0060 ***** 0706' ADD HL,DE ***** 0707' LD (T:010000),HL ***** 070A' LD HL,(J) ***** 070D' ADD HL,HL ***** 070E' ADD HL,HL ***** 070F' ADD HL,HL ***** 0710' LD DE,Y-0008 ***** 0713' ADD HL,DE ***** 0714' CALL $L3 ***** 0717' LD HL,(T:010000) ***** 071A' CALL $T3 ***** 071D' LD HL,(J) ***** 0720' INC HL ***** 0721' EX DE,HL ***** 0722' LD HL,(N) ***** 0725' LD A,(HL) ***** 0726' INC HL ***095" MAX 00A9" T:000002 00AB" K 00AC" KM1 00AE" I 00B0" J 00B2" ISCAN 00B4" JSCAN 00B6" T:000000 00B8" T:010000 00BA" T:000003 00BC" T:020000 00C4" T:030000 00C6" IROWK 00C8" JCOLK 00CA" AIJCK 00CC" IROWI 00D0" JCOLI 00D2" INTCH 00D4" NM1 00D6" IP1 00D8" JTEMP 00DA" IROWJ 00DC" JCOLJ 00DE" COMMON Length / /0418 (1048) Labels: 5L 0073' 200L 00E0" 18L 03D2' 11L 01D5' 9L 0142' 8L 0122' 13L 0224' 14L 025C' 17L 0356' 20L 0435' 22L 0527' 26L 05A4' **** 057F' LD HL,(INDIC) ***** 0582' LD A,(HL) ***** 0583' INC HL ***** 0584' LD H,(HL) ***** 0585' LD L,A ***** 0586' LD A,L ***** 0587' RLCA ***** 0588' OR L ***** 0589' AND 7F ***** 058B' OR H ***** 058C' DEC A ***** 058D' ADD 81 ***** 058F' SBC A ***** 0590' LD (T:000002),A 63 SIMUL = DETER ***** 0593' OR A ***** 0594' JP NZ,26L 64 RETURN ***** 0597' LD HL,DETER ***** 059A' CALL $L3 ***** 059D' LD HL,SIMUL ***** 05A0' CALL $T3 ***** 05A3' RET 65 26 DO 28 J=1,N 6 HL,0058 ***** 05ED' CALL $M9 ***** 05F0' EX DE,HL ***** 05F1' LD HL,(T:010000) ***** 05F4' ADD HL,DE ***** 05F5' LD DE,A-0060 ***** 05F8' ADD HL,DE ***** 05F9' CALL $L3 ***** 05FC' LD HL,(T:000000) ***** 05FF' CALL $T3 ***** 0602' LD HL,(I) ***** 0605' INC HL ***** 0606' EX DE,HL ***** 0607' LD HL,(N) ***** 060A' LD A,(HL) ***** 060B' INC HL ***** 060C' LD H,(HL) ***** 060D' LD L,A ***** 060E' EX DE,HL ***** 060F' LD A,E ***** 0610' SUB L ***** 0611' LD A,D ***** 0612' SBC H ***** 060667' LD A,(HL) ***** 0668' INC HL ***** 0669' LD H,(HL) ***** 066A' LD L,A ***** 066B' EX DE,HL ***** 066C' LD A,E ***** 066D' SUB L ***** 066E' LD A,D ***** 066F' SBC H ***** 0670' JP P,05A7' 73 DO 29 J=1,N ***** 0673' LD HL,0001 ***** 0676' LD (I),HL 74 IROWJ= IROW(J) ***** 0679' LD HL,0001 ***** 067C' LD (J),HL 75 JCOLJ = JCOL(J) 76 29 Y(IROWJ) = A(I,JCOLJ) ***** 067F' LD HL,(J) ***** 0682' ADD HL,HL ***** 0683' LD (T:000000),HL ***** 0686' LD DE,IROW-0002 HL,(N) ***** 06D9' LD A,(HL) ***** 06DA' INC HL ***** 06DB' LD H,(HL) ***** 06DC' LD L,A ***** 06DD' EX DE,HL ***** 06DE' LD A,E ***** 06DF' SUB L ***** 06E0' LD A,D ***** 06E1' SBC H ***** 06E2' JP P,067C' 78 30 A(I,J) = Y(J) ***** 06E5' LD HL,0001 ***** 06E8' LD (J),HL 79 RETURN ***** 06EB' LD HL,(I) ***** 06EE' ADD HL,HL ***** 06EF' ADD HL,HL ***** 06F0' ADD HL,HL ***** 06F1' LD (T:000000),HL ***** 06F4' LD HL,(J) ***** 06F7' EX DE,HL ***** 06F8' LD HL,0058 ***** 06FB' C** 074B' RET 80 200 FORMAT (10H0N TOO BIG ) 81 END ***** 074C' 00000081 ***** 0750' 00000000 ***** 0754' 0600 Program Unit Length=0756 (1878) Bytes Data Area Length=00EF (239) Bytes Subroutines Referenced: DABS $AT $W2 $ND $L1 $CK $T3 $M9 $SU $L3 $MU $DU $CY $T1 $NU $MR $D9 Variables: SIMUL 0001" N 0005" EPS 0007" INDIC 0009" NRC 000B" A / /+0000 X / /+03C8 DETER 000D" PIVOT 0015" Y 001D" IROW 006D" JCOL 0081" JORD 0 EX DE,HL ***** 0557' LD HL,(INTCH) ***** 055A' LD A,E ***** 055B' SUB L ***** 055C' LD L,A ***** 055D' LD A,D ***** 055E' SBC H ***** 055F' LD H,A ***** 0560' LD A,L ***** 0561' RLCA ***** 0562' OR L ***** 0563' AND 7F ***** 0565' OR H ***** 0566' ADD FF ***** 0568' SBC A ***** 0569' LD (T:000002),A ***** 056C' OR A ***** 056D' JP Z,0000 62 IF (INDIC.LE.0) GO TO 26 ***** 0570' LD HL,DETER ***** 0573' CALL $L3 ***** 0576' CALL $NU ***** 0579' LD HL,DETER ***** 057C' CALL $T3 *8' ADD HL,DE ***** 05C9' LD A,(HL) ***** 05CA' INC HL ***** 05CB' LD H,(HL) ***** 05CC' LD L,A ***** 05CD' LD (JCOLI),HL 70 DO 28 I = 1,N ***** 05D0' LD HL,(JCOLI) ***** 05D3' ADD HL,HL ***** 05D4' ADD HL,HL ***** 05D5' ADD HL,HL ***** 05D6' LD DE,Y-0008 ***** 05D9' ADD HL,DE ***** 05DA' LD (T:000000),HL ***** 05DD' LD HL,(IROWI) ***** 05E0' ADD HL,HL ***** 05E1' ADD HL,HL ***** 05E2' ADD HL,HL ***** 05E3' LD (T:010000),HL ***** 05E6' LD HL,(J) ***** 05E9' EX DE,HL ***** 05EA' LD** 063E' LD DE,Y-0008 ***** 0641' ADD HL,DE ***** 0642' CALL $L3 ***** 0645' LD HL,(T:010000) ***** 0648' CALL $T3 ***** 064B' LD HL,(I) ***** 064E' INC HL ***** 064F' EX DE,HL ***** 0650' LD HL,(N) ***** 0653' LD A,(HL) ***** 0654' INC HL ***** 0655' LD H,(HL) ***** 0656' LD L,A ***** 0657' EX DE,HL ***** 0658' LD A,E ***** 0659' SUB L ***** 065A' LD A,D ***** 065B' SBC H ***** 065C' JP P,0619' ***** 065F' LD HL,(J) ***** 0662' INC HL ***** 0663' EX DE,HL ***** 0664' LD HL,(N) ***** D (T:000000),HL ***** 06AC' LD HL,(I) ***** 06AF' ADD HL,HL ***** 06B0' ADD HL,HL ***** 06B1' ADD HL,HL ***** 06B2' LD (T:010000),HL ***** 06B5' LD HL,(JCOLJ) ***** 06B8' EX DE,HL ***** 06B9' LD HL,0058 ***** 06BC' CALL $M9 ***** 06BF' EX DE,HL ***** 06C0' LD HL,(T:010000) ***** 06C3' ADD HL,DE ***** 06C4' LD DE,A-0060 ***** 06C7' ADD HL,DE ***** 06C8' CALL $L3 ***** 06CB' LD HL,(T:000000) ***** 06CE' CALL $T3 ***** 06D1' LD HL,(J) ***** 06D4' INC HL ***** 06D5' EX DE,HL ***** 06D6' LD ** 0727' LD H,(HL) ***** 0728' LD L,A ***** 0729' EX DE,HL ***** 072A' LD A,E ***** 072B' SUB L ***** 072C' LD A,D ***** 072D' SBC H ***** 072E' JP P,06E8' ***** 0731' LD HL,(I) ***** 0734' INC HL ***** 0735' EX DE,HL ***** 0736' LD HL,(N) ***** 0739' LD A,(HL) ***** 073A' INC HL ***** 073B' LD H,(HL) ***** 073C' LD L,A ***** 073D' EX DE,HL ***** 073E' LD A,E ***** 073F' SUB L ***** 0740' LD A,D ***** 0741' SBC H ***** 0742' JP P,0676' ***** 0745' LD HL,SIMUL ***** 0748' CALL $L3 *** f*3#*,{z!!*.}20b! !!!"3I! !"5*3)))"7*5!X *7"9*9! >f*5#*,{zG !!*3#*,{z8 b.(10H0N TOO BIG)" " >! )(* ~#fo"I * ~#fo}?2K * ~#fo#"I * ~#fo}=Ɓ2K   !!!!h'}! ͙*X ͟*P #* ~#fo{z!"t * ~#fo+"v !"P *P #"x "R *R )3 "X *P )3 "Z *X ~#fo*Z ~#fo{ozg}?2K ·*R )3 "X ~#fo"z *P )3 "Z ~#fo*X s#r*z *Z s#r*t #"t *R #* ~#fo{z@*P #*v {z6*t !)*t {ozg}2K ! ͙H! ͟* ~#fo}=Ɓ2K 4! ͙! ͟!"R !"P *P )"X  ~#fo"p *X  ~#fo"r *r ))) "X *p )))"Z *R !X *Z ͙*X ͟*P #* ~#fo{z=!"P *P )))"X *R !X *X "Z *X  ͙*Z ͟*P #* ~#fo{z*R #* ~#fo{zt!0ͦh'} ͦtl!0y20!0/~q+6qWNsY+LB W~w+aZ!0^͡}h'!0w#ƒ͒}͙͛ð͙(!0F#^#V#N&Ͱ}!0h'>T.͑''!065g>(|>¿akh'H:HØ'͉'V͘'^=&=t|t'2!HC=e>O!0~+>wØ'QYq'Ø':0Ҝhڔ=:Ɓʜ͵!"JQYq'͘' '&|Yt&wz^Pc|u~r1q'T.>20!"=og͡S }}{(qn$!~00w# :0@*/DM*/"0*/ */_`i"/"/*/DM!"0*/*/ _`i"/<:}0*0+}|ҧ:0*0"/>; */_>; */"0*/_w#=b*/+"/|_<>G!/60#&!/4:00w&"{0~: w+w!/|}!/4#"/|/W}/_*{0"{0>2/:}0+%:/G:0Ix2/:0X2/[!:0y0ͩ!!/5ʊþ !/5G!:y0*!!{050!*/F#"/é!!i0~!0!/2y0y!!2//!:0-ĩ!!/5!.ͩ!6!Â!T. ͩ!!y05š!y!*/*/{z!#"/*/+p!w05!T.:x0:0$"!0~4!T. *{0"0*0"0$!"w0*u0*m0:02"~'#"w#"#~'""u0$~'I"Gͩ!*u0#"u0#2"*u0#"u0#~'8"$:0=p"!04!"u0"20$*0+"0|$*0"u0!04$!/6:0¸":k0¸"*/z±"*/6 "/%:/2/2/:0"!"/\#!"/$!0~"w"*/+"/|":0"* 4&go"0"0:.<2.!/6P#6#"f/6 :D'=G!(w#6'2(2- -----((((('s#r#q#pÑ'͉'t'͑'"0`i"0*0*0DM^#V#N#F#xyà':0:07#!#!#¹'###!/4|/G}/O!>' '7>{_zW}o|g='|g}o(CZQ( o-yOzW{_xG(`iN#F#q#p#=,(2/"0"0`i"/2/:/d(!0_x(:/*/}(_~#fo:/<2/W:/Ґ(W2 27 IF (INDIC.NE.0) GO TO 1 ***** 0288' LD HL,0001 *****0FORTDATÚ+3+ +**(*!(:0_~!(:*X "Z ͙! *Z ͟*R #*I {z *h )))"X *j !X *X "Z !h'}! *Z ͟!"P *P )))"X *j !X *X ͙͛!l Y'*P *h {ozg}2K b*P )))"X *j !X *X "Z !l h'}! H*Z ͟!"R *R *j {ozg}2K R*P )))"X *R !X "Z *X "d *h )))"f *Z *f ͙!l m*d H*d ͟*R #*I {z *P #* ~#fo{za *L #* ~#fo{z !"P *P )"X  ~#fo"p *X  ~#fo"r *p )3 "Z *r *Z s#r* ~#fo}?2K *r )))"X *p )))"Z *I !X *Z  ':0 '!0N#F)q<Gnlx %ç!0H4j0 TeB׳]h!I.k ͦh'}Ãͦt!0T.!0q#~++w+qžØ'DNn"~`35zr1{r1h!I͂ʘ'yO20͂ʘ':0l/))# =!0~7w?##wy7O*0*0"0"0 >i>i>i>2}0;(O(~#fo"g0O("/*g0"/:0ʔz"zT.>2020<2k0>2020<>20202k0~#fo"0!/"/~#+!"m0!9##"R:0 go"o0"q0"s08!o0>)(*o0~# go"o0*q0~#.go"q0*s0~#fo"s0202020202x02/go"u0"/'*m0|pj#*u0"0("T.2/2/2/2/<20:0>þ & š> $wÍ://<ʳ$=¬2/:}0%/ %$=&:/ ͩ!=2/:}0>$=:~0Gͩ!&è':}0==T.:~0og"~0go"0"020>20+%:/:/2/2/[!_ V V>00g:g020>ɯ2/2/2j0_ -…>20Ê+_ %*/+"/!w04& :0.µ>2/2/& :0ED_ +->2/_ y2j0*/+"/!w04F Ox  T.> G:/x /<2/G!/~wL:0& !&&&:0!0w#wK !/:/w&  _ :i02w0_ Fʁ T>2~0 %!~0~ʙ >T>Fw> #w#w#w>2}0> ø > ø >O}20O!D':0Қʚ=_^#V ^#V *q0|!T.*o0|!T.*:}0G*/~0%0w#%"/*0|$"!0~ʵ%x2/w:}0=T%=T%+>20>wk%!0~%/W+~/_s#r>207Ԩ'K&!/~wl%%&o%¨%!0~0 !&o%!~00 !!/4/G#+k%2/!~0w#½%!/4>þ *0*/"/yx$T. */"/$ ͩ!*0+"0|%$:}0G:i02/!/>w!/5#!#ajSXu& 'JS\E!0 '/F+N+=!0!0>q#p#=!0"/ͬ*/z ='"=~#h'=x͉'S-ͳP͉'ʘ':0q'c/<Ҙ'(g(!0~_#~W#~O&.(&G~_#~W#~O.&ͳò͉'ʘ'>OSͳ͉'t'͘'y+F+F+Fw`h|gy:0O|g}oxG-|}.&ElaO!0~Gx2ƀ 'w(w+ɷ& ''x{t'!065wgXzW >(|&{>20.&/2w0 ͩ!È*u0*m0DM*0+|"0"u0$*u0*m0Fͩ!*u0#"u0*0+"0|#$*/DM!w05*/+}| */ ~20`i"/>>  !w04:l0Ij:}0=2=9:/=T.!0NT.S#C:0%*~0|/g}/o#"~0%>2/:}0{T.!/:/:j0:/w~/g"~0|x!0~w+ ##!0~w+#!0##0!~0w# $6 j#"05$,$/ʏ")\"'!PC$:~020$|S$T.!"02l0!$l$# ]$T.yO!$:0|$!$ ^#V"0yҙ$*u0"{0!"w0j#£$T.}2i02/y $:0.$j#*~0}2/*u0+"u0"T.XH(AILEFGD%!kh kkkk%! -; QTURS#EU%0A`P`K@9- B8lb`" 1k3@! .ͣ@CHG6i`#J@4H9CP&m:f! ڌͪCVBx7 E" 8:5|f#uH='UX5C`tBx7 E" 8:6|f؀VoY{J+Rͳ`E :sX3@"<! H!YU"H8 BX3@X8B0M <:7|f0Vw Q8 ۬*1*>'W\%@2 -JQh9@hYn o6ŴBHdq0HdjHgd*])PSӠT44!`XK@Z!X$BX E" `XfV `A`3@V  Q8#ңqB 63i3ifmbmhmnffp!`V E!2`ȋ hCc&h* 1*>'VPK@&1 P pnEALB0:X2"ͯ * } 6EXDX<B0:X2 4m`mm&mV  udE"U@",J!bbpV%fV'f@)fѠV#f8***** 00F6' LD HL,[ 01 00] ***** 00F9' CALL $W2 21 2 CONTINUE ***** 00FC' LD HL,0001 ***** 00FF' LD (J),HL ***** 0102' LD HL,(I) ***** 0105' ADD HL,HL ***** 0106' ADD HL,HL ***** 0107' ADD HL,HL ***** 0108' LD (T:000000),HL ***** 010B' LD HL,(J) ***** 010E' EX DE,HL ***** 010F' LD HL,0058 ***** 0112' CALL $M9 ***** 0115' EX DE,HL ***** 0116' LD HL,(T:000000) ***** 0119' ADD HL,DE ***** 011A' LD DE,A-0060 ***** 011D' ADD HL,DE ***** 011E' LD (T:010000),HL ***** 0121' LD HL,(T:010000)  ***** 0181' LD (ITCON),HL 28 IF (DABS(XINC(I)).GT.EPS2) ITCON = 0 ***** 0184' LD HL,0001 ***** 0187' LD (I),HL ***** 018A' LD HL,(I) ***** 018D' ADD HL,HL ***** 018E' ADD HL,HL ***** 018F' ADD HL,HL ***** 0190' LD DE,XINC-0008 ***** 0193' ADD HL,DE ***** 0194' LD (T:000000),HL ***** 0197' LD HL,(T:000000) ***** 019A' CALL DABS ***** 019D' LD HL,EPS2 ***** 01A0' CALL $SU ***** 01A3' DEC A ***** 01A4' SUB 7F ***** 01A6' SBC A ***** 01A7' LD (T:000002),A ***** 01AA' OR A ***** 01AB**** 0204' LD A,02 ***** 0206' CALL $I3 ***** 0209' LD DE,N ***** 020C' LD HL,[ 01 00] ***** 020F' LD A,02 ***** 0211' CALL $I0 ***** 0214' CALL $ND ***** 0217' LD DE,306L ***** 021A' LD HL,[ 01 00] ***** 021D' CALL $W2 32 IF (ITCON.EQ.0) GO TO 9 ***** 0220' LD HL,0001 ***** 0223' LD (I),HL ***** 0226' LD HL,(I) ***** 0229' ADD HL,HL ***** 022A' ADD HL,HL ***** 022B' ADD HL,HL ***** 022C' LD DE,XOLD-0008 ***** 022F' ADD HL,DE ***** 0230' LD (T:000000),HL ***** 0233' LD HL,(T:00000HL ***** 0292' ADD HL,HL ***** 0293' ADD HL,HL ***** 0294' LD DE,XOLD-0008 ***** 0297' ADD HL,DE ***** 0298' LD (T:000000),HL ***** 029B' LD HL,(T:000000) ***** 029E' EX DE,HL ***** 029F' LD HL,[ 01 00] ***** 02A2' LD A,02 ***** 02A4' CALL $I3 ***** 02A7' LD HL,(I) ***** 02AA' INC HL ***** 02AB' EX DE,HL ***** 02AC' LD HL,(N) ***** 02AF' EX DE,HL ***** 02B0' LD A,E ***** 02B1' SUB L ***** 02B2' LD A,D ***** 02B3' SBC H ***** 02B4' JP P,028B' ***** 02B7' CALL $ND 36 9 CONTINUE ****JBz G "JR@= j/ @ O&H\.,ɄA@2 ILp ]1fK@! `S*dRLU&dqXR#I2(^1 a$f/ @ O&H\.,ɄA@2 ILp \1KE @N'dqXR#I2([.dKẁ"`X'!I LJ%RI0B)J L!I$NeRp ]( 9@E ȢI8U* RiRIPR $qLN)Kf*)ub0, aS" $T )%"$G$©'w3)n`P8"Fᠦ[(ҁ2D8$I0̀D%9)*@&HedȑHeh)iё0HsTi)p ***** 00CD' EX DE,HL ***** 00CE' LD A,E ***** 00CF' SUB L ***** 00D0' LD A,D ***** 00D1' SBC H ***** 00D2' JP P,00A9' ***** 00D5' CALL $ND 17 DO 9 ITER = 1,ITMAX ***** 00D8' LD HL,0001 ***** 00DB' LD (INDIC),HL 18 CALL CALCN(XOLD,NRC ) ***** 00DE' LD HL,0001 ***** 00E1' LD (ITER),HL 19 DO 2 I=1,7 ***** 00E4' LD DE,NRC ***** 00E7' LD HL,XOLD ***** 00EA' CALL CALCN 20 WRITE (1,306) (A(I,J),J=1,8) ***** 00ED' LD HL,0001 ***** 00F0' LD (I),HL ***** 00F3' LD DE,306L **** 014D' LD DE,EPS1 ***** 0150' LD HL,N ***** 0153' CALL SIMUL ***** 0156' CALL $CK ***** 0159' LD HL,DETER ***** 015C' CALL $T3 ***** 015F' LD HL,0000 ***** 0162' CALL $SQ ***** 0165' ADD FF ***** 0167' SBC A ***** 0168' LD (T:000002),A 24 WRITE (1,201) ***** 016B' OR A ***** 016C' JP NZ,3L ***** 016F' LD DE,201L ***** 0172' LD HL,[ 01 00] ***** 0175' CALL $W2 25 GO TO 1 ***** 0178' CALL $ND 26 3 ITCON = 1 ***** 017B' JP 1L 27 DO 5 I=1,N ***** 017E' LD HL,0001 HL,(T:010000) ***** 01D7' CALL $T3 ***** 01DA' LD HL,(I) ***** 01DD' INC HL ***** 01DE' EX DE,HL ***** 01DF' LD HL,(N) ***** 01E2' EX DE,HL ***** 01E3' LD A,E ***** 01E4' SUB L ***** 01E5' LD A,D ***** 01E6' SBC H ***** 01E7' JP P,0187' ***** 01EA' LD DE,202L ***** 01ED' LD HL,[ 01 00] ***** 01F0' CALL $W2 31 WRITE(1,306) (XOLD(I),I=1,N) ***** 01F3' LD DE,ITER ***** 01F6' LD HL,[ 01 00] ***** 01F9' LD A,02 ***** 01FB' CALL $I0 ***** 01FE' LD DE,DETER ***** 0201' LD HL,[ 01 00] *002),A 33 WRITE (1,203) ITER,N ***** 0261' OR A ***** 0262' JP NZ,9L ***** 0265' LD DE,203L ***** 0268' LD HL,[ 01 00] ***** 026B' CALL $W2 34 WRITE (1,306) (XOLD(I),I=1,N) ***** 026E' LD BC,N ***** 0271' LD DE,ITER ***** 0274' LD HL,[ 01 00] ***** 0277' LD A,03 ***** 0279' CALL $I0 ***** 027C' CALL $ND ***** 027F' LD DE,306L ***** 0282' LD HL,[ 01 00] ***** 0285' CALL $W2 35 GO TO 1 ***** 0288' LD HL,0001 ***** 028B' LD (I),HL ***** 028E' LD HL,(I) ***** 0291' ADD HL,fVoE!gVqfc|eBưU"H8wdE@Y#,U@<3j7b=OIǕͮ!!#6sLV >3o 6o3ofVoE!gVquh͠B#uH='T E@ZZF>2 v4 G9#@5|f &mPр! uf!")B;2" 0! O@YFfY*zN u!U@, ! C@0C@CB :X2ͦQYC0B:X2ͪND 6p(6p(66y 6BXYB3lH+6|@ 6ʂDXDZB:X2"ՠ* m ! L ! ۸ *  ! ^ *ՠ ۠ *` !` ! !L !bf@Q0fa``6@mfڐQfaVfҀafpa `6iX363i0j00udE U@"-r!bfaVfߠQfaVfaVf8' CALL $L1 ***** 02FB' CALL $CK ***** 02FE' LD HL,A+0018 ***** 0301' CALL $T3 ***** 0304' LD HL,[ 40 C0 87 92] ***** 0307' CALL $L1 ***** 030A' CALL $CK ***** 030D' LD HL,A+0070 ***** 0310' CALL $T3 ***** 0313' LD HL,[ 80 C2 98 91] ***** 0316' CALL $L1 ***** 0319' CALL $CK ***** 031C' LD HL,A+00C8 ***** 031F' CALL $T3 ***** 0322' LD HL,[ 00 DE 13 8F] ***** 0325' CALL $L1 ***** 0328' CALL $CK ***** 032B' LD HL,A+0120 ***** 032E' CALL $T3 ***** 0331' LD HL,[ 00 AC 03 8E] ***** 0334' CALL ***** 03A6' EX DE,HL ***** 03A7' LD HL,(XOLD) ***** 03AA' ADD HL,DE ***** 03AB' LD (T:0F0000),HL ***** 03AE' LD HL,0010 ***** 03B1' EX DE,HL ***** 03B2' LD HL,(XOLD) ***** 03B5' ADD HL,DE ***** 03B6' LD (T:100000),HL ***** 03B9' LD HL,0018 ***** 03BC' EX DE,HL ***** 03BD' LD HL,(XOLD) ***** 03C0' ADD HL,DE ***** 03C1' LD (T:110000),HL ***** 03C4' LD HL,0020 ***** 03C7' EX DE,HL ***** 03C8' LD HL,(XOLD) ***** 03CB' ADD HL,DE ***** 03CC' LD (T:120000),HL ***** 03CF' LD HL,0030 ***** 03D2' 0442' LD HL,(T:0F0000) ***** 0445' CALL $L3 ***** 0448' LD HL,[ 40 C0 07 92] ***** 044B' CALL $MR ***** 044E' LD HL,T:0C0003 ***** 0451' CALL $T3 ***** 0454' LD HL,[ 00 4A 61 8F] ***** 0457' CALL $L1 ***** 045A' CALL $CK ***** 045D' LD HL,(XOLD) ***** 0460' CALL $MU ***** 0463' LD HL,T:0C0003 ***** 0466' CALL $AU ***** 0469' LD HL,T:0B0003 ***** 046C' CALL $AU ***** 046F' LD HL,T:0A0003 ***** 0472' CALL $SU ***** 0475' LD HL,T:090003 ***** 0478' CALL $SU ***** 047B' LD HL,T:080003 ***** **** 0251' CALL $T3 ***** 0254' LD HL,[ 00 00 00 81] ***** 0257' CALL $L1 ***** 025A' CALL $CK ***** 025D' LD HL,A+0010 ***** 0260' CALL $T3 ***** 0263' LD HL,[ 00 00 00 81] ***** 0266' CALL $L1 ***** 0269' CALL $CK ***** 026C' LD HL,A+0068 ***** 026F' CALL $T3 ***** 0272' LD HL,[ 00 00 00 81] ***** 0275' CALL $L1 ***** 0278' CALL $CK ***** 027B' LD HL,A+0170 ***** 027E' CALL $T3 ***** 0281' LD HL,0030 ***** 0284' EX DE,HL ***** 0285' LD HL,(XOLD) ***** 0288' ADD HL,DE ***** 0289' CALL $bEٴh CY B3i06{ 6FXDY|B3i663j`iPudE@#U@,!@C U@ CSٴB3@cfՀ1(f01 `ȋ5BXDYB:X2"* qV Pueb@60 cYf9h@6H SٵH U C@ U@ CLٷ C@(U@CHٴXC@pU@B3i@ J3ip 53i fӠAfAͿ`l* C@$C@C@C@C@C@CHCs`V<b68dAPChCgAh?@jA6kCYCnpA>q T:0B0000),HL ***** 02C3' LD HL,0030 ***** 02C6' EX DE,HL ***** 02C7' LD HL,(XOLD) ***** 02CA' ADD HL,DE ***** 02CB' LD (T:0C0000),HL ***** 02CE' LD HL,[ 00 00 00 81] ***** 02D1' CALL $L1 ***** 02D4' CALL $CK ***** 02D7' LD HL,(T:0C0000) ***** 02DA' CALL $DU ***** 02DD' LD HL,(T:0B0000) ***** 02E0' CALL $SU ***** 02E3' LD HL,(T:0A0000) ***** 02E6' CALL $SU ***** 02E9' LD HL,(XOLD) ***** 02EC' CALL $SU ***** 02EF' LD HL,A+0278 ***** 02F2' CALL $T3 ***** 02F5' LD HL,[ 00 4A E1 8F] ***** 02F:0E0000),HL ***** 036E' LD HL,0030 ***** 0371' EX DE,HL ***** 0372' LD HL,(XOLD) ***** 0375' ADD HL,DE ***** 0376' CALL $L3 ***** 0379' LD HL,0002 ***** 037C' CALL $EQ ***** 037F' LD HL,T:060003 ***** 0382' CALL $T3 ***** 0385' LD HL,(T:0E0000) ***** 0388' CALL $L3 ***** 038B' LD HL,[ 00 08 27 8E] ***** 038E' CALL $MR ***** 0391' LD HL,[ 00 D0 52 8E] ***** 0394' CALL $SR ***** 0397' LD HL,T:060003 ***** 039A' CALL $DU ***** 039D' LD HL,A+0228 ***** 03A0' CALL $T3 ***** 03A3' LD HL,0008 T:080003 ***** 0409' CALL $T3 ***** 040C' LD HL,(T:120000) ***** 040F' CALL $L3 ***** 0412' LD HL,[ 00 AC 03 8E] ***** 0415' CALL $MR ***** 0418' LD HL,T:090003 ***** 041B' CALL $T3 ***** 041E' LD HL,(T:110000) ***** 0421' CALL $L3 ***** 0424' LD HL,[ 00 DE 13 8F] ***** 0427' CALL $MR ***** 042A' LD HL,T:0A0003 ***** 042D' CALL $T3 ***** 0430' LD HL,(T:100000) ***** 0433' CALL $L3 ***** 0436' LD HL,[ 80 C2 18 91] ***** 0439' CALL $MR ***** 043C' LD HL,T:0B0003 ***** 043F' CALL $T3 ***** 5' LD (T:090000),HL ***** 0218' LD HL,[ 00 00 00 82] ***** 021B' CALL $L1 ***** 021E' CALL $CK ***** 0221' LD HL,(T:090000) ***** 0224' CALL $DU ***** 0227' LD HL,T:040003 ***** 022A' CALL $T3 ***** 022D' LD HL,(T:080000) ***** 0230' CALL $L3 ***** 0233' LD HL,[ 00 00 00 82] ***** 0236' CALL $MR ***** 0239' LD HL,T:040003 ***** 023C' CALL $SU ***** 023F' LD HL,(T:070000) ***** 0242' CALL $AU ***** 0245' LD HL,(T:060000) ***** 0248' CALL $AU ***** 024B' CALL $NU ***** 024E' LD HL,A+0270 *VEXWD5YD8qH4bm|fhf֠ f@mm@mf٠f sf`mmV kf 9f fPmmR&mV 0udE"U@",!bdp fݰmdpmmdPmѐ&h! *` * چͱ`"ͥ`BFͣL>! ڤ!Lڰͫ`Bͬ BͬL!,!bf!͹"ͯBͮL! J!L,!be`ȋ-FXDYtCf@p gYti@YlO@Ydx@Xy= j 63n3ofߐ&mm ͡`b͢ @a Qf@af`afؠfL 3fL#13 p$d1FHcdiqfHg1zx$W22DIALUA >c%8QBt $f"ha#1%7h XGLy::chV\0$3 (dpfHdq`VHvdQVHdjfH djQ`fH5T44TSUS U4TȠ!`PXBX#3ZU@31Z@>2` vP*#j pfcxE @W1h=c S-VaB3@h!4hɘ@mgf֠ f"ՀX+j dfpmzPm@-`!Z$L1 ***** 0337' CALL $CK ***** 033A' LD HL,A+0178 ***** 033D' CALL $T3 ***** 0340' LD HL,0030 ***** 0343' EX DE,HL ***** 0344' LD HL,(XOLD) ***** 0347' ADD HL,DE ***** 0348' LD (T:0D0000),HL ***** 034B' LD HL,[ 00 08 27 8E] ***** 034E' CALL $L1 ***** 0351' CALL $CK ***** 0354' LD HL,(T:0D0000) ***** 0357' CALL $DU ***** 035A' CALL $NU ***** 035D' LD HL,A+01D0 ***** 0360' CALL $T3 ***** 0363' LD HL,0028 ***** 0366' EX DE,HL ***** 0367' LD HL,(XOLD) ***** 036A' ADD HL,DE ***** 036B' LD (TEX DE,HL ***** 03D3' LD HL,(XOLD) ***** 03D6' ADD HL,DE ***** 03D7' LD (T:130000),HL ***** 03DA' LD HL,0028 ***** 03DD' EX DE,HL ***** 03DE' LD HL,(XOLD) ***** 03E1' ADD HL,DE ***** 03E2' CALL $L3 ***** 03E5' LD HL,[ 00 08 27 8E] ***** 03E8' CALL $MR ***** 03EB' LD HL,(T:130000) ***** 03EE' CALL $DU ***** 03F1' LD HL,T:070003 ***** 03F4' CALL $T3 ***** 03F7' LD HL,[ 00 D0 52 8E] ***** 03FA' CALL $L1 ***** 03FD' CALL $CK ***** 0400' LD HL,(T:130000) ***** 0403' CALL $DU ***** 0406' LD HL,  0066' CALL $L1 ***** 0069' CALL $CK ***** 006C' LD HL,SIMUL ***** 006F' CALL $T3 ***** 0072' RET 12 5 DETER = 1. 13 DO 18 K=1,N ***** 0073' LD HL,[ 00 00 00 81] ***** 0076' CALL $L1 ***** 0079' CALL $CK ***** 007C' LD HL,DETER ***** 007F' CALL $T3 14 KM1=K-1 ***** 0082' LD HL,0001 ***** 0085' LD (K),HL 15 PIVOT = 0. 16 DO 11 I=1,N ***** 0088' LD HL,(K) ***** 008B' DEC HL ***** 008C' LD (KM1),HL ***** 008F' LD HL,[ 00 00 00 00] ***** 0092' CALL $L1 ****3' LD A,D ***** 00E4' SBC H ***** 00E5' LD H,A ***** 00E6' LD A,L ***** 00E7' RLCA ***** 00E8' OR L ***** 00E9' AND 7F ***** 00EB' OR H ***** 00EC' SUB 01 ***** 00EE' SBC A ***** 00EF' LD (T:000002),A 22 IF (J.EQ.JCOL(JSCAN)) GO TO 11 ***** 00F2' OR A ***** 00F3' JP NZ,11L ***** 00F6' LD HL,(JSCAN) ***** 00F9' ADD HL,HL ***** 00FA' LD DE,JCOL-0002 ***** 00FD' ADD HL,DE ***** 00FE' LD (T:000000),HL ***** 0101' LD HL,(J) ***** 0104' EX DE,HL ***** 0105' LD HL,(T:000000) ***** 0108HL ***** 014B' LD HL,(J) ***** 014E' EX DE,HL ***** 014F' LD HL,0058 ***** 0152' CALL $M9 ***** 0155' EX DE,HL ***** 0156' LD HL,(T:000000) ***** 0159' ADD HL,DE ***** 015A' LD DE,A-0060 ***** 015D' ADD HL,DE ***** 015E' LD (T:010000),HL ***** 0161' LD HL,PIVOT ***** 0164' CALL DABS ***** 0167' LD HL,T:000003 ***** 016A' CALL $T3 ***** 016D' LD HL,(T:010000) ***** 0170' CALL DABS ***** 0173' LD HL,T:000003 ***** 0176' CALL $SU ***** 0179' DEC A ***** 017A' ADD 81 ***** 017C' SBC A ****[ R)Y:%6ʀZ2fԀ!ftͪBd#uH='UdU"H\ PBX3khU!Ak tfр&mPmkm@-`*HE[(!, * k4h*V=MZF>2` v*HE[(!, * f3m6p63@Vf!"@dlTzNϠuc?-Ve`)DZUbJZ* d E"kXcB>2; B[A U8b*@7XT:iSqT"b`#u(Ff7iSqt@ZF`Pouҧ€@ kH\sXd#3[@U$L3 ***** 01A3' LD HL,PIVOT ***** 01A6' CALL $T3 ***** 01A9' LD HL,(K) ***** 01AC' ADD HL,HL ***** 01AD' LD (T:010000),HL ***** 01B0' LD DE,IROW-0002 ***** 01B3' ADD HL,DE ***** 01B4' LD (T:020000),HL ***** 01B7' LD HL,(I) ***** 01BA' EX DE,HL ***** 01BB' LD HL,(T:020000) ***** 01BE' LD (HL),E ***** 01BF' INC HL ***** 01C0' LD (HL),D ***** 01C1' LD HL,(T:010000) ***** 01C4' LD DE,JCOL-0002 ***** 01C7' ADD HL,DE ***** 01C8' LD (T:030000),HL ***** 01CB' LD HL,(J) ***** 01CE' EX DE,HL ****0 00] ***** 0217' CALL $L1 ***** 021A' CALL $CK ***** 021D' LD HL,SIMUL ***** 0220' CALL $T3 ***** 0223' RET 32 13 IROWK = IROW(K) 33 JCOLK = JCOL(K) 34 DETER = DETER*PIVOT 35 DO 14 J= 1, MAX ***** 0224' LD HL,(K) ***** 0227' ADD HL,HL ***** 0228' LD (T:000000),HL ***** 022B' LD DE,IROW-0002 ***** 022E' ADD HL,DE ***** 022F' LD A,(HL) ***** 0230' INC HL ***** 0231' LD H,(HL) ***** 0232' LD L,A ***** 0233' LD (IROWK),HL ***** 0236' LD HL,(T:000000) ***** 0239' LD,KM1 ***** 00BE' LD HL,0001 ***** 00C1' LD (ISCAN),HL 21 IF ( I.EQ.IROW(ISCAN)) GO TO 11 ***** 00C4' LD HL,0001 ***** 00C7' LD (JSCAN),HL ***** 00CA' LD HL,(ISCAN) ***** 00CD' ADD HL,HL ***** 00CE' LD DE,IROW-0002 ***** 00D1' ADD HL,DE ***** 00D2' LD (T:000000),HL ***** 00D5' LD HL,(I) ***** 00D8' EX DE,HL ***** 00D9' LD HL,(T:000000) ***** 00DC' LD A,(HL) ***** 00DD' INC HL ***** 00DE' LD H,(HL) ***** 00DF' LD L,A ***** 00E0' LD A,E ***** 00E1' SUB L ***** 00E2' LD L,A ***** 00E HL ***** 0126' EX DE,HL ***** 0127' LD HL,(KM1) ***** 012A' EX DE,HL ***** 012B' LD A,E ***** 012C' SUB L ***** 012D' LD A,D ***** 012E' SBC H ***** 012F' JP P,00C7' ***** 0132' LD HL,(ISCAN) ***** 0135' INC HL ***** 0136' EX DE,HL ***** 0137' LD HL,(KM1) ***** 013A' EX DE,HL ***** 013B' LD A,E ***** 013C' SUB L ***** 013D' LD A,D ***** 013E' SBC H ***** 013F' JP P,00C1' ***** 0142' LD HL,(I) ***** 0145' ADD HL,HL ***** 0146' ADD HL,HL ***** 0147' ADD HL,HL ***** 0148' LD (T:000000),Ug:Z5)["րEZR`kX\fcy*o=' :Z5)[VsdEZ*{JI9s֍` |e+*7XW:iSqqZFjY*zN<U"H\ BX3@X\B0M 3@fVfPf2>2` v kHR"d ڦkGY 3kXV DZ#"-*V9Btk Xcg)Yb ~TyZpVX#3g8 *ͮ "ͮ#ֱ>2` v! &Ͳʹ,V DZ#'"polf fͤ BE?O31E[HR)pdE[@R)Y:%6:Z2mkmkU@31{J+Z !ZU"H\ BX3op+U!AkVP66 *=OIǕk U@31{J+N !ZBY DZ#'"ۀpomVE!*kVE"-t* `ͥαV }ePZRAZF`Pouҧʾ@ kHR"d z kGDZU"HG Y$U<UnV?Y*zN<@0Xb ~TyW`a fPqfѰvLb-Rp #@ *$);# Q Nct3P  LD HL,(N) ***** 01F1' LD A,(HL) ***** 01F2' INC HL ***** 01F3' LD H,(HL) ***** 01F4' LD L,A ***** 01F5' EX DE,HL ***** 01F6' LD A,E ***** 01F7' SUB L ***** 01F8' LD A,D ***** 01F9' SBC H ***** 01FA' JP P,00A1' ***** 01FD' LD HL,PIVOT ***** 0200' CALL DABS ***** 0203' LD HL,(EPS) ***** 0206' CALL $SU ***** 0209' DEC A ***** 020A' SUB 7F ***** 020C' SBC A ***** 020D' LD (T:000002),A 30 SIMUL = 0. ***** 0210' OR A ***** 0211' JP NZ,13L 31 RETURN ***** 0214' LD HL,[ 00 00 0* 0095' CALL $CK ***** 0098' LD HL,PIVOT ***** 009B' CALL $T3 17 DO 11 J=1,N ***** 009E' LD HL,0001 ***** 00A1' LD (I),HL 18 IF (K.EQ.1) GO TO 9 ***** 00A4' LD HL,0001 ***** 00A7' LD (J),HL ***** 00AA' LD HL,(K) ***** 00AD' DEC HL ***** 00AE' LD A,L ***** 00AF' RLCA ***** 00B0' OR L ***** 00B1' AND 7F ***** 00B3' OR H ***** 00B4' SUB 01 ***** 00B6' SBC A ***** 00B7' LD (T:000002),A 19 DO 8 ISCAN = 1,KM1 ***** 00BA' OR A ***** 00BB' JP NZ,9L 20 DO 8 JSCAN = 1' LD A,(HL) ***** 0109' INC HL ***** 010A' LD H,(HL) ***** 010B' LD L,A ***** 010C' LD A,E ***** 010D' SUB L ***** 010E' LD L,A ***** 010F' LD A,D ***** 0110' SBC H ***** 0111' LD H,A ***** 0112' LD A,L ***** 0113' RLCA ***** 0114' OR L ***** 0115' AND 7F ***** 0117' OR H ***** 0118' SUB 01 ***** 011A' SBC A ***** 011B' LD (T:000002),A 23 8 CONTINUE ***** 011E' OR A ***** 011F' JP NZ,11L 24 9 IF (DABS(A(I,J)).LE.DABS(PIVOT)) GO TO 11 ***** 0122' LD HL,(JSCAN) ***** 0125' INC* 017D' LD (T:000002),A 25 PIVOT = A(I,J) ***** 0180' OR A ***** 0181' JP NZ,11L 26 IROW(K) = I 27 JCOL(K) = J 28 11 CONTINUE ***** 0184' LD HL,(I) ***** 0187' ADD HL,HL ***** 0188' ADD HL,HL ***** 0189' ADD HL,HL ***** 018A' LD (T:000000),HL ***** 018D' LD HL,(J) ***** 0190' EX DE,HL ***** 0191' LD HL,0058 ***** 0194' CALL $M9 ***** 0197' EX DE,HL ***** 0198' LD HL,(T:000000) ***** 019B' ADD HL,DE ***** 019C' LD DE,A-0060 ***** 019F' ADD HL,DE ***** 01A0' CALL 92~-* i:Z#9 HFf7@W1hZn[HR)~dE[@R)T:%6Z2mkm *=OIǕ-*+m`"`#mk HGIDZU!&kV?X]fcy*o=' :Z[VtdE31EZR`?X\g!VuFrm@jL{ZF`Pouҧ *TyZ`AVu@@h)u*o=' :Z1[ f`Fmm|HFf7@W1h=c S-Vaͮͯ(k "`)kk fcxh G* 01CF' LD HL,(T:030000) ***** 01D2' LD (HL),E ***** 01D3' INC HL ***** 01D4' LD (HL),D 29 IF (DABS(PIVOT) .GT. EPS) GO TO 13 ***** 01D5' LD HL,(J) ***** 01D8' INC HL ***** 01D9' EX DE,HL ***** 01DA' LD HL,(N) ***** 01DD' LD A,(HL) ***** 01DE' INC HL ***** 01DF' LD H,(HL) ***** 01E0' LD L,A ***** 01E1' EX DE,HL ***** 01E2' LD A,E ***** 01E3' SUB L ***** 01E4' LD A,D ***** 01E5' SBC H ***** 01E6' JP P,00A7' ***** 01E9' LD HL,(I) ***** 01EC' INC HL ***** 01ED' EX DE,HL ***** 01EE'  4BA' ADD HL,DE ***** 04BB' LD (T:000000),HL ***** 04BE' LD HL,(I) ***** 04C1' ADD HL,HL ***** 04C2' LD DE,JORD-0002 ***** 04C5' ADD HL,DE ***** 04C6' LD (T:010000),HL ***** 04C9' LD HL,(T:000000) ***** 04CC' LD A,(HL) ***** 04CD' INC HL ***** 04CE' LD H,(HL) ***** 04CF' LD L,A ***** 04D0' EX DE,HL ***** 04D1' LD HL,(T:010000) ***** 04D4' LD A,(HL) ***** 04D5' INC HL ***** 04D6' LD H,(HL) ***** 04D7' LD L,A ***** 04D8' LD A,E ***** 04D9' SUB L ***** 04DA' LD L,A ***** 04DB' LD A,D *****0000) ***** 051D' LD (HL),E ***** 051E' INC HL ***** 051F' LD (HL),D ***** 0520' LD HL,(INTCH) ***** 0523' INC HL ***** 0524' LD (INTCH),HL 61 IF ( INTCH/2*2.NE.INTCH ) DETER = -DETER ***** 0527' LD HL,(J) ***** 052A' INC HL ***** 052B' EX DE,HL ***** 052C' LD HL,(N) ***** 052F' LD A,(HL) ***** 0530' INC HL ***** 0531' LD H,(HL) ***** 0532' LD L,A ***** 0533' EX DE,HL ***** 0534' LD A,E ***** 0535' SUB L ***** 0536' LD A,D ***** 0537' SBC H ***** 0538' JP P,04B0' ***** 053B' LD 7F ***** 058B' OR H ***** 058C' DEC A ***** 058D' ADD 81 ***** 058F' SBC A ***** 0590' LD (T:000002),A 63 SIMUL = DETER ***** 0593' OR A ***** 0594' JP NZ,26L 64 RETURN ***** 0597' LD HL,DETER ***** 059A' CALL $L3 ***** 059D' LD HL,SIMUL ***** 05A0' CALL $T3 ***** 05A3' RET 65 26 DO 28 J=1,N 66 DO 27 I=1,N ***** 05A4' LD HL,0001 ***** 05A7' LD (J),HL 67 IROWI = IROW(I) ***** 05AA' LD HL,0001 ***** 05AD' LD (I),HL 68 JCOLI = JCOL(I) 69 27 Y(JC05FC' LD HL,(T:000000) ***** 05FF' CALL $T3 ***** 0602' LD HL,(I) ***** 0605' INC HL ***** 0606' EX DE,HL ***** 0607' LD HL,(N) ***** 060A' LD A,(HL) ***** 060B' INC HL ***** 060C' LD H,(HL) ***** 060D' LD L,A ***** 060E' EX DE,HL ***** 060F' LD A,E ***** 0610' SUB L ***** 0611' LD A,D ***** 0612' SBC H ***** 0613' JP P,05AD' 71 28 A(I,J) = Y(I) ***** 0616' LD HL,0001 ***** 0619' LD (I),HL 72 DO 30 I=1,N ***** 061C' LD HL,(I) ***** 061F' ADD HL,HL ***** 0620' ADD HL,HL **** 0670' JP P,05A7' 73 DO 29 J=1,N ***** 0673' LD HL,0001 ***** 0676' LD (I),HL 74 IROWJ= IROW(J) ***** 0679' LD HL,0001 ***** 067C' LD (J),HL 75 JCOLJ = JCOL(J) 76 29 Y(IROWJ) = A(I,JCOLJ) ***** 067F' LD HL,(J) ***** 0682' ADD HL,HL ***** 0683' LD (T:000000),HL ***** 0686' LD DE,IROW-0002 ***** 0689' ADD HL,DE ***** 068A' LD A,(HL) ***** 068B' INC HL ***** 068C' LD H,(HL) ***** 068D' LD L,A ***** 068E' LD (IROWJ),HL ***** 0691' LD HL,(T:000000) ***** 0694' LD DE,1' SBC H ***** 06E2' JP P,067C' 78 30 A(I,J) = Y(J) ***** 06E5' LD HL,0001 ***** 06E8' LD (J),HL 79 SIMUL=DETER ***** 06EB' LD HL,(I) ***** 06EE' ADD HL,HL ***** 06EF' ADD HL,HL ***** 06F0' ADD HL,HL ***** 06F1' LD (T:000000),HL ***** 06F4' LD HL,(J) ***** 06F7' EX DE,HL ***** 06F8' LD HL,0058 ***** 06FB' CALL $M9 ***** 06FE' EX DE,HL ***** 06FF' LD HL,(T:000000) ***** 0702' ADD HL,DE ***** 0703' LD DE,A-0060 ***** 0706' ADD HL,DE ***** 0707' LD (T:010000),HL ***** 070A' LD H 04F9' LD A,(HL) ***** 04FA' INC HL ***** 04FB' LD H,(HL) ***** 04FC' LD L,A ***** 04FD' LD (JTEMP),HL ***** 0500' LD HL,(I) ***** 0503' ADD HL,HL ***** 0504' LD DE,JORD-0002 ***** 0507' ADD HL,DE ***** 0508' LD (T:010000),HL ***** 050B' LD A,(HL) ***** 050C' INC HL ***** 050D' LD H,(HL) ***** 050E' LD L,A ***** 050F' EX DE,HL ***** 0510' LD HL,(T:000000) ***** 0513' LD (HL),E ***** 0514' INC HL ***** 0515' LD (HL),D ***** 0516' LD HL,(JTEMP) ***** 0519' EX DE,HL ***** 051A' LD HL,(T:0161' RLCA ***** 0562' OR L ***** 0563' AND 7F ***** 0565' OR H ***** 0566' ADD FF ***** 0568' SBC A ***** 0569' LD (T:000002),A ***** 056C' OR A ***** 056D' JP Z,0000 62 IF (INDIC.LE.0) GO TO 26 ***** 0570' LD HL,DETER ***** 0573' CALL $L3 ***** 0576' CALL $NU ***** 0579' LD HL,DETER ***** 057C' CALL $T3 ***** 057F' LD HL,(INDIC) ***** 0582' LD A,(HL) ***** 0583' INC HL ***** 0584' LD H,(HL) ***** 0585' LD L,A ***** 0586' LD A,L ***** 0587' RLCA ***** 0588' OR L ***** 0589' AND  05D3' ADD HL,HL ***** 05D4' ADD HL,HL ***** 05D5' ADD HL,HL ***** 05D6' LD DE,Y-0008 ***** 05D9' ADD HL,DE ***** 05DA' LD (T:000000),HL ***** 05DD' LD HL,(IROWI) ***** 05E0' ADD HL,HL ***** 05E1' ADD HL,HL ***** 05E2' ADD HL,HL ***** 05E3' LD (T:010000),HL ***** 05E6' LD HL,(J) ***** 05E9' EX DE,HL ***** 05EA' LD HL,0058 ***** 05ED' CALL $M9 ***** 05F0' EX DE,HL ***** 05F1' LD HL,(T:010000) ***** 05F4' ADD HL,DE ***** 05F5' LD DE,A-0060 ***** 05F8' ADD HL,DE ***** 05F9' CALL $L3 *****  ***** 0650' LD HL,(N) ***** 0653' LD A,(HL) ***** 0654' INC HL ***** 0655' LD H,(HL) ***** 0656' LD L,A ***** 0657' EX DE,HL ***** 0658' LD A,E ***** 0659' SUB L ***** 065A' LD A,D ***** 065B' SBC H ***** 065C' JP P,0619' ***** 065F' LD HL,(J) ***** 0662' INC HL ***** 0663' EX DE,HL ***** 0664' LD HL,(N) ***** 0667' LD A,(HL) ***** 0668' INC HL ***** 0669' LD H,(HL) ***** 066A' LD L,A ***** 066B' EX DE,HL ***** 066C' LD A,E ***** 066D' SUB L ***** 066E' LD A,D ***** 066F' SBC H ****HL ***** 06B9' LD HL,0058 ***** 06BC' CALL $M9 ***** 06BF' EX DE,HL ***** 06C0' LD HL,(T:010000) ***** 06C3' ADD HL,DE ***** 06C4' LD DE,A-0060 ***** 06C7' ADD HL,DE ***** 06C8' CALL $L3 ***** 06CB' LD HL,(T:000000) ***** 06CE' CALL $T3 ***** 06D1' LD HL,(J) ***** 06D4' INC HL ***** 06D5' EX DE,HL ***** 06D6' LD HL,(N) ***** 06D9' LD A,(HL) ***** 06DA' INC HL ***** 06DB' LD H,(HL) ***** 06DC' LD L,A ***** 06DD' EX DE,HL ***** 06DE' LD A,E ***** 06DF' SUB L ***** 06E0' LD A,D ***** 06E 04DC' SBC H ***** 04DD' LD H,A ***** 04DE' LD A,L ***** 04DF' RLCA ***** 04E0' OR L ***** 04E1' AND 7F ***** 04E3' OR H ***** 04E4' RLA ***** 04E5' CCF ***** 04E6' SBC A ***** 04E7' LD (T:000002),A 56 JTEMP = JORD(J) ***** 04EA' OR A ***** 04EB' JP NZ,22L 57 JORD(J)=JORD(I) 58 JORD(I)=JTEMP 59 INTCH = INTCH + 1 60 22 CONTINUE ***** 04EE' LD HL,(J) ***** 04F1' ADD HL,HL ***** 04F2' LD DE,JORD-0002 ***** 04F5' ADD HL,DE ***** 04F6' LD (T:000000),HL *****HL,(I) ***** 053E' INC HL ***** 053F' EX DE,HL ***** 0540' LD HL,(NM1) ***** 0543' EX DE,HL ***** 0544' LD A,E ***** 0545' SUB L ***** 0546' LD A,D ***** 0547' SBC H ***** 0548' JP P,04A6' ***** 054B' LD HL,(INTCH) ***** 054E' EX DE,HL ***** 054F' LD HL,0002 ***** 0552' CALL $D9 ***** 0555' ADD HL,HL ***** 0556' EX DE,HL ***** 0557' LD HL,(INTCH) ***** 055A' LD A,E ***** 055B' SUB L ***** 055C' LD L,A ***** 055D' LD A,D ***** 055E' SBC H ***** 055F' LD H,A ***** 0560' LD A,L ***** 05OLI) = A(IROWI,J) ***** 05B0' LD HL,(I) ***** 05B3' ADD HL,HL ***** 05B4' LD (T:000000),HL ***** 05B7' LD DE,IROW-0002 ***** 05BA' ADD HL,DE ***** 05BB' LD A,(HL) ***** 05BC' INC HL ***** 05BD' LD H,(HL) ***** 05BE' LD L,A ***** 05BF' LD (IROWI),HL ***** 05C2' LD HL,(T:000000) ***** 05C5' LD DE,JCOL-0002 ***** 05C8' ADD HL,DE ***** 05C9' LD A,(HL) ***** 05CA' INC HL ***** 05CB' LD H,(HL) ***** 05CC' LD L,A ***** 05CD' LD (JCOLI),HL 70 DO 28 I = 1,N ***** 05D0' LD HL,(JCOLI) ******* 0621' ADD HL,HL ***** 0622' LD (T:000000),HL ***** 0625' LD HL,(J) ***** 0628' EX DE,HL ***** 0629' LD HL,0058 ***** 062C' CALL $M9 ***** 062F' EX DE,HL ***** 0630' LD HL,(T:000000) ***** 0633' ADD HL,DE ***** 0634' LD DE,A-0060 ***** 0637' ADD HL,DE ***** 0638' LD (T:010000),HL ***** 063B' LD HL,(T:000000) ***** 063E' LD DE,Y-0008 ***** 0641' ADD HL,DE ***** 0642' CALL $L3 ***** 0645' LD HL,(T:010000) ***** 0648' CALL $T3 ***** 064B' LD HL,(I) ***** 064E' INC HL ***** 064F' EX DE,HL JCOL-0002 ***** 0697' ADD HL,DE ***** 0698' LD A,(HL) ***** 0699' INC HL ***** 069A' LD H,(HL) ***** 069B' LD L,A ***** 069C' LD (JCOLJ),HL 77 DO 30 J=1,N ***** 069F' LD HL,(IROWJ) ***** 06A2' ADD HL,HL ***** 06A3' ADD HL,HL ***** 06A4' ADD HL,HL ***** 06A5' LD DE,Y-0008 ***** 06A8' ADD HL,DE ***** 06A9' LD (T:000000),HL ***** 06AC' LD HL,(I) ***** 06AF' ADD HL,HL ***** 06B0' ADD HL,HL ***** 06B1' ADD HL,HL ***** 06B2' LD (T:010000),HL ***** 06B5' LD HL,(JCOLJ) ***** 06B8' EX DE,  K) JCOLK = JCOL(K) DETER = DETER*PIVOT DO 14 J= 1, MAX 14 A(IROWK,J) = A(IROWK,J)/PIVOT A(IROWK,JCOLK) = 1./PIVOT DO 18 I = 1,N AIJCK = A(I,JCOLK) IF ( I.EQ.IROWK) GO TO 18 A(I,JCOLK) = -AIJCK/PIVOT DO 17 J=1,MAX 17 IF (J.NE.JCOLK ) A(I,J) =A(I,J) - AIJCK*A(IROWK,J) 18 CONTINUE DO 20 I=1,N IROWI = IROW(I) JCOLI = JCOL(I) JORD(IROWI) = JCOLI 20 IF ( INDIC.GE.0 ) X(JCOLI) = A(IROWI,MAX)  DO 8 JSCAN = 1,KM1 IF ( I.EQ.IROW(ISCAN)) GO TO 11 IF (J.EQ.JCOL(JSCAN)) GO TO 11 8 CONTINUE 9 IF (DABS(A(I,J)).LE.DABS(PIVOT)) GO TO 11 PIVOT = A(I,J) IROW(K) = I JCOL(K) = J 11 CONTINUE IF (DABS(PIVOT) .GT. EPS) GO TO 13 SIMUL = 0. RETURN 13 IROWK = IROW(K) JCOLK = JCOL(K) DETER = DETER*PIVOT DO 14 J= 1, MAX 14 A(IROWK,J) = A(IROWK,J)/PIVOT A(IROWK,JCOLK) = 1./PIVOT DO 18 I DITIONED OR SINGULAR ) 202 FORMAT ( 10H0ITER = ,I8/ 10H DETER = , E18.5/ 2 26H XOLD(1)...XOLD(, I2, 1H) / (1H ,1P4E16.6) ) 203 FORMAT ( 24H0SUCCESSFUL CONVERGENCE / 10H0ITER = , I3/ 2 26H0 XOLD(1)...XOLD(, I2, 1H) / 1H / (1H ,1P4E16.6) ) 204 FORMAT ( 15H NO CONVERGENCE ) 300 FORMAT (5X,'CHEMICAL EQUILIBRIUM CALCULATION') 301 FORMAT (10X,'PLEASE INPUT NUMBER OF ITERATIONS') 305 FORMAT (10X,'PLEASE INPUT STARTING VALUES') 306 FORMAT (8E10.4)09.*XOLD(2) + 78213.*XOLD(3) 1 -18927.*XOLD(4) - 8427.*XOLD(5) - 13492./XOLD(7) + 10690. 2 *XOLD(6)/XOLD(7) A(5,1)= 1.0 A(5,2)=1.0 A(5,3)=1.0 A(5,4)=1.0 A(5,5)=1.0 A(5,8)= 1.0-XOLD(1)- XOLD(2)-XOLD(3)-XOLD(4)-XOLD(5) A(6,1)=P*P*XOLD(4)**3 A(6,3)= -1.78375E5*XOLD(5) A(6,4)= 3.0*P*P*XOLD(1)*XOLD(4)**2 A(6,5)= -1.7837E5*XOLD(3) A(6,8)= 1.7837E5*XOLD(5)*XOLD(3)-P*P*XOLD(1)*XOLD(4)**3 A(7,1)=XOLD(3) :<7=:PG!o~6xQ!^4!~ 32P:>//ͨUC<!6  6> /:k*&>2>2:2a2*("**"*,"*""W8@>2W;@>,!5:<2xa5-MgW͵#*"ͼb͵#*"ͼb͵#*"Ø':=@Q G'&y:=#Q'#y Eѯ{~0_:lk<͂#n!~wMAڊ7U>'b͑':=m>!l6 #=¤ɯ2Q<Ý';iU:ɯ222og"0"2͘ 2 OG:9!"#"! 2$͎+! 2$͎+! /$! >#͎+! 2$! >#͎+! 2$͎+! /$!"*)))"*! >#*#*{zi͎+! 2$!"*)))"*! >#*#*{z͎+!"!"!͝ !"! 2$!"*)))"*!Xz#*"*! >#*#>>͎+*#>>l!͂!| !ͅ2y! 2$͎+!"!"*)))"*!tT=2ʯ!"*)))""**d* *#*{z! 2$! >#|! >N 28 A(I,J) = Y(I) DO 30 I=1,N DO 29 J=1,N IROWJ= IROW(J) JCOLJ = JCOL(J) 29 Y(IROWJ) = A(I,JCOLJ) DO 30 J=1,N 30 A(I,J) = Y(J) SIMUL=DETER RETURN 200 FORMAT (10H0N TOO BIG ) END KM1=K-1 PIVOT = 0. DO 11 I=1,N DO 11 J=1,N IF (K.EQ.1) GO TO 9 DO 8 ISCAN = 1,KM1  (1,306) (A(I,J),J=1,8) 2 CONTINUE DETER = SIMUL(N,EPS1,INDIC,NRC) IF (DETER.NE.0) GO TO 3 WRITE (1,201) GO TO 1 3 ITCON = 1 DO 5 I=1,N IF (DABS(XINC(I)).GT.EPS2) ITCON = 0 5 XOLD(I) = XOLD(I) + XINC(I) WRITE(1,202) ITER,DETER,N,(XOLD(I),I=1,N) IF (ITCON.EQ.0) GO TO 9 WRITE (1,203) ITER,N,(XOLD(I),I=1,N) GO TO 1 9 CONTINUE WRITE (1,204) GO TO 1 100 FORMAT ( 2I3) 201 FORMAT ( 38H0MATRIX IS ILL-CON A(2,3)= 1.0 A(2,4)= 1.0 A(2,5)= 2.0 A(2,7)=2.0/XOLD(7)**2 A(2,8)=-XOLD(3)-XOLD(4)-2.0*XOLD(5) + 2.0/XOLD(7) A(3,1)= 1.0 A(3,2)=1.0 A(3,5)=1.0 A(3,7)= 1.0/XOLD(7)**2 A(3,8)=-XOLD(1)-XOLD(2)-XOLD(5) + 1.0/XOLD(7) A(4,1)= -28837. A(4,2)=-139009. A(4,3)= -78213. A(4,4)= 18927. A(4,5)= 8427. A(4,6)= -10690./XOLD(7) A(4,7)= (-13492. + 10690.*XOLD(6))/XOLD(7)**2 A(4,8)= 28837.*XOLD(1) + 1390R = 1. DO 18 K=1,N KM1=K-1 PIVOT = 0. DO 11 I=1,N DO 11 J=1,N IF (K.EQ.1) GO TO 9 DO 8 ISCAN = 1,KM1 DO 8 JSCAN = 1,KM1 IF ( I.EQ.IROW(ISCAN)) GO TO 11 IF (J.EQ.JCOL(JSCAN)) GO TO 11 8 CONTINUE 9 IF (DABS(A(I,J)).LE.DABS(PIVOT)) GO TO 11 PIVOT = A(I,J) IROW(K) = I JCOL(K) = J 11 CONTINUE IF (DABS(PIVOT) .GT. EPS) GO TO 13 SIMUL = 0. RETURN 13 IROWK = IROW(K) !>22>2=2<2ñ .!!xe!yH!2>2=2*|*!*}!*|!":0-´ mô[_GG'p!!!w:P0!!wg":(2I3)(38H0MATRIX IS ILL-CONDITIONED OR SINGULAR)(10H0ITER = ,I8/10H DETER = ,E18.5/26H XOLD(1)...XOLD(,I2,1H))(24H0SUCCESSFUL CONVERGENCE /10H0ITER = ,I3/26H0 XOLD(1)...XOLD(,I2,1H))(15H NO CONVERGENCE)(5X,'CHEMICAL EQUILIBRIUM CALCULATION')(10X,'PLEASE INPUT NUMBER OF ITERATIONS')(10X,'PLEASE INPUT STARTING VALUES')(8E10.4)Q INTCH = 0 NM1 = N-1 DO 22 I=1, NM1 IP1 = I+1 DO 22 J=IP1,N IF (JORD(J).GE.JORD(I)) GO TO 22 JTEMP = JORD(J) JORD(J)=JORD(I) JORD(I)=JTEMP INTCH = INTCH + 1 22 CONTINUE IF ( INTCH/2*2.NE.INTCH ) DETER = -DETER IF (INDIC.LE.0) GO TO 26 SIMUL = DETER RETURN 26 DO 28 J=1,N DO 27 I=1,N IROWI = IROW(I) JCOLI = JCOL(I) 27 Y(JCOLI) = A(IROWI,J) DO 28 I = 1, PROGRAM EQUIL C SEE CARNAHAN P325 DOUBLE PRECISION XOLD,XINC,A,EPS1,EPS2,DETER DIMENSION XOLD(10) COMMON //A(11,11),XINC(10) DATA EPS1,EPS2 /1.D-10,1.D-5/ N=7 NRC=N+1 WRITE (1,300) 1 WRITE (1,301) READ (1,100) ITMAX WRITE (1,100) ITMAX WRITE (1,305) READ (1,306) (XOLD(I),I=1,N) WRITE (1,306) (XOLD(I),I=1,N) INDIC=1 DO 9 ITER = 1,ITMAX CALL CALCN(XOLD,NRC ) DO 2 I=1,7 WRITE END SUBROUTINE CALCN( XOLD, NRC) DOUBLE PRECISION XOLD,A,X DIMENSION XOLD(10) COMMON //A(11,11),X(10) DATA P / 20. / DO 1 I =1,7 DO 1 J=1,8 1 A(I,J) = 0. C COMPUTE NONZERO ELEMENTS OF DERIVATIVE MATRIX A(1,1) = 0.5 A(1,2)=1.0 A(1,3)=0.5 WRITE (1,300) (XOLD(I),I=1,7) 300 FORMAT (7E10.4) A(1,6)= -1.0/XOLD(7) A(1,7)= XOLD(6)/XOLD(7)**2 A(1,8)=-XOLD(1)/2.-XOLD(2)-XOLD(3)/2. + XOLD(6)/XOLD(7)  A(7,2)= -2.6058*XOLD(4) A(7,3) = XOLD(1) A(7,4) = -2.6058*XOLD(2) A(7,8) = 2.6058*XOLD(4)*XOLD(2)-XOLD(1)*XOLD(3) RETURN END FUNCTION SIMUL(N,EPS,INDIC,NRC) C SEE PAGE 290 OF CARNAHAN REAL*8 A,X,EPS,SIMUL, DETER, PIVOT,Y DIMENSION IROW(10), JCOL(10), JORD(10), Y(10) COMMON //A(11,11),X(10) MAX=N IF (INDIC.GE.0) MAX=N+1 IF (N.LE.50) GO TO 5 WRITE (6,200) SIMUL=0. RETURN 5 DETE 22<2>2y.ʯ!:"y' 2=2>2:0 :8 y+ʼ -ʼ xڱ yZH Xz 'w K*|*!*"2<2<2ñ 2o!6yIڋ Oҋ 5>2!yʸ w+Җ $ʖ >2[B:u20:02P2oñ !6:!24͑#ͪ :P02o#`#>2xڐ":PE.!D$!>2e!:]h>2>2g"-]"+a":!2o:!Q!.!> 22>2*|*!*}!*|!":0-*B":ʱ 2o͎Bô  ڐ">2:!<2o:PNʩ!>s!3 "!3   "QY12y1H1|Yt&wz^Pc|u~r119!8>2;!V!͏!ͪ#og#! }}{(qnz ͪ#15"͏!5"ͪ#~#1=x2"1!Ú! #ý!12:;1!/O! #2"112͇"y+F+F+Fw`h|"gyc":;O|g}oxG-|W"}R"Û0ElaOQ"!;~Gx"ƀx1w͕#w+ɷH1y11x{1!;6͢#wg"z" >r2|k1{>2;͛0/2{0G,}0G,::})͘/:<:G!=::;w_'-&2?::<:2:::=:x'2?:!::w+:=:/'*:":!":ͣ+!=:4’'*:":ͣ+â'::})͘/:=:&!<:/<&Fw!:~w2:!;~w<2:K'2;2:2<:2:>2: ,'*:*:DM*;|(+| (";":Å.*:*:F,*:#":*;+";|#.(Å.*$:DM!:5r(*":+}|w(* : ~2;`i"$:>h(> h(K(ʅ( !:4|(::I(::=ʟ(=¦(:<:ʪ(8!:ʻ(8(#°(:;p/*:|/g}/o#":p/>2:::(8!<::>:(::(:B:w~/<::%-::%-*$:z-* :6 "":50::2$:2%:::=-!"":-!"$:Å.!;~V-we-*:+":|e-:;,*;|ʅ.+";:B:2<:2;2=:2:go":":":":::!:w::ʩ-4*:::G*::»-:~#»-":å-::++!"":ɯ2:go":*:*:~# ---: -":2;G0!. !.*:T])))_-*:|,.#x::=.o>g":|x!:~w+ G.A.!:~w+Z.!;W.e.:!:w#y.6 -";¢.,ʅ./,),'h,P°.::2;Å.|.8!";2:!>/.# .8yO!\/::.!H/ ^#V":y/*:":!D2.*"B~^6M~ʐ^ʐ^ 'Pڐ^xoP<*;FORTDAT$65z5H5>5L35!3: ;_~!&3: ;_^#V2<32=3]34@~94: ;_!36*g9 7|©3}ʬ3"?3"g9i3*:|3+)7|! w#w#3ҝ4d7<3 #w#G3w#8:O9<2O9i9i9IDF0MPIRFWITEXDOMLDZLGSQIBTLOBDEISBEINOVCNGLGSSNA2IODTBIRCEF ¯2O92P9!f9"g9*+`i~9!q9*g9> 9> 9P 'ʆVͦV:'<†Vn1:X> eV> eV#>OU> _V!X V4~²V5> 429_*9~ 9#~"9 9!9^#6 #> 9wɁ*w:aG:!+w+*]r++*pr+s+w+ 9W:<LW>2b:*r+s*]++r+s:==‚W:}W>B>}B>B>B>B>4C>BW!~µW4Q}Q*"2!0w W5-M*͒ ~!iY*|W:0 W: %> p.w$:=:/< % %p.=%2=:::p// p/p.=.%r0:=:S%S% ,=G%2=:::e%>p.=^%::G,p.e%!:;w# {%::ʭ%*$:DM*":";* : *:%`i"$:":*":DM!";*:* : %`i"":é%::*;+}|&::%*;"":>++* :%>++*":";* :%w#=%*:+":|%<5&>G!&:60#̀0!=:4::0w>1;&":~:o& w+_&w!&:|„&}„&!=:4#"@:|/W}/_*:":>2;:::})͘/:<:G:;&x2::-&:;&2?:+:;&0,!::5&<&:::&&ͣ+&E:;:'+>T>Fw> #w#w#w>2:8%> %+> %+>O}2 ;O!1: ;$$=_^#V ^#VQ+a+*:|m+8*:|m+8*#!;:w#|+!:;z+:;!:w-::>++!::5´+::*,!:50,*@:F#"@:,!:~!;!::2:+,2::/,:;-,!?:5+.,ͣ++8 ,!:5,+*":*$:{z7,#"$:* :+p!:5;,8:::;ʅ.I-!:~4Y,8 *:":*;":Å.!":*:*:::ʟ,~'ʐ,|(w#},#~'ʃ,":Å.~'ʶ,G,*:#":#ß,*:#":#~'ʥ,Å.::=,!;4!":V-2:….*:+":|ʅ.*:":!:4Å.!:61ch|$1{_zW}o|g1":":x!:~w#41-1@1a18J!;qx!;k1F#~怩O1 4L1go";";:P9<2P9!96P#6#"96 :1=G!3w#£12>327 087080808O3O3O3O3O31s#r#q#p1111";`i";*;*;DM^#V#N#F#xy 2:;:;7T.t+T.r%t+c.r%&2T.T.T.!=:4|/G}/O!>Y2 X27>{_zW}o|g=O2|g}oځ2CZQt2 o-yOzW{_xGÄ2`iN#F#q#p#=™22:":":`i":2:::2!:_2::*:2_~#fo::<2:W::2>3]3ww!3~{2 ;5{ 7~ L32 ;̓1]34i3 w R7#wF7: ;_!1O3p+qq7: ;_!2^#VT*7777777:77> 8/827:7*7> 8> 8*g97|7}7"7"g9$:7'8:":* :=> 8~+ 817> 8 8> 8~0 8> 8#~#8=8_>27͸7C8r8@8@8>8>8>8$2":9 G8*":&* :w"": }G88:":* :=> 9~+ʩ81”8> 9é8> 9~ ʩ80©8> 9#~#9=ë8B!6*[#~#O!9 ^#V"8!86 #6 #6*#6*###6*#6*#6 #6 !8 ~9#+~+~+LUX)!";]#Û0ͪ##q#zQ#zBU#|G||U#>2U#8O_yW!;~/woG}_}W}Oͪ#1||DM!>))Ґ# =ˆ#!;~7w?##wy7O*;*;";";:>#>#>#>2:ͨ2ͼ2~#fo":ͼ2":*:"::;$ͅ%#e-#8>2:2:<2:H$>2:7$2:<7$>2:2:2:~#fo" ;!C:" :~#+V$!":!9##"#::w$go":":":å$!:>͖2*:~#$go":*:~#›$go":*:~#fo":2:2;2;2;2:2;:go":":̓1*:|$-*:":(0-82$:2%:2":2#:<2;::>++r0K(%G::x )!;2=:͸0::p/!:~w:=:h)!:::9)!:~#~+w#;)M)!:4::o:;^)}o::g":p/!)02::=ʁ)=ʌ)8::og":go":":2:>2:͘/:=:-&:=:2::2?:+*)|( )>00):)02;>ɯ2<:2>:2:*-)>2;)+)*p/*$:+"$:!:4͓*:;."*>2>:2=:͓*:;(E2*D7**+F*-K*>2:*y2:*$:+"$:!:4ͳ)|*Ox q*8> G\*::xʆ*/<2=:G!<:~wÌ(͹)::“*t+-1-1r%-1:;!;w#w¸*r%!<::=:wÓ*|( *::2:*})F*T})>2:|(*p/!:~":-/8}2:2B:y 0/:;.:/-*:}2B:*:+":I-8XH(AILEFGD:0'L,$)*))))\0&(L,8%r)*&@''&::G*::‚/:w#‚/":*;|ʅ.I-!:~"0x2:w::=/=/+>2;>w/!:~0/W+~/_s#r>2;72͸0!:~w//̀0/0! ;:z+̀0/!: ;z+!=:4/GA.+/2:!:w#*0!=:4>++*;*$:"$:yx҅.8 *":"$:Å. ,*;+";|]0Å.::G::2=:!:>w!=:5A.t+A.r%ajSX0y1JS\E!;y1/wT] 6]3L36*56i3d7! w]3@5ɯx3(w* :*":&h5w#_5i377x38w* :i3>2":?ɯx32>3:":* :~~4w:":~+515> ͣ650> ̣6> ģ6#=5~ͣ65> ͣ6:>37:=36]36]32<3<2=37x3 w2":2>3g6\6e6*":&* :w"": \6}.68:>3?ɷ!3: ;_~̊65!2/Ɓ_#>Wd7i32>3>w7!3: ;_~64!2_#>Wɇ6>ͣ6!3: ;_~6q7o&66#<6+d7i32W:QW*͒ +X:ڵB>B>4C>BW!~µW4Q}Q*"2!0w W5-M*͒ ~!iY*|W:0 W:   SUBROUTINE EQKON(TX) COMMON /AB2/ T,P,R,XK(8) DIMENSION VAR(6,5) DATA VAR/-7.65, 3.87E-3, 3.43E5, 71959.3, -93.68, 2 -4.8, 2.21E-3, 4.5E5, 53742.2, -66.41, 3 -2.1, 0.6E-3, 5.42E5, 33417.3, -42.264, 4 4.36, -4.93E-3, 3.2E5, 11992.4, 5.42, 5 -1.9, 2.18E-3, 0.7E5, -5095.62, -7.66 6 3.36, -2.62E-3, 5.55E5, -6323.94, 0.551/ T2=ALOG(T) T3=2.*T*T T4=T/2 CR=1.987 DO 10 I=1,6 XK(I)=VAR(I,4)/T + VAR(I,5) -VAR(I,1)*T2 -V IF(T.EQ.0) GO TO 1000 ***** 0027' LD DE,T ***** 002A' LD HL,[ 01 00] ***** 002D' LD A,02 ***** 002F' CALL $I1 ***** 0032' CALL $ND ***** 0035' LD HL,T ***** 0038' CALL $L1 ***** 003B' LD HL,0000 ***** 003E' CALL $SA ***** 0041' SUB 01 ***** 0043' SBC A ***** 0044' LD (T:000002),A 7 300 FORMAT (5X,'SILICON EQUILIBRIUM CONSTANTS') ***** 0047' OR A ***** 0048' JP NZ,1000L 8 310 FORMAT (5X,'PLEASE INPUT TEMPERATURE IN KELVIN') 9 320 FORMAT (F7.1) 10 TX=1. 11 P=.001 12  PROGRAM SIKIL COMMON /AB2/ T,P,R,XK(8) WRITE (1,300) 300 FORMAT (5X,'SILICON EQUILIBRIUM CONSTANTS') TX=1. T=1400. P=.001 R= .01 CALL EQKON(TX) DO 10 I=1,6 WRITE (1,301) I, XK(I) 10 CONTINUE 301 FORMAT(5X,I5,E12.4) END SUBROUTINE EQKON(TX) DIMENSION VAR(6,5) COMMON /AB2/ T,P,R,XK(8) DATA VAR/-7.65, -4.8, -2.1, 4.36, -1.9, 3.36, 2 3.87E-3, 2.21E-3, .6E-3,-4.93E-3, 2.18E-3, -2.62E-3, 3AR(I,2)*T4 -VAR(I,3)/T3 XK(I)=EXP(XK(I)/CR) 10 CONTINUE XK(7)=R XK(8)=P RETURN END Cromemco FORTRAN 1978 v. 03.21 -- Bytes: 24247 Created: 28-Feb-79 1 PROGRAM SIKIL 2 COMMON /AB2/ T,P,R,XK(8) ***** 0000' LD BC,$$L ***** 0003' JP $INIT 3 WRITE (1,300) ***** 0006' LD DE,300L ***** 0009' LD HL,[ 01 00] ***** 000C' CALL $W2 4 5 WRITE (1,310) ***** 000F' CALL $ND ***** 0012' LD DE,310L ***** 0015' LD HL,[ 01 00] ***** 0018' CALL $W2 5 READ (1,320) T ***** 001B' CALL $ND ***** 001E' LD DE,320L ***** 0021' LD HL,[ 01 00] ***** 0024' CALL $R2 6 007E' LD HL,[ 01 00] ***** 0081' CALL $W2 16 10 CONTINUE ***** 0084' LD DE,I ***** 0087' LD HL,[ 01 00] ***** 008A' LD A,02 ***** 008C' CALL $I0 ***** 008F' LD HL,(I) ***** 0092' ADD HL,HL ***** 0093' ADD HL,HL ***** 0094' LD DE,XK-0004 ***** 0097' ADD HL,DE ***** 0098' LD (T:000000),HL ***** 009B' LD HL,(T:000000) ***** 009E' EX DE,HL ***** 009F' LD HL,[ 01 00] ***** 00A2' LD A,02 ***** 00A4' CALL $I1 ***** 00A7' CALL $ND 17 301 FORMAT(5X,I5,E12.4) ***** 00AA' LD HL,(I) ***** 00ADN END AR(I,2)*T4 -VAR(I,3)/T3 XK(I)=EXP(-XK(I)/CR) 10 CONTINUE XK(7)=R XK(8)=P RETURN END  SUBROUTINE EQKON(TX) COMMON /AB2/ T,P,R,XK(8) DIMENSION VAR(6,5) DATA VAR/-7.65, 3.87E-3, 3.43E5, 71959.3, -93.68, 2 -4.8, 2.21E-3, 4.5E5, 53742.2, -66.41, 3 -2.1, 0.6E-3, 5.42E5, 33417.3, -42.264, 4 4.36, -4.93E-3, 3.2E5, 11992.4, 5.42, 5 -1.9, 2.18E-3, 0.7E5, -5095.62, -7.66 6 3.36, -2.62E-3, 5.55E5, -6323.94, 0.551/ T2=ALOG(T) T3=2.*T*T T4=T/2 CR=1.987 DO 10 I=1,6 XK(I)=VAR(I,4)/T + VAR(I,5) -VAR(I,1)*T2 -V R= .01 13 CALL EQKON(TX) 14 DO 10 I=1,6 ***** 004B' LD HL,[ 00 00 00 81] ***** 004E' CALL $L1 ***** 0051' LD HL,TX ***** 0054' CALL $T1 ***** 0057' LD HL,[ 6E 12 03 77] ***** 005A' CALL $L1 ***** 005D' LD HL,P ***** 0060' CALL $T1 ***** 0063' LD HL,[ 0A D7 23 7A] ***** 0066' CALL $L1 ***** 0069' LD HL,R ***** 006C' CALL $T1 ***** 006F' LD HL,TX ***** 0072' CALL EQKON 15 WRITE (1,301) I, XK(I) ***** 0075' LD HL,0001 ***** 0078' LD (I),HL ***** 007B' LD DE,301L *****  3.43E5, 4.5E5, 5.42E5, 3.2E5,.7E5, 5.55E5, 4 71959.3, 53742.2, 33417.3, 11992.4, -5095.62, -6323.94, 5 -93.68, -66.41, -42.264, -5.42, -7.66, 0.551/ T2=ALOG(T) T3=2.*T*T T4=T/2 CR=1.987 DO 10 I=1,6 XK(I)=VAR(I,4)/T + VAR(I,5) -VAR(I,1)*T2 -VAR(I,2)*T4 -VAR(I,3)/T3 XK(I)=EXP(XK(I)/CR) 10 CONTINUE XK(7)=R XK(8)=P DO 100 I=1,6 WRITE (1,300) (VAR(I,J),J=1,5) 100 CONTINUE 300 FORMAT (5X, 5(E14.7, 5X)) RETUR   PROGRAM SIKIL COMMON /AB2/ T,P,R,XK(8) WRITE (1,300) 5 WRITE (1,310) READ (1,320) T IF(T.EQ.0) GO TO 1000 300 FORMAT (5X,'SILICON EQUILIBRIUM CONSTANTS') 310 FORMAT (5X,'PLEASE INPUT TEMPERATURE IN KELVIN') 320 FORMAT (F7.1) TX=1. P=.001 R= .01 CALL EQKON(TX) DO 10 I=1,6 WRITE (1,301) I, XK(I) 10 CONTINUE 301 FORMAT(5X,I5,E12.4) GO TO 5 1000 CONTINUE END SUBROUTINE EQKON(TX) DIMENSIT #2`%X3Iӑc/A!38Ṣ4A%<32I#2AA3TQ92 L,@563Uc)93Dc3 MV2TP0 aQ:QUSL!dDs t"&'3@ a@R-dDD49`t"!@`d"($FA" d(")@@!"j2lSbDE4QM%:ki)0 Hb 4f 3@h!4zpfѐhzR >a1~JaIZX]eXNWN rS (BB `&s3S(syMZ;qDhtZ1sTfBm(8$MFMbDDc%8$MU͢@rå@ |>N0Z~&r~&ZH~$$} vc0QTRTBDd$MMFBD38x$DSQR2RDDX$MAFBDD( $DIVF2DXɀ$DPLXbDD3(ɰ$SUbDED$cP$MUbDDs(ʘ$ERR2D49QRSaTDԔXt&@0$ATIbDDH$NUQPV `bDDU5`j`+1@ H G66hԡ` : uA@>@+$ATRbDd$x $DAC0RD$s($ATIBDdȈ$SAFbDD4 T bDD4U =5BjtĈi H` inʦ@ 6'5V$DTP12T2EC8A%<$DPIO6Iӑ`RDTh 0$DOVF BDDcHA&$DLOG2̠BDD(MEH$DPIT2DY 1UQ2kzd*ݠH{djӰ,HĈ*(ـ*HyĈiɑp,HĈ*(qP*HĈIȈ1p*HdjQ&Hd*Q&Hċ*H(q` H-ɱ0H:(HFQԐ HaN 1@ He 1ְ H|j( Hdݐ HbeĈj*Fq0HgĈ )Q HmȱpH-dqH;Q԰HiCĈJ 1ְHecĈfQ0He[ĈTDEit"&"+F@8`V d@)" 9d FB1""1IbDE4$(M!I6DSHRG@)""Bt"))!GB1""Bi6 2Ȭ*¡ :!qXuH@thB,:<-(fS2D4 j$CX`2D4M92DBLERQS!dD$4t!B@"\tX!GG"t!M"Tt)#GA""b,tl$*3@ f`"B`Vdh3h4axFV`6Y1H( qHa5IȑH (qpH%ɳPSRS!TԔ Ht3@h# Q&2RS # 1Rp 10 CONTINUE XK(7)=R XK(8)=P RETURN END 00A9" COMMON Length /AB2/002C (44) Labels: 10L 00FE' 100L 0172' 300L 00AB" ubroutines Referenced: $I1 A@Ɍ2DX$MFMUSdDD)@h,r 4z*4zh-*^!6<_Ta!T3hxfH wdT)" 9e E1""zjdt"$#FO!""Ld"($FH1"" r e,")F@!"j2le& #M"rd)F1" e<- @2" qQUSTDDQ9P5#5> Bm>p> `V>EiQ8c>Ӑ?*iFc1L? RT?@3@:gٰ+4 ڔ 4eW0@+4 k4hC@`f!u2 P5 V2iNg z@ҮР^iD LFBHT8@@ 2fi[!l4 uUUz@h!T@43@h! BY4f4f3k̀3@"`04} 03l Qc- 5V3( P.@A!2I+ 2Ic, 9V2 /A1b2 0A1f3c153"I6`=93<Uc4 ME3NUT5I3r Q7i30IU65U12IQ7=YpDDyxt"&'@UR\ B3@h!h6j4 B3@.3@)t"#&)@p±@>\$+@z@h<<Hq*2c"@ 2 T#M!J28c$ 2f P$9I522T QI RDEŘY 1A1fjx P H"4Pmf 0D4(8\/!Ѡhr@$2D $DPHFBE5hȀ$DPPARDTh$MMFBDdɐ$MU2DYU`RDE$H MV$DNRMLUPQ RDDDhM!2$DSHFCUPQRDU58I9$SUBFԠ2E5 J$AQTPbDd$x56$MMF PQBDd9J$DNGR0I`RDEhAA$DPPF QBEDh '@:d |z@h4 Dz¢3B!0" 562 P%520I#`A!2<IU$ AA2HIQ$2TUc&AA2 T(52vITS!DDIpt"&5Y2h3@3@h BDd0$MVARGBE4h`$DPHFRDEȐ$DIVF2DX$DFLRRDE$MULFBDhH$DPPFRE5T$iQQӡTE4tt)5Y c tH2 Qc 5YI2 P̣"MpPDD%98t" )RDEhh$DFLRRDEɰ$NURDDDhx$DF42E5X$SAFBDDc($DPLX2DYU`2DU( FDBLEXPGG"*t"GC"*8Y4a`+68Y4 B2hʬ@W!U8Y` ! yF";^>` T}B:*x :*x :x:@s i"Xe|٠Á-V,CFAOꦂVCDA@qUBgฌ 0+ `(@oPU*΂^xĂV4k5<0u[G =e\Vm$M8).3*o^AtZ QeTY\,uXYuwFV elmfаvf۰G@Dæk Vn?xʔX(ʠhX;Z¦ ɂ!+v; YՕqb+X: *R`VP9EqAW8EtADWPuB6 :+vU  HꩃW#r_F HFH HUBK%^2&sxKUl'M eRpO4`(t)"C+ Gk U) iVU *E"ȤXtp(UTBVk #jeWzT>' *|Z BCu°@V dl Px? ̅2BAͲAV dB ᣸&a[f` a^qU;  YT %'eZpA{‡V(eRAXw)|[S6ʄ (4-TvtZɅYl*2e0!]a a@"ŻĀ\} ρs8]@@D1Hq5Ċi F1 HcdF1HgIp HaQPH(fqH8TQ bE4CM!)H$SHJR4MF<[quA,@*&ɎbE4C$SHJR1bE4#H$J1ՑBDdh=Y$RND0  RE$C@#F@&B<hhH3'(0t))GE)"BT2DHQ&$ATRMZfI8FF"3=aP0H`4xdZF4!`eS& j$UeTu%ݙi@ݐ*Q5e[*Ā?,*>3lĀL.~ V$| fp?ٶ\=aذ+/(d#,ͦ $BF2D$P$CLSFLbDU%($OPNFLd  2DTh $FL TBBDT B$FWѠBBD$$IM`2D( ($J4Q2E( A2$PRRDTdx M$SG2ED8 QN$DTBFAIbDED$cQ $DTBF5bDED$c$DAC0IT BD$xI$ARG3`2Dh ($J8Iѥ2 p \.K%` pxD\J.&K%  qX\\./K%Ơ` q].aK%ؠp vH,].hK@%`P v\].mK%ۀ v|].qK@% wX].xK%`@ w].~K@% p \.2D$A22D80$AD2D$hh$BL UT CwƠaꪄ (F":!䀆%7,2`Ja@ꆄ-%UްQʝڪf v.fڐp#$ fP 0C #UdiSq]b*Z @Ȭ!V3j@9REu_쫈ñ"~-V`!]" ""w(Bx@:U2: vU??쫞jm#Ewv &#??쫄Z ??Grj@W!#aXW!#??øC=yQ @ aP a\@UIjʸB8"":4a@ GV!n¬#aV@Ueul*V" :#*2V uPBV*!BU|%B P*; H VX$|-Y ;H: d*qH*1pHXdH1HY dHљ0H dh&HJ (qH9iћH d`fH *ђHdP H.dћHdѐH]dјH dq@Ht dH}dHad&1Hed&q8H d)` Hĉ))1PHc dFQЖHgdFHk dF񝀶HqlĉʈQHĉH*HvH QHOd豚 H5 dQH' d QHndFQHk dHq vHHqЀH dhёPH Ċi)0HcĊi JF10,Hi djHdpY"fHdUHIxFVdan<$T1BD$$MRM1BDx$L1BD$$MAR0BDdH8$ACBD3($MARRE4cQBE4hM$SAF1MD@ݐ:݀c  2 PM: IQ#MpUBDDi@Ԥ3@mf @m 63h6@&: U$@M!)02RT")I2$S #$ pQBDEi@/#E^vzA-*4'|+B\n3CLt$+NBE4(M!H$SHR1MI!Ce8à@ `I^-dC|{;G@!"BIAeEPݐ*m5W-QUBq{ fpO[*?,+$>3k?S2^SDcUz hm`c nz݅vVms, FX Œ7Oi*2 Hڣ o*EWLu`ew W eQ3ڽ`L$dݐ+vH X_ ʎwd2/c P++UHAf@,|hG nrD8LuJd~Z1C(b|G!' %`c6pn62\71сܭVm!f`pU 69gv^wd@  Sǁ. D&{W4 ~T$euڬVhͣ@ wl`i0 ;a*,! T &K ĈihɑH)* H8VU DUH b$EXITѓ%@Q` ph@U#k4,2ca%R:V#a%R<ѓ#"QQe=R0SPRSbEEEQQe=R$LNPTRPbDĔ$iBrA"4 xZ,KUS2,dZPG3@/Ȍl LwuX,p@ *` x$TTYINbEEEH$LNPTRbD5$h$LINBFԑU 2D $A3Q 2D$h 2$POS2ET $AC P̠BD38 $POS2 P`2D8 $CD 2DD8 $DE 2DEJ$DV2E$8 $REC3V§ }ZU0 HuCaFTpeyAYHuS𥳐" "036/[*ƒ?+ >+6'T aXǃj(0 bo+n)l>p.#0- / 1CFA3o+uvA\f0vhͼ@bͼê !'ux#;XWl@D7(uqvf۰GB#H23*+x >+ ͷ`*bõb^~[(8#;ìc?uY*eT`%YK2#dڮՠd jµV \Յkٶt `˭! v0n"n;r:@"x w[H-3o@(¾`㫸[6 8QeUBT,V"l;ͩh" PA$b9SKMh5\Qpk6YA,XEw W;uW$*aA0XF0pBV,UBVU a؀ ?-U1xVyulǬ=ePXH;`؅*/+L#9#+7jZVmP @yQ0 L1FUP0Bm[6܃m`]F40#MCpF{GfHl¢*>!`AFaV P 2G&j\зXDU0RJ3kX0UD„V1qa6 `EhV p P󅹁qQ H c$CXD1Co` 63m8ͰL;b (1иHHњxHHL&1HHeE(fH*HPHgA#=Y:IՑ#"` 2  T I9:" #$ 5J:FI g֐PRE3iP5 p"1PH(fQ HaH(sSSU bD@Ըt-bl(tH^wXTd+@&2P$M-J2: %=1:SSc1%9 2&#1U9Q 2V'SSbDTD)@TPpP D@#! IX:S"`M-IZ2PU RDC)@,#3gd@)"J 2EC5I$MRM1 PS`2D5I6$MAR0 U BD)@ͤ#8 :pͦVm °JVm!f@Ƌ Ԁ(xʡ L* գ*Gc@ ¹`wՉnJ=3VnfѠDuWzGU,pdxvKeWHeU Xͤ`2y_2 ͷ@L۬p/ՇO@"d@ʅ lmCHy^p;i08 ͤ VhhAݒͣ@b*m5s3ƒr@ á%;tp  3@m&YPݐ7djgqUl(DE%iQRDE%c9@T D(p K Wh?L.QP0*/M"`k @64 +#bd H>@Vʯ  ڞï S6Uy U| fp?ٵ=aװȈ$IOERR2D$@$BFBDU%($TTYINbEEEH$DRV3TBDU%)@P8~Vw<  IxFVu 1(F6¡F6¡#8 b #aS`<T Q2?ɯB2:*~~4w:~+ʭ1Ÿ> mí0> m> m#=ʻ~mí> m:7:'Ę'2<27B w221&/*&*w" &}͋:?ɷ!:_~T5!/Ɓ_#>W.32>w7!:_~4!_#>Wɇ>m!:_~6;o&6#<ð+.32'ww!~{2{ ~ 2MRN END VAR/-7.65, -4.8, -2.1, 4.36, -1.9, 3.36, 2 3.87E-3, 2.21E-3, .6E-3,-4.93E-3, 2.18E-3, -2.62E-3, 3 3.43E5, 4.5E5, 5.42E5, 3.2E5,.7E5, 5.55E5, 4 71959.3, 53742.2, 33417.3, 11992.4, -5095.62, -6323.94, 5 -93.68, -66.41, -42.264, +5.42, -7.66, 0.551/ T2=ALOG(T) T3=2.*T*T T4=T/2 CR=1.987 DO 10 I=1,6 XK(I)=VAR(I,4)/T + VAR(I,5) -VAR(I,1)*T2 -VAR(I,2)*T4 -VAR(I,3)/T3 XK(I)=EXP(XK(I)/CR) 10 CONTINUE XK(7)=R XK(8)=P RETEINOVCNGLGSSNA2IODTBIRCEFAs22!0"1*+`iH!;*1> ͟> ͟P this StatemenMissing Do TerminatioCode Output in BLOCK DATUndefined Labels HavK_*~ ̿#~" I!J^#6 #> ͟wand UsagFunction with no ParameteHexConstant OverfloDivision by ZerArray Name ExpecteIllegal Argument to ENCODE/DECODcV:'ˆ;|c;}c;G:'Ÿ;xz;xz;0:eVeVÓ;)))z; Ž;`Ñ;*O;>'eV*O;>"ò;<>N#F#;|c;}c;G:'Ÿ;xAR/-7.65, -4.8, -2.1, 4.36, -1.9, 3.36, 2 3.87E-3, 2.21E-3, .6E-3,-4.93E-3, 2.18E-3, -2.62E-3, 3 3.43E5, 4.5E5, 5.42E5, 3.2E5,.7E5, 5.55E5, 4 71959.3, 53742.2, 33417.3, 11992.4, -5095.62, -6323.94, 5 -93.68, -66.41, -42.264, +5.42, -7.66, 0.551/ T2=ALOG(T) T3=2.*T*T T4=T/2 CR=1.987 DO 10 I=1,6 XK(I)=VAR(I,4)/T + VAR(I,5) -VAR(I,1)*T2 -VAR(I,2)*T4 -VAR(I,3)/T3 XK(I)=EXP(XK(I)/CR) 10 CONTINUE XK(7)=R XK(8)=P RETU PROGRAM SIKIL COMMON /AB2/ T,P,R,XK(8) CALL OPEN (7,SILICON1DAT,0) T=1000. DT=20. TX=1. P=.001 R= .01 DO 30 FOR I=1,26 CALL EQKON(TX) DO 10 I=1,6 WRITE (1,301) I, XK(I) WRITE (7,302) T,I,XK(I) 10 CONTINUE T=T+DT 30 CONTINUE 301 FORMAT(5X,I5,E12.4) 302 FORMAT (F7.1,I2,E12.4) ENDFILE 7 END SUBROUTINE EQKON(TX) DIMENSION VAR(6,5) COMMON /AB2/ T,P,R,XK(8) DATA 'y3 w #w:_!{p+q;:_!^#Vre[:Jʚ> 2J:J*K> > *1o|’}ʕ"K"1:J:*=> ~+1> > ~0> #~#=_>2J͂ <  2ͮ *&*w" }͋:*=> ͟~+s1^> ͟s> ͟~ s0s> ͟#~#͟=uNot Allowe~#O! ^#V"!6 #6 #6*#6*###6*#6*#6 #6 ! ~͟#:<233IDF0MPIRFWITEXDOMLDZLGSQIBTLOBDEISB PROGRAM SIKIL COMMON /AB2/ T,P,R,XK(8) CALL OPEN(7,SILICON1DAT,0) T=1000. DT=20. TX=1. P=.001 R= .01 DO 30 FOR I=1,26 CALL EQKON(TX) DO 10 I=1,6 WRITE (1,301) I, XK(I) WRITE (7,302) T,I,XK(I) 10 CONTINUE T=T+DT 30 CONTINUE 301 FORMAT(5X,I5,E12.4) 302 FORMAT (F7.1,I2,E12.4) ENDFILE 7 END SUBROUTINE EQKON(TX) DIMENSION VAR(6,5) COMMON /AB2/ T,P,R,XK(8) DATA VURN END                   ! " # $ % & ' ( ) * + , - . / 0 1 2 3 4 5 6 7 8 9 : ; < = > ? @ A B C D E F G H I J K L