mirror of
https://github.com/pkimpel/retro-220.git
synced 2026-02-28 17:19:46 +00:00
1. Correct missing relocation pseudo-op in MONIT Library routine. 2. Rebuild Generator and Compiler tapes. 3. Recompile and generate corrected listings for example programs. 4. Add Simpson's Rule example from Algol-58 Preliminary Report.
2454 lines
148 KiB
Plaintext
2454 lines
148 KiB
Plaintext
0200 BAC-220 STANDARD VERSION 2/1/62
|
||
|
||
0200 COMMENT FOURTH EXAMPLE PROGRAM FROM BALGOL MANUAL, MARCH 1963.
|
||
|
||
0200 NOTE THAT THIS REQUIRES A SYSTEM WITH AT LEAST 8000 WORDS.
|
||
|
||
0200
|
||
|
||
0200 THE PROGRAM BELOW HAS BEEN WRITTEN BY G. FORSYTHE, OF STANFORD
|
||
|
||
0200 UNIVERSITY. IT SOLVES A SET OF LINEAR EQUATIONS OF THE FORM AY = B,
|
||
|
||
0200 USING CROUT@S METHOD WITH INTERCHANGES$
|
||
|
||
0200
|
||
|
||
0200 COMMENT FORSYTHE PROGRAM$
|
||
0200 0 0000 30
|
||
0201 0 9999 01 0201
|
||
0202 0 0000 40 9995 TEMP
|
||
0203 0 0000 30
|
||
|
||
0202 PROCEDURE PRODUCT ($ N, A(), P, E)$
|
||
|
||
0202 BEGIN COMMENT THIS FORMS THE PRODUCT OF ARBITRARY FLOATING NUMBERS
|
||
|
||
0202 A(I), FOR I=(1,1,N). EXPONENT OVERFLOW OR UNDERFLOW IS
|
||
|
||
0202 PREVENTED. THE ANSWER IS P TIMES 10*E WHERE E IS 0 IF POSSIBLE.
|
||
|
||
0202 IF E NEQ 0, THEN WE NORMALIZE P SO THAT 0.1 LEQ ABS(P) LSS 1.0$
|
||
|
||
0202 INTEGER E, F, I, K, N$
|
||
|
||
0202
|
||
0204 0 0000 10 9994 CONST
|
||
0205 0 0000 40 9993 Q
|
||
|
||
0204 Q = 1.0**-10$ F = 10$
|
||
0206 0 0000 10 9992 CONST
|
||
0207 0 0000 40 9991 F
|
||
0208 0 0000 10 9990 CONST
|
||
0209 0 0000 40 9989 I
|
||
0210 0 0000 44
|
||
0211 0 0000 30
|
||
0212 0 0000 10 9990 CONST
|
||
0213 0 0000 12 9989 I
|
||
0214 0 0000 40 9989 I
|
||
0211 0215
|
||
0215 0 0000 13 9999 N
|
||
0216 0 1111 31 0218
|
||
0217 0 0000 36
|
||
0218 0 0001 33
|
||
|
||
0208 FOR I = (1,1,N)$
|
||
0219 0 0000 30
|
||
0218 0220
|
||
0217 0220
|
||
|
||
0219 BEGIN
|
||
0220 0 0000 10 9989 I
|
||
0221 0 0000 14 9997
|
||
0222 0 0001 49 0010
|
||
0222 0 0001 40 9988 TEMP
|
||
|
||
0223 IF A(I) EQL 0.0$
|
||
0223 0 0000 42 9988 TEMP
|
||
0224 1 0000 10 9998 A
|
||
0225 0 0000 36 0227
|
||
0226 0 0000 30
|
||
|
||
0223 BEGIN
|
||
|
||
0223 P = 0.0$
|
||
0227 0 0000 46 9996 P
|
||
|
||
0228 E = 0$
|
||
0228 0 0000 46 9995 E
|
||
|
||
0229 RETURN
|
||
0229 0 0000 42 0201
|
||
0230 1 0000 30 0000
|
||
|
||
0231 END$
|
||
|
||
0231
|
||
0226 0231
|
||
0231 0 0000 10 9989 I
|
||
0232 0 0000 14 9997
|
||
0233 0 0001 49 0010
|
||
0233 0 0001 40 9988 TEMP
|
||
0234 0 0000 42 9988 TEMP
|
||
0235 1 0001 10 9998 A
|
||
|
||
0231 IF ABS(A(I)) LEQ 1.0$
|
||
0236 0 0000 13 9987 CONST
|
||
0237 0 1111 31 0239
|
||
0238 0 0000 36 0240
|
||
0239 0 0000 33
|
||
|
||
0231 BEGIN
|
||
|
||
0231 F = F-20$
|
||
0240 0 0000 10 9991 F
|
||
0241 0 0000 13 9986 CONST
|
||
0242 0 0000 40 9991 F
|
||
|
||
0243 Q = Q.(10.0*20)
|
||
0243 0 0000 41 9986 CONST
|
||
0244 0 0000 10 9985 CONST
|
||
0245 0 0000 44 9947 FL*FX
|
||
0246 0 0000 30 9947
|
||
0247 0 0000 24 9993 Q
|
||
0248 0 0000 40 9993 Q
|
||
|
||
0247 END$
|
||
|
||
0247
|
||
0239 0249
|
||
0249 0 0000 10 9989 I
|
||
0250 0 0000 14 9997
|
||
0251 0 0001 49 0010
|
||
|
||
0249 Q = Q.A(I)$
|
||
0251 0 0001 40 9988 TEMP
|
||
0252 0 0000 10 9993 Q
|
||
0253 0 0000 42 9988 TEMP
|
||
0254 1 0000 24 9998 A
|
||
0255 0 0000 40 9993 Q
|
||
|
||
0256 X = ABS(Q)$
|
||
0256 0 0001 10 9993 Q
|
||
0257 0 0000 40 9946 X
|
||
0258 0 0000 11 9992 CONST
|
||
0259 0 0000 40 9945 K
|
||
0260 0 0000 44
|
||
0261 0 0000 30
|
||
0262 0 0000 10 9990 CONST
|
||
0263 0 0000 12 9945 K
|
||
0264 0 0000 40 9945 K
|
||
0261 0265
|
||
0265 0 0000 13 9992 CONST
|
||
0266 0 1111 31 0268
|
||
0267 0 0000 36
|
||
0268 0 0001 33
|
||
0269 0 0000 11 9944 CONST
|
||
0270 0 0000 40 9945 K
|
||
0271 0 0000 44
|
||
0272 0 0000 30
|
||
0273 0 0000 10 9945 K
|
||
0274 0 0000 13 9990 CONST
|
||
0275 0 0000 40 9945 K
|
||
0272 0276
|
||
0276 0 0000 12 9943 CONST
|
||
0277 0 1111 31 0279
|
||
0278 0 0000 36
|
||
0279 0 0000 33
|
||
0280 0 0000 10 9944 CONST
|
||
0281 0 0000 40 9945 K
|
||
0282 0 0000 44
|
||
0283 0 0000 30
|
||
0284 0 0000 10 9990 CONST
|
||
0285 0 0000 12 9945 K
|
||
0286 0 0000 40 9945 K
|
||
0283 0287
|
||
0287 0 0000 13 9943 CONST
|
||
0288 0 1111 31 0290
|
||
0289 0 0000 36
|
||
0290 0 0001 33
|
||
|
||
0258 FOR K = (-10,1,10), (-11,-1,-41), (11,1,41)$
|
||
0291 0 0000 30
|
||
0290 0292
|
||
0289 0292
|
||
0279 0292
|
||
0278 0292
|
||
0268 0292
|
||
0267 0292
|
||
0292 0 0000 41 9945 K
|
||
0293 0 0000 10 9985 CONST
|
||
0294 0 0000 44 9947 FL*FX
|
||
0295 0 0000 30 9947
|
||
0296 0 0000 13 9946 X
|
||
0297 0 1111 31 0299
|
||
0298 0 0000 36 0302
|
||
0299 0 0001 33 0302
|
||
0300 0 0001 45 0002
|
||
0301 0 0000 30 0303
|
||
0302 0 0000 10 9990 CONST
|
||
0303 0 0000 40 9988 TEMP
|
||
0304 0 0000 10 9990 CONST
|
||
0305 0 0000 12 9945 K
|
||
0306 0 0001 48 0010
|
||
0307 0 0000 10 9985 CONST
|
||
0308 0 0000 44 9947 FL*FX
|
||
0309 0 0000 30 9947
|
||
0310 0 0000 13 9946 X
|
||
0311 0 1111 31 0313
|
||
0312 0 0000 36 0316
|
||
0313 0 0001 33 0316
|
||
0314 0 0000 10 9990 CONST
|
||
0315 0 0000 30 0317
|
||
0316 0 0001 45 0002
|
||
0317 0 0000 17 9988 TEMP
|
||
|
||
0291 IF ((10.0*K LEQ X) AND (X LSS 10.0*(K+1)))$
|
||
0318 0 0000 36
|
||
|
||
0291 BEGIN
|
||
0319 0 0000 11 9992 CONST
|
||
0320 0 0000 13 9945 K
|
||
0321 0 0001 48 0010
|
||
0322 0 0000 10 9985 CONST
|
||
0323 0 0000 44 9947 FL*FX
|
||
0324 0 0000 30 9947
|
||
|
||
0321 Q = Q.(10.0*(-10-K))$
|
||
0325 0 0000 24 9993 Q
|
||
0326 0 0000 40 9993 Q
|
||
|
||
0327 F = F + K + 10$
|
||
0327 0 0000 10 9992 CONST
|
||
0328 0 0000 12 9945 K
|
||
0329 0 0000 12 9991 F
|
||
0330 0 0000 40 9991 F
|
||
|
||
0331 GO TO 1
|
||
0331 0 0000 30 0000
|
||
|
||
0332 END$
|
||
0318 0332
|
||
0291 0333
|
||
0282 0332
|
||
0271 0332
|
||
0260 0332
|
||
0332 0 0000 30
|
||
0331 0333
|
||
|
||
0333 1.. END$
|
||
|
||
0333
|
||
0219 0334
|
||
0210 0333
|
||
0333 0 0000 30
|
||
0334 0 0000 11 9942 CONST
|
||
0335 0 0000 13 9991 F
|
||
0336 0 1111 31 0338
|
||
0337 0 0000 36 0341
|
||
0338 0 0001 33 0341
|
||
0339 0 0001 45 0002
|
||
0340 0 0000 30 0342
|
||
0341 0 0000 10 9990 CONST
|
||
0342 0 0000 40 9988 TEMP
|
||
0343 0 0000 10 9991 F
|
||
0344 0 0000 13 9941 CONST
|
||
0345 0 1111 31 0347
|
||
0346 0 0000 36 0350
|
||
0347 0 0001 33 0350
|
||
0348 0 0001 45 0002
|
||
0349 0 0000 30 0351
|
||
0350 0 0000 10 9990 CONST
|
||
0351 0 0000 17 9988 TEMP
|
||
|
||
0334 IF (((-40) LEQ F) AND (F LEQ 58))$
|
||
0352 0 0000 36
|
||
|
||
0334 BEGIN
|
||
0353 0 0000 41 9940 CONST
|
||
0354 0 0000 10 9985 CONST
|
||
0355 0 0000 44 9947 FL*FX
|
||
0356 0 0000 30 9947
|
||
0357 0 0000 24 9993 Q
|
||
0358 0 0000 40 9988 TEMP
|
||
0359 0 0000 10 9991 F
|
||
0360 0 0000 13 9940 CONST
|
||
0361 0 0001 48 0010
|
||
0362 0 0000 10 9985 CONST
|
||
0363 0 0000 44 9947 FL*FX
|
||
0364 0 0000 30 9947
|
||
|
||
0357 P = (Q.(10.0*9)).(10.0*(F-9))$
|
||
0365 0 0000 24 9988 TEMP
|
||
0366 0 0000 40 9996 P
|
||
|
||
0367 E = 0$
|
||
0367 0 0000 46 9995 E
|
||
|
||
0368 RETURN
|
||
0368 0 0000 42 0201
|
||
0369 1 0000 30 0000
|
||
|
||
0370 END$
|
||
|
||
0370
|
||
0352 0370
|
||
0370 0 0000 41 9940 CONST
|
||
0371 0 0000 10 9985 CONST
|
||
0372 0 0000 44 9947 FL*FX
|
||
0373 0 0000 30 9947
|
||
|
||
0370 P = Q.(10.0*9)$
|
||
0374 0 0000 24 9993 Q
|
||
0375 0 0000 40 9996 P
|
||
|
||
0376 E = F - 9$
|
||
0376 0 0000 10 9991 F
|
||
0377 0 0000 13 9940 CONST
|
||
0378 0 0000 40 9995 E
|
||
|
||
0379 RETURN
|
||
0379 0 0000 42 0201
|
||
0380 1 0000 30 0000
|
||
|
||
0381 END PRODUCT()$
|
||
|
||
0381
|
||
0203 0381
|
||
0381 0 0410 40 0378
|
||
0382 0 0410 40 0367
|
||
0383 0 0410 40 0228
|
||
0384 0 0000 10 9996
|
||
0385 0 0410 40 0375
|
||
0386 0 0410 40 0366
|
||
0387 0 0410 40 0227
|
||
0388 0 0000 10 9998
|
||
0389 0 0410 40 0254
|
||
0390 0 0410 40 0235
|
||
0391 0 0410 40 0224
|
||
0392 0 0000 10 9999
|
||
0393 0 0410 40 0215
|
||
0394 0 0000 30 0204
|
||
0200 0395
|
||
0395 0 0000 30
|
||
0396 0 9939 01 0396
|
||
0397 0 0000 40 9934 TEMP
|
||
0398 0 0000 30
|
||
|
||
0395 PROCEDURE INNERPRODUCT (S, F, U(), V())$
|
||
|
||
0395 BEGIN COMMENT THIS FORMS THE INNER PRODUCT OF THE VECTORS
|
||
|
||
0395 U(I) AND V(I) FOR I = (S,1,F)$
|
||
|
||
0395 INTEGER I, S, F$
|
||
|
||
0395
|
||
|
||
0399 SUM = 0.0$
|
||
0399 0 0000 46 9933 SUM
|
||
0400 0 0000 10 9939 S
|
||
0401 0 0000 40 9932 I
|
||
0402 0 0000 44
|
||
0403 0 0000 30
|
||
0404 0 0000 10 9990 CONST
|
||
0405 0 0000 12 9932 I
|
||
0406 0 0000 40 9932 I
|
||
0403 0407
|
||
0407 0 0000 13 9938 F
|
||
0408 0 1111 31 0410
|
||
0409 0 0000 36
|
||
0410 0 0001 33
|
||
|
||
0400 FOR I = (S,1,F)$
|
||
0411 0 0000 30
|
||
0410 0412
|
||
0409 0412
|
||
0412 0 0000 10 9932 I
|
||
0413 0 0000 14 9936
|
||
0414 0 0001 49 0010
|
||
0414 0 0001 40 9931 TEMP
|
||
0415 0 0000 10 9932 I
|
||
0416 0 0000 14 9934
|
||
0417 0 0001 49 0010
|
||
|
||
0411 SUM = SUM + U(I).V(I)$
|
||
0417 0 0001 40 9930 TEMP
|
||
0418 0 0000 42 9931 TEMP
|
||
0419 1 0000 10 9937 U
|
||
0420 0 0000 42 9930 TEMP
|
||
0421 1 0000 24 9935 V
|
||
0422 0 0000 22 9933 SUM
|
||
0423 0 0000 40 9933 SUM
|
||
0411 0425
|
||
0402 0424
|
||
0424 0 0000 30
|
||
|
||
0425 INNERPRODUCT() = SUM$
|
||
0425 0 0000 10 9933 SUM
|
||
|
||
0426 RETURN
|
||
0426 0 0000 42 0396
|
||
0427 1 0000 30 0000
|
||
|
||
0428 END INNERPRODUCT()$
|
||
|
||
0428
|
||
0398 0428
|
||
0428 0 0000 10 9935
|
||
0429 0 0410 40 0421
|
||
0430 0 0000 10 9937
|
||
0431 0 0410 40 0419
|
||
0432 0 0000 30 0399
|
||
0395 0433
|
||
0433 0 0000 30
|
||
0434 0 9929 01 0434
|
||
|
||
0433 PROCEDURE CROUT4 ($ N, A(,), B(), Y(), PIVOT(), DET, EX7$
|
||
0435 0 0000 40 9916 TEMP
|
||
0436 0 0000 30
|
||
|
||
0437 SINGULAR, IP())$
|
||
|
||
0437 BEGIN COMMENT THIS IS CROUTS METHOD WITH INTERCHANGES, TO SOLVE
|
||
|
||
0437 AY = B AND OBTAIN THE TRIANGULAR DECOMPOSITION. IP() STANDS FOR
|
||
|
||
0437 AN INNERPRODUCT ROUTINE THAT MUST BE AVAILABLE WHEN CROUT4() IS
|
||
|
||
0437 CALLED. ALSO, PRODUCT() MUST BE AVAILABLE. THE DETERMINANT OF A
|
||
|
||
0437 IS COMPUTED IN THE FORM DET TIMES 10*EX7, WHERE EX7 IS 0 IF
|
||
|
||
0437 POSSIBLE. IF EX7 NEQ 0, THEN WE NORMALIZE DET WITH 0.1 LEQ
|
||
|
||
0437 ABS(DET) LSS 1$
|
||
|
||
0437 INTEGER K, I, J, IMAX, N, PIVOT$
|
||
|
||
0437 INTEGER EX7$
|
||
|
||
0437
|
||
|
||
0437 INT = 1.0$
|
||
0437 0 0000 10 9987 CONST
|
||
0438 0 0000 40 9915 INT
|
||
0439 0 0000 10 9990 CONST
|
||
0440 0 0000 40 9914 K
|
||
0441 0 0000 44
|
||
0442 0 0000 30
|
||
0443 0 0000 10 9990 CONST
|
||
0444 0 0000 12 9914 K
|
||
0445 0 0000 40 9914 K
|
||
0442 0446
|
||
0446 0 0000 13 9929 N
|
||
0447 0 1111 31 0449
|
||
0448 0 0000 36
|
||
0449 0 0001 33
|
||
|
||
0439 FOR K = (1,1,N)$
|
||
0450 0 0000 30
|
||
0449 0451
|
||
0448 0451
|
||
|
||
0450 BEGIN
|
||
|
||
0450 TEMP = 0$
|
||
0451 0 0000 46 9913 TEMP
|
||
0452 0 0000 10 9914 K
|
||
0453 0 0000 40 9912 I
|
||
0454 0 0000 44
|
||
0455 0 0000 30
|
||
0456 0 0000 10 9990 CONST
|
||
0457 0 0000 12 9912 I
|
||
0458 0 0000 40 9912 I
|
||
0455 0459
|
||
0459 0 0000 13 9929 N
|
||
0460 0 1111 31 0462
|
||
0461 0 0000 36
|
||
0462 0 0001 33
|
||
|
||
0452 FOR I = (K,1,N)$
|
||
0463 0 0000 30
|
||
0462 0464
|
||
0461 0464
|
||
|
||
0463 BEGIN
|
||
0464 0 0000 10 9912 I
|
||
0465 0 0000 14 9927
|
||
0466 0 0001 49 0010
|
||
0467 0 0000 12 9914 K
|
||
0468 0 0000 14 9926
|
||
0469 0 0001 49 0010
|
||
0469 0 0001 40 9911 TEMP
|
||
0470 0 0000 10 9912 I
|
||
0471 0 0000 14 9927
|
||
0472 0 0001 49 0010
|
||
0473 0 0000 12 9914 K
|
||
0474 0 0000 14 9926
|
||
0475 0 0001 49 0010
|
||
0475 0 0001 40 9910 TEMP
|
||
0476 0 0000 10 9990 CONST
|
||
0477 0 4400 28 9916
|
||
0478 1 0000 40 0000
|
||
0479 0 0000 10 9914 K
|
||
0480 0 0000 13 9990 CONST
|
||
0481 1 0000 40 9999
|
||
0482 0 0000 10 9912 I
|
||
0483 0 0000 14 9927
|
||
0484 0 0001 49 0010
|
||
0485 0 0000 14 9926
|
||
0486 0 0001 49 0010
|
||
0487 0 0000 12 0488
|
||
0488 0 0000 01 9928 A
|
||
0489 1 0000 40 9998
|
||
0490 0 0000 10 9926
|
||
0491 1 0000 40 9997
|
||
0492 0 0000 10 9914 K
|
||
0493 0 0000 14 9926
|
||
0494 0 0001 49 0010
|
||
0495 0 0000 12 0496
|
||
0496 0 0000 01 9928 A
|
||
0497 1 0000 40 9996
|
||
0498 0 0000 10 9927
|
||
0499 0 0000 14 9926
|
||
0500 0 0001 49 0010
|
||
|
||
0464 A(I,K) = A(I,K) - IP(1, K-1, A(I,), A(,K))$
|
||
0501 0 0000 44 9916 IP
|
||
0502 0 0500 30 9916
|
||
0503 0 0000 42 9910 TEMP
|
||
0504 1 0000 23 9928 A
|
||
0505 0 0000 40 9910 TEMP
|
||
0506 0 0000 11 9910 TEMP
|
||
0507 0 0000 42 9911 TEMP
|
||
0508 1 0000 40 9928 A
|
||
0509 0 0000 10 9912 I
|
||
0510 0 0000 14 9927
|
||
0511 0 0001 49 0010
|
||
0512 0 0000 12 9914 K
|
||
0513 0 0000 14 9926
|
||
0514 0 0001 49 0010
|
||
0514 0 0001 40 9911 TEMP
|
||
0515 0 0000 42 9911 TEMP
|
||
0516 1 0001 10 9928 A
|
||
|
||
0509 IF ABS(A(I,K)) GTR TEMP$
|
||
0517 0 0000 13 9913 TEMP
|
||
0518 0 1111 31 0520
|
||
0519 0 0000 36
|
||
0520 0 0001 33
|
||
|
||
0509 BEGIN
|
||
0521 0 0000 10 9912 I
|
||
0522 0 0000 14 9927
|
||
0523 0 0001 49 0010
|
||
0524 0 0000 12 9914 K
|
||
0525 0 0000 14 9926
|
||
0526 0 0001 49 0010
|
||
0526 0 0001 40 9911 TEMP
|
||
|
||
0521 TEMP = ABS(A(I,K))$
|
||
0527 0 0000 42 9911 TEMP
|
||
0528 1 0001 10 9928 A
|
||
0529 0 0000 40 9913 TEMP
|
||
|
||
0530 IMAX = I
|
||
|
||
0530 END
|
||
0530 0 0000 10 9912 I
|
||
0531 0 0000 40 9909 IMAX
|
||
0520 0532
|
||
0519 0532
|
||
|
||
0532 END$
|
||
|
||
0532
|
||
0463 0533
|
||
0454 0532
|
||
0532 0 0000 30
|
||
0533 0 0000 10 9914 K
|
||
0534 0 0000 14 9920
|
||
0535 0 0001 49 0010
|
||
0535 0 0001 40 9911 TEMP
|
||
|
||
0533 PIVOT(K) = IMAX$
|
||
0536 0 0000 10 9909 IMAX
|
||
0537 0 0000 42 9911 TEMP
|
||
0538 1 0000 40 9921 PIVOT
|
||
|
||
0539 COMMENT WE HAVE FOUND THAT A(IMAX,K) IS THE LARGEST PIVOT IN COL
|
||
|
||
0539 K. NOW WE INTERCHANGE ROWS K AND IMAX$
|
||
|
||
0539 IF IMAX NEQ K$
|
||
0539 0 0000 10 9909 IMAX
|
||
0540 0 0000 13 9914 K
|
||
0541 0 1111 31 0543
|
||
0542 0 0000 36
|
||
|
||
0539 BEGIN
|
||
|
||
0539 INT = -INT$
|
||
0543 0 0000 11 9915 INT
|
||
0544 0 0000 40 9915 INT
|
||
0545 0 0000 10 9990 CONST
|
||
0546 0 0000 40 9908 J
|
||
0547 0 0000 44
|
||
0548 0 0000 30
|
||
0549 0 0000 10 9990 CONST
|
||
0550 0 0000 12 9908 J
|
||
0551 0 0000 40 9908 J
|
||
0548 0552
|
||
0552 0 0000 13 9929 N
|
||
0553 0 1111 31 0555
|
||
0554 0 0000 36
|
||
0555 0 0001 33
|
||
|
||
0545 FOR J = (1,1,N)$
|
||
0556 0 0000 30
|
||
0555 0557
|
||
0554 0557
|
||
|
||
0556 BEGIN
|
||
0557 0 0000 10 9914 K
|
||
0558 0 0000 14 9927
|
||
0559 0 0001 49 0010
|
||
0560 0 0000 12 9908 J
|
||
0561 0 0000 14 9926
|
||
0562 0 0001 49 0010
|
||
|
||
0557 TEMP = A(K,J)$
|
||
0562 0 0001 40 9911 TEMP
|
||
0563 0 0000 42 9911 TEMP
|
||
0564 1 0000 10 9928 A
|
||
0565 0 0000 40 9913 TEMP
|
||
0566 0 0000 10 9914 K
|
||
0567 0 0000 14 9927
|
||
0568 0 0001 49 0010
|
||
0569 0 0000 12 9908 J
|
||
0570 0 0000 14 9926
|
||
0571 0 0001 49 0010
|
||
0571 0 0001 40 9911 TEMP
|
||
0572 0 0000 10 9909 IMAX
|
||
0573 0 0000 14 9927
|
||
0574 0 0001 49 0010
|
||
0575 0 0000 12 9908 J
|
||
0576 0 0000 14 9926
|
||
0577 0 0001 49 0010
|
||
|
||
0566 A(K,J) = A(IMAX,J)$
|
||
0577 0 0001 40 9910 TEMP
|
||
0578 0 0000 42 9910 TEMP
|
||
0579 1 0000 10 9928 A
|
||
0580 0 0000 42 9911 TEMP
|
||
0581 1 0000 40 9928 A
|
||
0582 0 0000 10 9909 IMAX
|
||
0583 0 0000 14 9927
|
||
0584 0 0001 49 0010
|
||
0585 0 0000 12 9908 J
|
||
0586 0 0000 14 9926
|
||
0587 0 0001 49 0010
|
||
0587 0 0001 40 9911 TEMP
|
||
|
||
0582 A(IMAX,J) = TEMP
|
||
0588 0 0000 10 9913 TEMP
|
||
0589 0 0000 42 9911 TEMP
|
||
0590 1 0000 40 9928 A
|
||
|
||
0591 END$
|
||
|
||
0591
|
||
0556 0592
|
||
0547 0591
|
||
0591 0 0000 30
|
||
0592 0 0000 10 9914 K
|
||
0593 0 0000 14 9924
|
||
0594 0 0001 49 0010
|
||
|
||
0592 TEMP = B(K)$
|
||
0594 0 0001 40 9911 TEMP
|
||
0595 0 0000 42 9911 TEMP
|
||
0596 1 0000 10 9925 B
|
||
0597 0 0000 40 9913 TEMP
|
||
0598 0 0000 10 9914 K
|
||
0599 0 0000 14 9924
|
||
0600 0 0001 49 0010
|
||
0600 0 0001 40 9911 TEMP
|
||
0601 0 0000 10 9909 IMAX
|
||
0602 0 0000 14 9924
|
||
0603 0 0001 49 0010
|
||
|
||
0598 B(K) = B(IMAX)$
|
||
0603 0 0001 40 9910 TEMP
|
||
0604 0 0000 42 9910 TEMP
|
||
0605 1 0000 10 9925 B
|
||
0606 0 0000 42 9911 TEMP
|
||
0607 1 0000 40 9925 B
|
||
0608 0 0000 10 9909 IMAX
|
||
0609 0 0000 14 9924
|
||
0610 0 0001 49 0010
|
||
0610 0 0001 40 9911 TEMP
|
||
|
||
0608 B(IMAX) = TEMP
|
||
0611 0 0000 10 9913 TEMP
|
||
0612 0 0000 42 9911 TEMP
|
||
0613 1 0000 40 9925 B
|
||
|
||
0614 END$
|
||
|
||
0614
|
||
0542 0614
|
||
|
||
0614 COMMENT NOW FOR THE ELIMINATION$
|
||
0614 0 0000 10 9914 K
|
||
0615 0 0000 14 9927
|
||
0616 0 0001 49 0010
|
||
0617 0 0000 12 9914 K
|
||
0618 0 0000 14 9926
|
||
0619 0 0001 49 0010
|
||
0619 0 0001 40 9911 TEMP
|
||
|
||
0614 IF A(K,K) EQL 0$
|
||
0620 0 0000 42 9911 TEMP
|
||
0621 1 0000 10 9928 A
|
||
0622 0 0000 36 0624
|
||
0623 0 0000 30
|
||
|
||
0614 BEGIN
|
||
|
||
0614 DET = 0.0$
|
||
0624 0 0000 46 9919 DET
|
||
|
||
0625 EX7 = 0$
|
||
0625 0 0000 46 9918 EX7
|
||
|
||
0626 GO TO SINGULAR$
|
||
0626 0 0000 30 9917 SINGU
|
||
|
||
0627 END$
|
||
|
||
0627
|
||
0623 0627
|
||
0627 0 0000 10 9990 CONST
|
||
0628 0 0000 12 9914 K
|
||
0629 0 0000 40 9912 I
|
||
0630 0 0000 44
|
||
0631 0 0000 30
|
||
0632 0 0000 10 9990 CONST
|
||
0633 0 0000 12 9912 I
|
||
0634 0 0000 40 9912 I
|
||
0631 0635
|
||
0635 0 0000 13 9929 N
|
||
0636 0 1111 31 0638
|
||
0637 0 0000 36
|
||
0638 0 0001 33
|
||
|
||
0627 FOR I = (K+1,1,N)$
|
||
0639 0 0000 30
|
||
0638 0640
|
||
0637 0640
|
||
|
||
0639 BEGIN
|
||
0640 0 0000 10 9912 I
|
||
0641 0 0000 14 9927
|
||
0642 0 0001 49 0010
|
||
0643 0 0000 12 9914 K
|
||
0644 0 0000 14 9926
|
||
0645 0 0001 49 0010
|
||
|
||
0640 XX = A(I,K)$
|
||
0645 0 0001 40 9911 TEMP
|
||
0646 0 0000 42 9911 TEMP
|
||
0647 1 0000 10 9928 A
|
||
0648 0 0000 40 9907 XX
|
||
0649 0 0000 10 9914 K
|
||
0650 0 0000 14 9927
|
||
0651 0 0001 49 0010
|
||
0652 0 0000 12 9914 K
|
||
0653 0 0000 14 9926
|
||
0654 0 0001 49 0010
|
||
|
||
0649 XY = A(K,K)$
|
||
0654 0 0001 40 9911 TEMP
|
||
0655 0 0000 42 9911 TEMP
|
||
0656 1 0000 10 9928 A
|
||
0657 0 0000 40 9906 XY
|
||
|
||
0658 X = 1.0$
|
||
0658 0 0000 10 9987 CONST
|
||
0659 0 0000 40 9905 X
|
||
|
||
0660 X = X.X$
|
||
0660 0 0000 10 9905 X
|
||
0661 0 0000 24 9905 X
|
||
0662 0 0000 40 9905 X
|
||
0663 0 0000 10 9912 I
|
||
0664 0 0000 14 9927
|
||
0665 0 0001 49 0010
|
||
0666 0 0000 12 9914 K
|
||
0667 0 0000 14 9926
|
||
0668 0 0001 49 0010
|
||
0668 0 0001 40 9911 TEMP
|
||
|
||
0663 A(I,K) = XX/XY
|
||
0669 0 0000 10 9907 XX
|
||
0670 0 0002 45 0000
|
||
0671 0 0000 25 9906 XY
|
||
0672 0 0000 42 9911 TEMP
|
||
0673 1 0000 40 9928 A
|
||
|
||
0674 END$
|
||
|
||
0674
|
||
0639 0675
|
||
0630 0674
|
||
0674 0 0000 30
|
||
0675 0 0000 10 9990 CONST
|
||
0676 0 0000 12 9914 K
|
||
0677 0 0000 40 9908 J
|
||
0678 0 0000 44
|
||
0679 0 0000 30
|
||
0680 0 0000 10 9990 CONST
|
||
0681 0 0000 12 9908 J
|
||
0682 0 0000 40 9908 J
|
||
0679 0683
|
||
0683 0 0000 13 9929 N
|
||
0684 0 1111 31 0686
|
||
0685 0 0000 36
|
||
0686 0 0001 33
|
||
|
||
0675 FOR J = (K+1,1,N)$
|
||
0687 0 0000 30
|
||
0686 0688
|
||
0685 0688
|
||
0688 0 0000 10 9914 K
|
||
0689 0 0000 14 9927
|
||
0690 0 0001 49 0010
|
||
0691 0 0000 12 9908 J
|
||
0692 0 0000 14 9926
|
||
0693 0 0001 49 0010
|
||
0693 0 0001 40 9911 TEMP
|
||
0694 0 0000 10 9914 K
|
||
0695 0 0000 14 9927
|
||
0696 0 0001 49 0010
|
||
0697 0 0000 12 9908 J
|
||
0698 0 0000 14 9926
|
||
0699 0 0001 49 0010
|
||
0699 0 0001 40 9910 TEMP
|
||
0700 0 0000 10 9990 CONST
|
||
0701 0 4400 28 9916
|
||
0702 1 0000 40 0000
|
||
0703 0 0000 10 9914 K
|
||
0704 0 0000 13 9990 CONST
|
||
0705 1 0000 40 9999
|
||
0706 0 0000 10 9914 K
|
||
0707 0 0000 14 9927
|
||
0708 0 0001 49 0010
|
||
0709 0 0000 14 9926
|
||
0710 0 0001 49 0010
|
||
0711 0 0000 12 0712
|
||
0712 0 0000 01 9928 A
|
||
0713 1 0000 40 9998
|
||
0714 0 0000 10 9926
|
||
0715 1 0000 40 9997
|
||
0716 0 0000 10 9908 J
|
||
0717 0 0000 14 9926
|
||
0718 0 0001 49 0010
|
||
0719 0 0000 12 0720
|
||
0720 0 0000 01 9928 A
|
||
0721 1 0000 40 9996
|
||
0722 0 0000 10 9927
|
||
0723 0 0000 14 9926
|
||
0724 0 0001 49 0010
|
||
|
||
0687 A(K,J) = A(K,J) - IP(1, K-1, A(K,), A(,J))$
|
||
|
||
0687
|
||
0725 0 0000 44 9916 IP
|
||
0726 0 0500 30 9916
|
||
0727 0 0000 42 9910 TEMP
|
||
0728 1 0000 23 9928 A
|
||
0729 0 0000 40 9910 TEMP
|
||
0730 0 0000 11 9910 TEMP
|
||
0731 0 0000 42 9911 TEMP
|
||
0732 1 0000 40 9928 A
|
||
0687 0734
|
||
0678 0733
|
||
0733 0 0000 30
|
||
0734 0 0000 10 9914 K
|
||
0735 0 0000 14 9924
|
||
0736 0 0001 49 0010
|
||
0736 0 0001 40 9911 TEMP
|
||
0737 0 0000 10 9914 K
|
||
0738 0 0000 14 9924
|
||
0739 0 0001 49 0010
|
||
0739 0 0001 40 9910 TEMP
|
||
0740 0 0000 10 9990 CONST
|
||
0741 0 4400 28 9916
|
||
0742 1 0000 40 0000
|
||
0743 0 0000 10 9914 K
|
||
0744 0 0000 13 9990 CONST
|
||
0745 1 0000 40 9999
|
||
0746 0 0000 10 9914 K
|
||
0747 0 0000 14 9927
|
||
0748 0 0001 49 0010
|
||
0749 0 0000 14 9926
|
||
0750 0 0001 49 0010
|
||
0751 0 0000 12 0752
|
||
0752 0 0000 01 9928 A
|
||
0753 1 0000 40 9998
|
||
0754 0 0000 10 9926
|
||
0755 1 0000 40 9997
|
||
0756 0 0000 10 0757
|
||
0757 0 0000 01 9925 B
|
||
0758 1 0000 40 9996
|
||
|
||
0734 B(K) = B(K) - IP(1, K-1, A(K,), B())
|
||
0759 0 0000 10 9924
|
||
0760 0 0000 44 9916 IP
|
||
0761 0 0500 30 9916
|
||
0762 0 0000 42 9910 TEMP
|
||
0763 1 0000 23 9925 B
|
||
0764 0 0000 40 9910 TEMP
|
||
0765 0 0000 11 9910 TEMP
|
||
0766 0 0000 42 9911 TEMP
|
||
0767 1 0000 40 9925 B
|
||
|
||
0759 END$
|
||
|
||
0759
|
||
0450 0769
|
||
0441 0768
|
||
0768 0 0000 30
|
||
0769 0 0000 10 9990 CONST
|
||
0770 0 0000 40 9912 I
|
||
0771 0 0000 44
|
||
0772 0 0000 30
|
||
0773 0 0000 10 9990 CONST
|
||
0774 0 0000 12 9912 I
|
||
0775 0 0000 40 9912 I
|
||
0772 0776
|
||
0776 0 0000 13 9929 N
|
||
0777 0 1111 31 0779
|
||
0778 0 0000 36
|
||
0779 0 0001 33
|
||
|
||
0769 FOR I = (1,1,N)$
|
||
0780 0 0000 30
|
||
0779 0781
|
||
0778 0781
|
||
0781 0 0000 10 9912 I
|
||
0782 0 0000 14 9922
|
||
0783 0 0001 49 0010
|
||
0783 0 0001 40 9911 TEMP
|
||
0784 0 0000 10 9912 I
|
||
0785 0 0000 14 9927
|
||
0786 0 0001 49 0010
|
||
0787 0 0000 12 9912 I
|
||
0788 0 0000 14 9926
|
||
0789 0 0001 49 0010
|
||
|
||
0780 Y(I) = A(I,I)$
|
||
0789 0 0001 40 9910 TEMP
|
||
0790 0 0000 42 9910 TEMP
|
||
0791 1 0000 10 9928 A
|
||
0792 0 0000 42 9911 TEMP
|
||
0793 1 0000 40 9923 Y
|
||
0780 0795
|
||
0771 0794
|
||
0794 0 0000 30
|
||
0795 0 0000 10 0796
|
||
0796 0 0000 01 9929 N
|
||
0797 0 4400 28 0201
|
||
0798 1 0000 40 0000
|
||
0799 0 0000 10 0800
|
||
0800 0 0000 01 9923 Y
|
||
0801 1 0000 40 9999
|
||
0802 0 0000 10 9922
|
||
0803 1 0000 40 9998
|
||
0804 0 0000 10 0805
|
||
0805 0 0000 01 9919 DET
|
||
0806 1 0000 40 9997
|
||
|
||
0795 PRODUCT ($ N, Y(), DET, EX7)$
|
||
0807 0 0000 10 0808
|
||
0808 0 0000 01 9918 EX7
|
||
0809 0 0000 44 0201 PRODU
|
||
0810 0 0400 30 0201
|
||
|
||
0811 DET = INT.DET$
|
||
|
||
0811
|
||
0811 0 0000 10 9915 INT
|
||
0812 0 0000 24 9919 DET
|
||
0813 0 0000 40 9919 DET
|
||
|
||
0814 COMMENT NOW FOR THE BACK SUBSTITUTION$
|
||
0814 0 0000 10 9929 N
|
||
0815 0 0000 40 9914 K
|
||
0816 0 0000 44
|
||
0817 0 0000 30
|
||
0818 0 0000 10 9914 K
|
||
0819 0 0000 13 9990 CONST
|
||
0820 0 0000 40 9914 K
|
||
0817 0821
|
||
0821 0 0000 13 9990 CONST
|
||
0822 0 1111 31 0824
|
||
0823 0 0000 36
|
||
0824 0 0000 33
|
||
|
||
0814 FOR K = (N,-1,1)$
|
||
0825 0 0000 30
|
||
0824 0826
|
||
0823 0826
|
||
|
||
0825 BEGIN
|
||
0826 0 0000 10 9914 K
|
||
0827 0 0000 14 9924
|
||
0828 0 0001 49 0010
|
||
0828 0 0001 40 9911 TEMP
|
||
0829 0 0000 10 9990 CONST
|
||
0830 0 0000 12 9914 K
|
||
0831 0 4400 28 9916
|
||
0832 1 0000 40 0000
|
||
0833 0 0000 10 9929 N
|
||
0834 1 0000 40 9999
|
||
0835 0 0000 10 9914 K
|
||
0836 0 0000 14 9927
|
||
0837 0 0001 49 0010
|
||
0838 0 0000 14 9926
|
||
0839 0 0001 49 0010
|
||
0840 0 0000 12 0841
|
||
0841 0 0000 01 9928 A
|
||
0842 1 0000 40 9998
|
||
0843 0 0000 10 9926
|
||
0844 1 0000 40 9997
|
||
0845 0 0000 10 0846
|
||
0846 0 0000 01 9923 Y
|
||
0847 1 0000 40 9996
|
||
|
||
0829 XX = B(K) - IP(K+1, N, A(K,), Y())$
|
||
0848 0 0000 10 9922
|
||
0849 0 0000 44 9916 IP
|
||
0850 0 0500 30 9916
|
||
0851 0 0000 42 9911 TEMP
|
||
0852 1 0000 23 9925 B
|
||
0853 0 0000 40 9911 TEMP
|
||
0854 0 0000 11 9911 TEMP
|
||
0855 0 0000 40 9907 XX
|
||
0856 0 0000 10 9914 K
|
||
0857 0 0000 14 9927
|
||
0858 0 0001 49 0010
|
||
0859 0 0000 12 9914 K
|
||
0860 0 0000 14 9926
|
||
0861 0 0001 49 0010
|
||
|
||
0856 XY = A(K,K)$
|
||
0861 0 0001 40 9911 TEMP
|
||
0862 0 0000 42 9911 TEMP
|
||
0863 1 0000 10 9928 A
|
||
0864 0 0000 40 9906 XY
|
||
|
||
0865 X = 1.0$
|
||
0865 0 0000 10 9987 CONST
|
||
0866 0 0000 40 9905 X
|
||
|
||
0867 X = X.X$
|
||
0867 0 0000 10 9905 X
|
||
0868 0 0000 24 9905 X
|
||
0869 0 0000 40 9905 X
|
||
0870 0 0000 10 9914 K
|
||
0871 0 0000 14 9922
|
||
0872 0 0001 49 0010
|
||
0872 0 0001 40 9911 TEMP
|
||
|
||
0870 Y(K) = XX/XY
|
||
0873 0 0000 10 9907 XX
|
||
0874 0 0002 45 0000
|
||
0875 0 0000 25 9906 XY
|
||
0876 0 0000 42 9911 TEMP
|
||
0877 1 0000 40 9923 Y
|
||
|
||
0878 END$
|
||
|
||
0878
|
||
0825 0879
|
||
0816 0878
|
||
0878 0 0000 30
|
||
0879 0 0000 42 0434
|
||
0880 1 0000 30 0000
|
||
|
||
0879 RETURN$
|
||
|
||
0881 END CROUT4()$
|
||
|
||
0881
|
||
0436 0881
|
||
0881 0 0410 40 0850
|
||
0882 0 0410 40 0849
|
||
0883 0 0410 40 0831
|
||
0884 0 0410 40 0761
|
||
0885 0 0410 40 0760
|
||
0886 0 0410 40 0741
|
||
0887 0 0410 40 0726
|
||
0888 0 0410 40 0725
|
||
0889 0 0410 40 0701
|
||
0890 0 0410 40 0502
|
||
0891 0 0410 40 0501
|
||
0892 0 0410 40 0477
|
||
0893 0 0000 10 9917
|
||
0894 0 0410 40 0626
|
||
0895 0 0000 10 9918
|
||
0896 0 0410 40 0808
|
||
0897 0 0410 40 0625
|
||
0898 0 0000 10 9919
|
||
0899 0 0410 40 0813
|
||
0900 0 0410 40 0812
|
||
0901 0 0410 40 0805
|
||
0902 0 0410 40 0624
|
||
0903 0 0000 10 9921
|
||
0904 0 0410 40 0538
|
||
0905 0 0000 10 9923
|
||
0906 0 0410 40 0877
|
||
0907 0 0410 40 0846
|
||
0908 0 0410 40 0800
|
||
0909 0 0410 40 0793
|
||
0910 0 0000 10 9925
|
||
0911 0 0410 40 0852
|
||
0912 0 0410 40 0767
|
||
0913 0 0410 40 0763
|
||
0914 0 0410 40 0757
|
||
0915 0 0410 40 0613
|
||
0916 0 0410 40 0607
|
||
0917 0 0410 40 0605
|
||
0918 0 0410 40 0596
|
||
0919 0 0000 10 9928
|
||
0920 0 0410 40 0863
|
||
0921 0 0410 40 0841
|
||
0922 0 0410 40 0791
|
||
0923 0 0410 40 0752
|
||
0924 0 0410 40 0732
|
||
0925 0 0410 40 0728
|
||
0926 0 0410 40 0720
|
||
0927 0 0410 40 0712
|
||
0928 0 0410 40 0673
|
||
0929 0 0410 40 0656
|
||
0930 0 0410 40 0647
|
||
0931 0 0410 40 0621
|
||
0932 0 0410 40 0590
|
||
0933 0 0410 40 0581
|
||
0934 0 0410 40 0579
|
||
0935 0 0410 40 0564
|
||
0936 0 0410 40 0528
|
||
0937 0 0410 40 0516
|
||
0938 0 0410 40 0508
|
||
0939 0 0410 40 0504
|
||
0940 0 0410 40 0496
|
||
0941 0 0410 40 0488
|
||
0942 0 0000 10 9929
|
||
0943 0 0410 40 0833
|
||
0944 0 0410 40 0814
|
||
0945 0 0410 40 0796
|
||
0946 0 0410 40 0776
|
||
0947 0 0410 40 0683
|
||
0948 0 0410 40 0635
|
||
0949 0 0410 40 0552
|
||
0950 0 0410 40 0459
|
||
0951 0 0410 40 0446
|
||
0952 0 0000 30 0437
|
||
0433 0953
|
||
0953 0 0000 30
|
||
0954 0 9904 01 0954
|
||
0955 0 0000 40 9894 TEMP
|
||
0956 0 0000 30
|
||
|
||
0953 PROCEDURE SOLV2 ($ N, B(,), C(), PIVOT(), Z()$ IP())$
|
||
|
||
0953 BEGIN COMMENT IT IS ASSUMED THAT A MATRIX A HAS ALREADY BEEN
|
||
|
||
0953 TRANSFORMED INTO B BY CROUT, BUT THAT A NEW COLUMN C HAS
|
||
|
||
0953 NOT BEEN PROCESSED. SOLV2() SOLVES THE SYSTEM BZ = C.
|
||
|
||
0953 AN INNERPRODUCT PROCEDURE MUST BE USED WITH SOLV2()$
|
||
|
||
0953 INTEGER K, N, PIVOT$
|
||
|
||
0953
|
||
0957 0 0000 10 9990 CONST
|
||
0958 0 0000 40 9893 K
|
||
0959 0 0000 44
|
||
0960 0 0000 30
|
||
0961 0 0000 10 9990 CONST
|
||
0962 0 0000 12 9893 K
|
||
0963 0 0000 40 9893 K
|
||
0960 0964
|
||
0964 0 0000 13 9904 N
|
||
0965 0 1111 31 0967
|
||
0966 0 0000 36
|
||
0967 0 0001 33
|
||
|
||
0957 FOR K = (1,1,N)$
|
||
0968 0 0000 30
|
||
0967 0969
|
||
0966 0969
|
||
|
||
0968 BEGIN
|
||
0969 0 0000 10 9893 K
|
||
0970 0 0000 14 9897
|
||
0971 0 0001 49 0010
|
||
0971 0 0001 40 9892 TEMP
|
||
0972 0 0000 42 9892 TEMP
|
||
0973 1 0000 10 9898 PIVOT
|
||
0974 0 0000 14 9899
|
||
0975 0 0001 49 0010
|
||
|
||
0972 TEMP = C(PIVOT(K))$
|
||
0975 0 0001 40 9892 TEMP
|
||
0976 0 0000 42 9892 TEMP
|
||
0977 1 0000 10 9900 C
|
||
0978 0 0000 40 9891 TEMP
|
||
0979 0 0000 10 9893 K
|
||
0980 0 0000 14 9897
|
||
0981 0 0001 49 0010
|
||
0981 0 0001 40 9892 TEMP
|
||
0982 0 0000 42 9892 TEMP
|
||
0983 1 0000 10 9898 PIVOT
|
||
0984 0 0000 14 9899
|
||
0985 0 0001 49 0010
|
||
0985 0 0001 40 9892 TEMP
|
||
0986 0 0000 10 9893 K
|
||
0987 0 0000 14 9899
|
||
0988 0 0001 49 0010
|
||
|
||
0979 C(PIVOT(K)) = C(K)$
|
||
0988 0 0001 40 9890 TEMP
|
||
0989 0 0000 42 9890 TEMP
|
||
0990 1 0000 10 9900 C
|
||
0991 0 0000 42 9892 TEMP
|
||
0992 1 0000 40 9900 C
|
||
0993 0 0000 10 9893 K
|
||
0994 0 0000 14 9899
|
||
0995 0 0001 49 0010
|
||
0995 0 0001 40 9892 TEMP
|
||
|
||
0993 C(K) = TEMP$
|
||
0996 0 0000 10 9891 TEMP
|
||
0997 0 0000 42 9892 TEMP
|
||
0998 1 0000 40 9900 C
|
||
0999 0 0000 10 9893 K
|
||
1000 0 0000 14 9899
|
||
1001 0 0001 49 0010
|
||
1001 0 0001 40 9892 TEMP
|
||
1002 0 0000 10 9893 K
|
||
1003 0 0000 14 9899
|
||
1004 0 0001 49 0010
|
||
1004 0 0001 40 9890 TEMP
|
||
1005 0 0000 10 9990 CONST
|
||
1006 0 4400 28 9894
|
||
1007 1 0000 40 0000
|
||
1008 0 0000 10 9893 K
|
||
1009 0 0000 13 9990 CONST
|
||
1010 1 0000 40 9999
|
||
1011 0 0000 10 9893 K
|
||
1012 0 0000 14 9902
|
||
1013 0 0001 49 0010
|
||
1014 0 0000 14 9901
|
||
1015 0 0001 49 0010
|
||
1016 0 0000 12 1017
|
||
1017 0 0000 01 9903 B
|
||
1018 1 0000 40 9998
|
||
1019 0 0000 10 9901
|
||
1020 1 0000 40 9997
|
||
1021 0 0000 10 1022
|
||
1022 0 0000 01 9900 C
|
||
1023 1 0000 40 9996
|
||
|
||
0999 C(K) = C(K) - IP(1, K-1, B(K,), C())
|
||
1024 0 0000 10 9899
|
||
1025 0 0000 44 9894 IP
|
||
1026 0 0500 30 9894
|
||
1027 0 0000 42 9890 TEMP
|
||
1028 1 0000 23 9900 C
|
||
1029 0 0000 40 9890 TEMP
|
||
1030 0 0000 11 9890 TEMP
|
||
1031 0 0000 42 9892 TEMP
|
||
1032 1 0000 40 9900 C
|
||
|
||
1024 END$
|
||
|
||
1024
|
||
0968 1034
|
||
0959 1033
|
||
1033 0 0000 30
|
||
1034 0 0000 10 9904 N
|
||
1035 0 0000 40 9893 K
|
||
1036 0 0000 44
|
||
1037 0 0000 30
|
||
1038 0 0000 10 9893 K
|
||
1039 0 0000 13 9990 CONST
|
||
1040 0 0000 40 9893 K
|
||
1037 1041
|
||
1041 0 0000 13 9990 CONST
|
||
1042 0 1111 31 1044
|
||
1043 0 0000 36
|
||
1044 0 0000 33
|
||
|
||
1034 FOR K = (N,-1,1)$
|
||
1045 0 0000 30
|
||
1044 1046
|
||
1043 1046
|
||
1046 0 0000 10 9893 K
|
||
1047 0 0000 14 9895
|
||
1048 0 0001 49 0010
|
||
1048 0 0001 40 9892 TEMP
|
||
1049 0 0000 10 9893 K
|
||
1050 0 0000 14 9899
|
||
1051 0 0001 49 0010
|
||
1051 0 0001 40 9890 TEMP
|
||
1052 0 0000 10 9990 CONST
|
||
1053 0 0000 12 9893 K
|
||
1054 0 4400 28 9894
|
||
1055 1 0000 40 0000
|
||
1056 0 0000 10 9904 N
|
||
1057 1 0000 40 9999
|
||
1058 0 0000 10 9893 K
|
||
1059 0 0000 14 9902
|
||
1060 0 0001 49 0010
|
||
1061 0 0000 14 9901
|
||
1062 0 0001 49 0010
|
||
1063 0 0000 12 1064
|
||
1064 0 0000 01 9903 B
|
||
1065 1 0000 40 9998
|
||
1066 0 0000 10 9901
|
||
1067 1 0000 40 9997
|
||
1068 0 0000 10 1069
|
||
1069 0 0000 01 9896 Z
|
||
1070 1 0000 40 9996
|
||
1071 0 0000 10 9895
|
||
1072 0 0000 44 9894 IP
|
||
1073 0 0500 30 9894
|
||
1074 0 0000 42 9890 TEMP
|
||
1075 1 0000 23 9900 C
|
||
1076 0 0000 40 9890 TEMP
|
||
1077 0 0000 10 9893 K
|
||
1078 0 0000 14 9902
|
||
1079 0 0001 49 0010
|
||
1080 0 0000 12 9893 K
|
||
1081 0 0000 14 9901
|
||
1082 0 0001 49 0010
|
||
|
||
1045 Z(K) = (C(K) - IP(K+1, N, B(K,), Z()))/B(K,K)$
|
||
1082 0 0001 40 9889 TEMP
|
||
1083 0 0000 11 9890 TEMP
|
||
1084 0 0002 45 0000
|
||
1085 0 0000 42 9889 TEMP
|
||
1086 1 0000 25 9903 B
|
||
1087 0 0000 42 9892 TEMP
|
||
1088 1 0000 40 9896 Z
|
||
1045 1090
|
||
1036 1089
|
||
1089 0 0000 30
|
||
|
||
1090 RETURN
|
||
1090 0 0000 42 0954
|
||
1091 1 0000 30 0000
|
||
|
||
1092 END SOLV2()$
|
||
|
||
1092
|
||
0956 1092
|
||
1092 0 0410 40 1073
|
||
1093 0 0410 40 1072
|
||
1094 0 0410 40 1054
|
||
1095 0 0410 40 1026
|
||
1096 0 0410 40 1025
|
||
1097 0 0410 40 1006
|
||
1098 0 0000 10 9896
|
||
1099 0 0410 40 1088
|
||
1100 0 0410 40 1069
|
||
1101 0 0000 10 9898
|
||
1102 0 0410 40 0983
|
||
1103 0 0410 40 0973
|
||
1104 0 0000 10 9900
|
||
1105 0 0410 40 1075
|
||
1106 0 0410 40 1032
|
||
1107 0 0410 40 1028
|
||
1108 0 0410 40 1022
|
||
1109 0 0410 40 0998
|
||
1110 0 0410 40 0992
|
||
1111 0 0410 40 0990
|
||
1112 0 0410 40 0977
|
||
1113 0 0000 10 9903
|
||
1114 0 0410 40 1086
|
||
1115 0 0410 40 1064
|
||
1116 0 0410 40 1017
|
||
1117 0 0000 10 9904
|
||
1118 0 0410 40 1056
|
||
1119 0 0410 40 1034
|
||
1120 0 0410 40 0964
|
||
1121 0 0000 30 0957
|
||
0953 1122
|
||
|
||
1122 COMMENT FORSYTHE TEST CROUT US169 EXT 2274$
|
||
1122 0 0000 30
|
||
1123 0 0006 60 0000
|
||
1124 0 0064 60 1908
|
||
1125 0 0006 60 0000
|
||
1126 1 0000 00 1124
|
||
1127 1 0000 00 1123
|
||
|
||
1122 FORMAT FRMTFL (W0, (6F19.8, W0))$
|
||
1122 1128
|
||
1128 0 0000 30
|
||
1129 0 0006 60 0000
|
||
1130 0 0064 90 1900
|
||
1131 0 0006 60 0000
|
||
1132 1 0000 00 1130
|
||
1133 1 0000 00 1129
|
||
|
||
1128 FORMAT FRMTFX (W0, (6I19, W0))$
|
||
1128 1134
|
||
|
||
1134 INTEGER PIVOT()$
|
||
|
||
1134 INTEGER EX$
|
||
|
||
1134 INTEGER I, J, N$
|
||
|
||
1134 ARRAY A(70,70), B(70), Y(70), C(70), PIVOT(70)$
|
||
1134 0 0000 30
|
||
1135 0 0000 30
|
||
1136 0 0100 42 1136
|
||
1137 0 0002 20 1135
|
||
1138 0 0000 40 4708 N
|
||
1139 0 0000 10 9990 CONST
|
||
1140 0 0000 40 4707 I
|
||
1141 0 0000 44
|
||
1142 0 0000 30
|
||
1143 0 0000 10 9990 CONST
|
||
1144 0 0000 12 4707 I
|
||
1145 0 0000 40 4707 I
|
||
1142 1146
|
||
1146 0 0000 13 4708 N
|
||
1147 0 1111 31 1149
|
||
1148 0 0000 36
|
||
1149 0 0001 33
|
||
1150 0 0000 30
|
||
1149 1151
|
||
1148 1151
|
||
1151 0 0000 10 9990 CONST
|
||
1152 0 0000 40 4706 J
|
||
1153 0 0000 44
|
||
1154 0 0000 30
|
||
1155 0 0000 10 9990 CONST
|
||
1156 0 0000 12 4706 J
|
||
1157 0 0000 40 4706 J
|
||
1154 1158
|
||
1158 0 0000 13 4708 N
|
||
1159 0 1111 31 1161
|
||
1160 0 0000 36
|
||
1161 0 0001 33
|
||
1162 0 0000 30
|
||
1161 1163
|
||
1160 1163
|
||
1163 0 0000 10 4707 I
|
||
1164 0 0000 14 4705 CONST
|
||
1165 0 0001 49 0010
|
||
1166 0 0000 12 4706 J
|
||
1167 0 0000 40 4704 TEMP
|
||
1168 0 0000 42 1168
|
||
1169 0 0002 20 1135
|
||
1170 0 0000 42 4704 TEMP
|
||
1171 1 0000 40 4918 A
|
||
1162 1173
|
||
1153 1172
|
||
1172 0 0000 30
|
||
1173 0 0000 42 1173
|
||
1174 0 0002 20 1135
|
||
1175 0 0000 42 4707 I
|
||
1176 1 0000 40 4918 B
|
||
1150 1178
|
||
1141 1177
|
||
1177 0 0000 30
|
||
|
||
1134 INPUT DATA (N, FOR I=(1,1,N)$ (FOR J=(1,1,N)$ A(I,J), B(I)))$
|
||
1178 0 0009 43 0000
|
||
1179 0 0000 30 1135
|
||
1134 1180
|
||
1180 0 0000 30
|
||
1181 0 0000 30
|
||
1182 0 0100 42 1182
|
||
1183 0 0002 20 1181
|
||
1184 0 0000 40 4708 N
|
||
1185 0 0000 10 9990 CONST
|
||
1186 0 0000 40 4707 I
|
||
1187 0 0000 44
|
||
1188 0 0000 30
|
||
1189 0 0000 10 9990 CONST
|
||
1190 0 0000 12 4707 I
|
||
1191 0 0000 40 4707 I
|
||
1188 1192
|
||
1192 0 0000 13 4708 N
|
||
1193 0 1111 31 1195
|
||
1194 0 0000 36
|
||
1195 0 0001 33
|
||
1196 0 0000 30
|
||
1195 1197
|
||
1194 1197
|
||
1197 0 0000 42 1197
|
||
1198 0 0002 20 1181
|
||
1199 0 0000 42 4707 I
|
||
1200 1 0000 40 4778 C
|
||
1196 1202
|
||
1187 1201
|
||
1201 0 0000 30
|
||
|
||
1180 INPUT VECTOR (N, FOR I=(1,1,N)$ C(I))$
|
||
|
||
1180
|
||
1202 0 0009 43 0000
|
||
1203 0 0000 30 1181
|
||
1180 1204
|
||
|
||
1204 START..
|
||
|
||
1204 READ ($$ DATA)$
|
||
1204 0 0000 10 1205
|
||
1205 0 0000 01 1135 DATA
|
||
1206 0 0000 44 4565 READ
|
||
1207 0 0000 30 4565
|
||
|
||
1208 READ ($$ VECTOR)$
|
||
1208 0 0000 10 1209
|
||
1209 0 0000 01 1181 VECTO
|
||
1210 0 0000 44 4565 READ
|
||
1211 0 0000 30 4565
|
||
1212 0 0000 30
|
||
1213 0 0000 30
|
||
|
||
1212 OUTPUT ORDER (N)$
|
||
1214 0 0000 10 4708 N
|
||
1215 0 0100 42 1215
|
||
1216 0 0002 20 1213
|
||
1217 0 0009 43 0000
|
||
1218 0 0000 30 1213
|
||
1212 1219
|
||
1219 0 0000 30
|
||
1220 0 0000 30
|
||
1221 0 0000 10 9990 CONST
|
||
1222 0 0000 40 4707 I
|
||
1223 0 0000 44
|
||
1224 0 0000 30
|
||
1225 0 0000 10 9990 CONST
|
||
1226 0 0000 12 4707 I
|
||
1227 0 0000 40 4707 I
|
||
1224 1228
|
||
1228 0 0000 13 4708 N
|
||
1229 0 1111 31 1231
|
||
1230 0 0000 36
|
||
1231 0 0001 33
|
||
1232 0 0000 30
|
||
1231 1233
|
||
1230 1233
|
||
1233 0 0000 10 9990 CONST
|
||
1234 0 0000 40 4706 J
|
||
1235 0 0000 44
|
||
1236 0 0000 30
|
||
1237 0 0000 10 9990 CONST
|
||
1238 0 0000 12 4706 J
|
||
1239 0 0000 40 4706 J
|
||
1236 1240
|
||
1240 0 0000 13 4708 N
|
||
1241 0 1111 31 1243
|
||
1242 0 0000 36
|
||
1243 0 0001 33
|
||
1244 0 0000 30
|
||
1243 1245
|
||
1242 1245
|
||
1245 0 0000 10 4707 I
|
||
1246 0 0000 14 4705 CONST
|
||
1247 0 0001 49 0010
|
||
1248 0 0000 12 4706 J
|
||
1249 0 0000 40 4704 TEMP
|
||
1250 0 0000 42 4704 TEMP
|
||
1251 1 0000 10 4918 A
|
||
1252 0 0000 42 1252
|
||
1253 0 0002 20 1220
|
||
1244 1255
|
||
1235 1254
|
||
1254 0 0000 30
|
||
1255 0 0000 42 4707 I
|
||
1256 1 0000 10 4918 B
|
||
1257 0 0000 42 1257
|
||
1258 0 0002 20 1220
|
||
1232 1260
|
||
1223 1259
|
||
1259 0 0000 30
|
||
|
||
1219 OUTPUT DATAO (FOR I=(1,1,N)$ (FOR J=(1,1,N)$ A(I,J), B(I)))$
|
||
1260 0 0009 43 0000
|
||
1261 0 0000 30 1220
|
||
1219 1262
|
||
1262 0 0000 30
|
||
1263 0 0000 30
|
||
1264 0 0000 10 9990 CONST
|
||
1265 0 0000 40 4707 I
|
||
1266 0 0000 44
|
||
1267 0 0000 30
|
||
1268 0 0000 10 9990 CONST
|
||
1269 0 0000 12 4707 I
|
||
1270 0 0000 40 4707 I
|
||
1267 1271
|
||
1271 0 0000 13 4708 N
|
||
1272 0 1111 31 1274
|
||
1273 0 0000 36
|
||
1274 0 0001 33
|
||
1275 0 0000 30
|
||
1274 1276
|
||
1273 1276
|
||
1276 0 0000 42 4707 I
|
||
1277 1 0000 10 4778 C
|
||
1278 0 0000 42 1278
|
||
1279 0 0002 20 1263
|
||
1275 1281
|
||
1266 1280
|
||
1280 0 0000 30
|
||
|
||
1262 OUTPUT VECTORO (FOR I=(1,1,N)$ C(I))$
|
||
1281 0 0009 43 0000
|
||
1282 0 0000 30 1263
|
||
1262 1283
|
||
1283 0 0000 10 1284
|
||
1284 0 0000 01 1213 ORDER
|
||
1285 0 4400 28 4325
|
||
1286 1 0000 40 0000
|
||
|
||
1283 WRITE ($$ ORDER, FRMTFX)$
|
||
1287 0 0000 10 1288
|
||
1288 0 0000 01 1129 FRMTF
|
||
1289 0 0000 44 4325 WRITE
|
||
1290 0 0100 30 4325
|
||
1291 0 0000 10 1292
|
||
1292 0 0000 01 1220 DATAO
|
||
1293 0 4400 28 4325
|
||
1294 1 0000 40 0000
|
||
|
||
1291 WRITE ($$ DATAO, FRMTFL)$
|
||
1295 0 0000 10 1296
|
||
1296 0 0000 01 1123 FRMTF
|
||
1297 0 0000 44 4325 WRITE
|
||
1298 0 0100 30 4325
|
||
1299 0 0000 10 1300
|
||
1300 0 0000 01 1213 ORDER
|
||
1301 0 4400 28 4325
|
||
1302 1 0000 40 0000
|
||
|
||
1299 WRITE ($$ ORDER, FRMTFX)$
|
||
1303 0 0000 10 1304
|
||
1304 0 0000 01 1129 FRMTF
|
||
1305 0 0000 44 4325 WRITE
|
||
1306 0 0100 30 4325
|
||
1307 0 0000 10 1308
|
||
1308 0 0000 01 1263 VECTO
|
||
1309 0 4400 28 4325
|
||
1310 1 0000 40 0000
|
||
|
||
1307 WRITE ($$ VECTORO, FRMTFL)$
|
||
1311 0 0000 10 1312
|
||
1312 0 0000 01 1123 FRMTF
|
||
1313 0 0000 44 4325 WRITE
|
||
1314 0 0100 30 4325
|
||
1315 0 0000 10 1316
|
||
1316 0 0000 01 4708 N
|
||
1317 0 4400 28 0434
|
||
1318 1 0000 40 0000
|
||
1319 0 0000 10 1320
|
||
1320 0 0000 01 4918 A
|
||
1321 1 0000 40 9999
|
||
1322 0 0000 10 4705 CONST
|
||
1323 1 0000 40 9998
|
||
1324 0 0000 10 9990 CONST
|
||
1325 1 0000 40 9997
|
||
1326 0 0000 10 1327
|
||
1327 0 0000 01 4918 B
|
||
1328 1 0000 40 9996
|
||
1329 0 0000 10 9990 CONST
|
||
1330 1 0000 40 9995
|
||
1331 0 0000 10 1332
|
||
1332 0 0000 01 4848 Y
|
||
1333 1 0000 40 9994
|
||
1334 0 0000 10 9990 CONST
|
||
1335 1 0000 40 9993
|
||
1336 0 0000 10 1337
|
||
1337 0 0000 01 4708 PIVOT
|
||
1338 1 0000 40 9992
|
||
1339 0 0000 10 9990 CONST
|
||
1340 1 0000 40 9991
|
||
1341 0 0000 10 1342
|
||
1342 0 0000 01 4324 DET
|
||
1343 1 0000 40 9990
|
||
|
||
1315 CROUT4 ($ N, A(,), B(), Y(), PIVOT(), DET, EX$
|
||
1344 0 0000 10 1345
|
||
1345 0 0000 01 4323 EX
|
||
1346 1 0000 40 9989
|
||
1347 0 0000 10 1348
|
||
1348 0 0000 01 0000 SINGU
|
||
1349 1 0000 40 9988
|
||
|
||
1347 SINGULAR, INNERPRODUCT())$
|
||
1350 0 0000 10 1351
|
||
1351 0 0000 01 0396
|
||
1352 0 0000 44 0434 CROUT
|
||
1353 0 1300 30 0434
|
||
1354 0 0000 10 1355
|
||
1355 0 0000 01 1220 DATAO
|
||
1356 0 4400 28 4325
|
||
1357 1 0000 40 0000
|
||
|
||
1354 WRITE ($$ DATAO, FRMTFL)$
|
||
1358 0 0000 10 1359
|
||
1359 0 0000 01 1123 FRMTF
|
||
1360 0 0000 44 4325 WRITE
|
||
1361 0 0100 30 4325
|
||
1362 0 0000 30
|
||
1363 0 0000 30
|
||
1364 0 0000 10 9990 CONST
|
||
1365 0 0000 40 4707 I
|
||
1366 0 0000 44
|
||
1367 0 0000 30
|
||
1368 0 0000 10 9990 CONST
|
||
1369 0 0000 12 4707 I
|
||
1370 0 0000 40 4707 I
|
||
1367 1371
|
||
1371 0 0000 13 4708 N
|
||
1372 0 1111 31 1374
|
||
1373 0 0000 36
|
||
1374 0 0001 33
|
||
1375 0 0000 30
|
||
1374 1376
|
||
1373 1376
|
||
1376 0 0000 42 4707 I
|
||
1377 1 0000 10 4848 Y
|
||
1378 0 0000 42 1378
|
||
1379 0 0002 20 1363
|
||
1375 1381
|
||
1366 1380
|
||
1380 0 0000 30
|
||
|
||
1362 OUTPUT ANSWER (FOR I=(1,1,N)$ Y(I))$
|
||
1381 0 0009 43 0000
|
||
1382 0 0000 30 1363
|
||
1362 1383
|
||
1383 0 0000 30
|
||
1384 0 0000 30
|
||
1385 0 0000 10 4708 N
|
||
1386 0 0100 42 1386
|
||
1387 0 0002 20 1384
|
||
1388 0 0000 10 9990 CONST
|
||
1389 0 0000 40 4707 I
|
||
1390 0 0000 44
|
||
1391 0 0000 30
|
||
1392 0 0000 10 9990 CONST
|
||
1393 0 0000 12 4707 I
|
||
1394 0 0000 40 4707 I
|
||
1391 1395
|
||
1395 0 0000 13 4708 N
|
||
1396 0 1111 31 1398
|
||
1397 0 0000 36
|
||
1398 0 0001 33
|
||
1399 0 0000 30
|
||
1398 1400
|
||
1397 1400
|
||
1400 0 0000 42 4707 I
|
||
1401 1 0000 10 4708 PIVOT
|
||
1402 0 0100 42 1402
|
||
1403 0 0002 20 1384
|
||
1399 1405
|
||
1390 1404
|
||
1404 0 0000 30
|
||
|
||
1383 OUTPUT PIVOTO (N, FOR I=(1,1,N)$ PIVOT(I))$
|
||
1405 0 0009 43 0000
|
||
1406 0 0000 30 1384
|
||
1383 1407
|
||
1407 0 0000 30
|
||
1408 0 0000 30
|
||
|
||
1407 OUTPUT DETO (DET)$
|
||
1409 0 0000 10 4324 DET
|
||
1410 0 0000 42 1410
|
||
1411 0 0002 20 1408
|
||
1412 0 0009 43 0000
|
||
1413 0 0000 30 1408
|
||
1407 1414
|
||
1414 0 0000 30
|
||
1415 0 0000 30
|
||
|
||
1414 OUTPUT EXPO (EX)$
|
||
1416 0 0000 10 4323 EX
|
||
1417 0 0100 42 1417
|
||
1418 0 0002 20 1415
|
||
1419 0 0009 43 0000
|
||
1420 0 0000 30 1415
|
||
1414 1421
|
||
1421 0 0000 10 1422
|
||
1422 0 0000 01 1384 PIVOT
|
||
1423 0 4400 28 4325
|
||
1424 1 0000 40 0000
|
||
|
||
1421 WRITE ($$ PIVOTO, FRMTFX)$
|
||
1425 0 0000 10 1426
|
||
1426 0 0000 01 1129 FRMTF
|
||
1427 0 0000 44 4325 WRITE
|
||
1428 0 0100 30 4325
|
||
1429 0 0000 10 1430
|
||
1430 0 0000 01 1363 ANSWE
|
||
1431 0 4400 28 4325
|
||
1432 1 0000 40 0000
|
||
|
||
1429 WRITE ($$ ANSWER, FRMTFL)$
|
||
1433 0 0000 10 1434
|
||
1434 0 0000 01 1123 FRMTF
|
||
1435 0 0000 44 4325 WRITE
|
||
1436 0 0100 30 4325
|
||
1437 0 0000 10 1438
|
||
1438 0 0000 01 1408 DETO
|
||
1439 0 4400 28 4325
|
||
1440 1 0000 40 0000
|
||
|
||
1437 WRITE ($$ DETO, FRMTFL)$
|
||
1441 0 0000 10 1442
|
||
1442 0 0000 01 1123 FRMTF
|
||
1443 0 0000 44 4325 WRITE
|
||
1444 0 0100 30 4325
|
||
1445 0 0000 10 1446
|
||
1446 0 0000 01 1415 EXPO
|
||
1447 0 4400 28 4325
|
||
1448 1 0000 40 0000
|
||
|
||
1445 WRITE ($$ EXPO, FRMTFX)$
|
||
1449 0 0000 10 1450
|
||
1450 0 0000 01 1129 FRMTF
|
||
1451 0 0000 44 4325 WRITE
|
||
1452 0 0100 30 4325
|
||
1453 0 0000 10 1454
|
||
1454 0 0000 01 4708 N
|
||
1455 0 4400 28 0954
|
||
1456 1 0000 40 0000
|
||
1457 0 0000 10 1458
|
||
1458 0 0000 01 4918 A
|
||
1459 1 0000 40 9999
|
||
1460 0 0000 10 4705 CONST
|
||
1461 1 0000 40 9998
|
||
1462 0 0000 10 9990 CONST
|
||
1463 1 0000 40 9997
|
||
1464 0 0000 10 1465
|
||
1465 0 0000 01 4778 C
|
||
1466 1 0000 40 9996
|
||
1467 0 0000 10 9990 CONST
|
||
1468 1 0000 40 9995
|
||
1469 0 0000 10 1470
|
||
1470 0 0000 01 4708 PIVOT
|
||
1471 1 0000 40 9994
|
||
1472 0 0000 10 9990 CONST
|
||
1473 1 0000 40 9993
|
||
1474 0 0000 10 1475
|
||
1475 0 0000 01 4848 Y
|
||
1476 1 0000 40 9992
|
||
1477 0 0000 10 9990 CONST
|
||
1478 1 0000 40 9991
|
||
|
||
1453 SOLV2 ($ N, A(,), C(), PIVOT(), Y()$ INNERPRODUCT())$
|
||
1479 0 0000 10 1480
|
||
1480 0 0000 01 0396
|
||
1481 0 0000 44 0954 SOLV2
|
||
1482 0 1000 30 0954
|
||
1483 0 0000 10 1484
|
||
1484 0 0000 01 1263 VECTO
|
||
1485 0 4400 28 4325
|
||
1486 1 0000 40 0000
|
||
|
||
1483 WRITE ($$ VECTORO, FRMTFL)$
|
||
1487 0 0000 10 1488
|
||
1488 0 0000 01 1123 FRMTF
|
||
1489 0 0000 44 4325 WRITE
|
||
1490 0 0100 30 4325
|
||
1491 0 0000 10 1492
|
||
1492 0 0000 01 1363 ANSWE
|
||
1493 0 4400 28 4325
|
||
1494 1 0000 40 0000
|
||
|
||
1491 WRITE ($$ ANSWER, FRMTFL)$
|
||
1495 0 0000 10 1496
|
||
1496 0 0000 01 1123 FRMTF
|
||
1497 0 0000 44 4325 WRITE
|
||
1498 0 0100 30 4325
|
||
|
||
1499 GO TO START$
|
||
|
||
1499
|
||
1499 0 0000 30 1204 START
|
||
|
||
1500 SINGULAR..
|
||
1348 1500
|
||
|
||
1500 WRITE ($$ FRMTSI)$
|
||
1500 0 0000 10 1501
|
||
1501 0 0000 01 0000 FRMTS
|
||
1502 0 0000 44 4325 WRITE
|
||
1503 0 0000 30 4325
|
||
1504 0 0000 30
|
||
1501 1505
|
||
1505 0 0006 60 0000
|
||
1506 2 6249 55 4764
|
||
1507 3 5341 59 1400
|
||
1508 0 0006 60 0000
|
||
1509 1 0000 00 1505
|
||
|
||
1504 FORMAT FRMTSI (W0, *SINGULAR*, W0)$
|
||
1504 1510
|
||
|
||
1510 GO TO START$
|
||
1510 0 0000 30 1204 START
|
||
|
||
1511 FINISH$
|
||
1511 0 9669 00 9669
|
||
1512 0 1000 60 0000
|
||
9990 0 0000 00 0001 POOL
|
||
9992 0 0000 00 0010
|
||
9986 0 0000 00 0020
|
||
9944 0 0000 00 0011
|
||
9943 0 0000 00 0041
|
||
9942 0 0000 00 0040
|
||
9941 0 0000 00 0058
|
||
9940 0 0000 00 0009
|
||
4705 0 0000 00 0070
|
||
9987 0 5110 00 0000
|
||
9994 0 4110 00 0000
|
||
9985 0 5210 00 0000
|
||
4325 0 0100 01 0000
|
||
4326 0 0410 40 4361
|
||
4327 0 0000 42 4329
|
||
4328 0 0010 29 4561
|
||
4329 0 0300 29 0101
|
||
4330 0 0000 42 4325
|
||
4331 0 0412 40 4410
|
||
4332 1 0000 41 9999
|
||
4333 0 0009 43 7557
|
||
4334 0 2200 37 4338
|
||
4335 0 0000 42 0100
|
||
4336 1 0000 44 0000
|
||
4337 1 2201 37 0001
|
||
4338 0 0000 40 0132
|
||
4339 0 0412 40 4392
|
||
4340 0 0000 41 0130
|
||
4341 0 0000 37 4359
|
||
4342 0 0000 10 0134
|
||
4343 0 0000 46 0133
|
||
4344 0 5246 37 4521
|
||
4345 0 5267 37 4497
|
||
4346 0 5249 37 4435
|
||
4347 0 5262 37 4548
|
||
4348 0 5241 37 4350
|
||
4349 0 0000 30 4536
|
||
4350 0 0000 41 4368
|
||
4351 0 0001 27 0129
|
||
4352 0 0000 44 4493
|
||
4353 0 0000 30 4484
|
||
4354 0 0000 10 0129
|
||
4355 0 0000 36 4375
|
||
4356 0 2299 37 4389
|
||
4357 0 0001 49 0010
|
||
4358 0 0000 30 4351
|
||
4359 0 0000 42 4326
|
||
4360 0 0000 31 4277 ERROR
|
||
4361 0 0000 10 9999
|
||
4362 0 0401 26 4361
|
||
4363 0 0000 33 4377
|
||
4364 0 0000 41 4368
|
||
4365 0 0001 33 4393
|
||
4366 0 2214 36 4361
|
||
4367 0 0000 44 4493
|
||
4368 0 9900 30 4484
|
||
4369 0 2299 37 4361
|
||
4370 0 0001 49 0010
|
||
4371 0 0000 30 4366
|
||
4372 0 0000 10 0129
|
||
4373 0 0000 19 0125
|
||
4374 0 0000 19 0125
|
||
4375 0 0000 10 0130
|
||
4376 0 3300 36 4359
|
||
4377 0 3300 36 4379
|
||
4378 0 0000 13 4562
|
||
4379 0 0000 40 0130
|
||
4380 0 0210 40 0131
|
||
4381 0 0001 48 0002
|
||
4382 0 0310 40 0129
|
||
4383 0 7242 36 4372
|
||
4384 0 7266 36 4413
|
||
4385 0 7257 36 4421
|
||
4386 0 7263 36 4423
|
||
4387 0 7243 36 4412
|
||
4388 0 0300 36 4361
|
||
4389 0 0000 10 0132
|
||
4390 0 0000 40 0134
|
||
4391 0 0009 33 4372
|
||
4392 0 0000 30 9999
|
||
4393 0 0000 42 4361
|
||
4394 0 6600 36 4409
|
||
4395 0 0001 49 0003
|
||
4396 0 0300 37 4403
|
||
4397 0 0000 13 4562
|
||
4398 0 3300 36 4405
|
||
4399 1 6301 27 9999
|
||
4400 1 0000 41 9999
|
||
4401 0 0411 40 4361
|
||
4402 0 0000 30 4361
|
||
4403 1 3310 40 9999
|
||
4404 0 0000 30 4397
|
||
4405 1 0000 10 9999
|
||
4406 0 0001 48 0003
|
||
4407 1 6610 40 9999
|
||
4408 0 0000 30 4361
|
||
4409 0 0000 10 0132
|
||
4410 0 0009 33 9999
|
||
4411 0 0000 30 4400
|
||
4412 0 1008 61 0124
|
||
4413 0 0001 48 0003
|
||
4414 0 3111 40 4415
|
||
4415 0 2408 44 4310 REED
|
||
4416 0 0101 30 4310 REED
|
||
4417 0 0000 42 4419
|
||
4418 0 0010 29 4561
|
||
4419 0 0240 29 0101
|
||
4420 0 0000 30 4375
|
||
4421 0 1008 61 0124
|
||
4422 0 0000 30 4417
|
||
4423 0 0300 36 4428
|
||
4424 0 0010 09 4427
|
||
4425 0 0001 12 4405
|
||
4426 0 0000 30 4423
|
||
4427 2 1602 02 0202
|
||
4428 0 0000 10 0125
|
||
4429 0 0000 36 4375
|
||
4430 0 0000 12 4498
|
||
4431 0 0001 48 0004
|
||
4432 0 3211 40 4433
|
||
4433 0 0990 09 0101
|
||
4434 0 0000 30 4417
|
||
4435 0 0000 42 4445
|
||
4436 0 0000 36 4494
|
||
4437 0 1100 36 4495
|
||
4438 0 0000 40 0134
|
||
4439 0 0003 45 0080
|
||
4440 0 0002 40 0128
|
||
4441 0 0000 46 0126
|
||
4442 0 0410 40 0126
|
||
4443 0 0001 40 0127
|
||
4444 0 0000 10 0134
|
||
4445 0 0002 48 0010
|
||
4446 0 0001 43 0023
|
||
4447 0 0000 12 0129
|
||
4448 0 0000 13 0128
|
||
4449 0 0001 33 4536
|
||
4450 0 0000 19 0125
|
||
4451 0 5000 19 0125
|
||
4452 0 0000 10 0134
|
||
4453 0 0000 41 4489
|
||
4454 0 0000 44 4493
|
||
4455 0 0001 33 4483
|
||
4456 0 0001 27 0128
|
||
4457 0 0000 32 4469
|
||
4458 0 0000 41 0127
|
||
4459 0 0000 37 4375
|
||
4460 0 0000 41 4446
|
||
4461 0 0000 44 4493
|
||
4462 0 0000 30 4483
|
||
4463 0 0000 10 0127
|
||
4464 0 0000 13 4451
|
||
4465 0 0000 42 4381
|
||
4466 0 0000 46 0129
|
||
4467 0 0003 26 0129
|
||
4468 0 0000 30 4438
|
||
4469 0 0001 27 0126
|
||
4470 0 0000 41 0126
|
||
4471 0 0000 37 4477
|
||
4472 0 0000 41 0133
|
||
4473 0 0000 37 4479
|
||
4474 0 0001 27 0133
|
||
4475 0 0000 41 4439
|
||
4476 0 0000 30 4483
|
||
4477 0 0000 41 4406
|
||
4478 0 0000 30 4483
|
||
4479 0 0000 10 0134
|
||
4480 0 0000 41 4498
|
||
4481 0 0001 49 0001
|
||
4482 0 0000 40 0134
|
||
4483 0 0001 49 0018
|
||
4484 0 0001 48 0008
|
||
4485 0 0001 42 0125
|
||
4486 1 0000 49 0008
|
||
4487 0 9400 28 0125
|
||
4488 0 0024 21 4493
|
||
4489 0 0000 43 0020
|
||
4490 1 0000 12 0125
|
||
4491 1 0010 40 0125
|
||
4492 0 0002 26 0125
|
||
4493 0 0000 30 9999
|
||
4494 0 0009 21 4439
|
||
4495 0 0000 49 0051
|
||
4496 0 0001 21 4437
|
||
4497 0 0002 45 7557
|
||
4498 0 0001 48 0008
|
||
4499 0 0410 18 4563
|
||
4500 0 0001 40 0134
|
||
4501 0 0001 34 4511
|
||
4502 0 0410 40 4509
|
||
4503 0 0449 27 4509
|
||
4504 0 0000 42 4509
|
||
4505 0 0000 10 0130
|
||
4506 0 0001 48 0004
|
||
4507 0 4211 40 4510
|
||
4508 0 0000 10 4509
|
||
4509 0 0002 45 9999
|
||
4510 0 0000 20 4440
|
||
4511 0 0000 43 7557
|
||
4512 0 0000 13 4563
|
||
4513 0 0210 40 0133
|
||
4514 0 0211 18 0130
|
||
4515 0 0000 35 4518
|
||
4516 0 0000 10 4563
|
||
4517 0 0000 30 4502
|
||
4518 0 0000 10 4495
|
||
4519 0 1110 40 0134
|
||
4520 0 0000 30 4502
|
||
4521 0 0304 27 0129
|
||
4522 0 0000 32 4525
|
||
4523 0 0000 30 4535
|
||
4524 0 0000 10 4564
|
||
4525 0 0000 36 4524
|
||
4526 0 0002 45 7557
|
||
4527 0 2210 40 4450
|
||
4528 0 0001 48 0008
|
||
4529 0 0001 40 0134
|
||
4530 0 0002 48 0010
|
||
4531 0 0000 12 0130
|
||
4532 0 0000 17 4495
|
||
4533 0 0310 18 0129
|
||
4534 0 0001 34 4544
|
||
4535 0 0004 26 0129
|
||
4536 0 0000 10 0129
|
||
4537 0 0000 19 0125
|
||
4538 0 0000 19 0125
|
||
4539 0 0002 27 0125
|
||
4540 0 0000 41 4551
|
||
4541 0 0000 44 4493
|
||
4542 0 0000 30 4483
|
||
4543 0 0000 30 4375
|
||
4544 0 0000 41 4450
|
||
4545 0 0000 10 4481
|
||
4546 0 0000 42 0131
|
||
4547 0 0001 20 4440
|
||
4548 0 0002 45 7557
|
||
4549 0 0001 48 0008
|
||
4550 0 0001 40 0134
|
||
4551 0 0000 43 0014
|
||
4552 0 0000 13 4563
|
||
4553 0 0002 45 7557
|
||
4554 0 0001 33 4559
|
||
4555 0 0210 18 0130
|
||
4556 0 0000 34 4536
|
||
4557 0 0000 12 4481
|
||
4558 0 0000 30 4546
|
||
4559 0 0010 40 0133
|
||
4560 0 0000 30 4545
|
||
4561 2 0000 00 0000
|
||
4562 0 0010 00 0000
|
||
4563 0 0000 00 0050
|
||
4564 0 5000 00 0000
|
||
4565 0 4693 01 0000
|
||
4566 0 0410 40 4655
|
||
4567 0 0410 40 4656
|
||
4568 0 0401 26 4656
|
||
4569 0 0000 44 4662
|
||
4570 0 0000 30 4655
|
||
4571 0 0000 42 4693
|
||
4572 1 0000 46 0000
|
||
4573 0 0000 46 4701
|
||
4574 0 0000 44 4253 RITE
|
||
4575 0 0117 30 4253 RITE
|
||
4576 0 0000 46 4703
|
||
4577 0 0502 26 4703
|
||
4578 0 0000 42 4693
|
||
4579 0 9999 20 4583
|
||
4580 0 6400 28 4701
|
||
4581 0 0000 10 0101
|
||
4582 0 9999 20 4677
|
||
4583 0 9400 28 4703
|
||
4584 1 0000 10 0101
|
||
4585 0 0000 42 4703
|
||
4586 1 0000 49 0000
|
||
4587 0 0502 26 4703
|
||
4588 0 0000 41 4701
|
||
4589 0 0000 37 4603
|
||
4590 0 0160 21 4574
|
||
4591 0 4400 28 4701
|
||
4592 0 2213 36 4618
|
||
4593 0 0000 48 0008
|
||
4594 1 0000 49 0000
|
||
4595 0 0000 19 4700
|
||
4596 0 4102 27 4701
|
||
4597 0 0000 32 4583
|
||
4598 0 0000 10 4700
|
||
4599 0 0000 44 4662
|
||
4600 0 0000 30 4655
|
||
4601 0 5212 27 4701
|
||
4602 0 0000 30 4583
|
||
4603 0 0162 21 4574
|
||
4604 0 1108 36 4634
|
||
4605 0 2203 36 4630
|
||
4606 0 2220 36 4628
|
||
4607 0 2234 36 4628
|
||
4608 0 2223 36 4623
|
||
4609 0 2213 36 4685
|
||
4610 0 2214 36 4616
|
||
4611 0 0000 42 4697
|
||
4612 0 9999 20 4583
|
||
4613 0 0000 44 4662
|
||
4614 0 0000 30 4641
|
||
4615 0 0000 30 4583
|
||
4616 0 9416 26 4703
|
||
4617 0 0000 30 4611
|
||
4618 0 0000 46 4701
|
||
4619 0 0000 10 4700
|
||
4620 0 0000 44 4662
|
||
4621 0 9992 20 4655
|
||
4622 0 0000 30 4583
|
||
4623 0 0001 26 4699
|
||
4624 0 0000 44 4662
|
||
4625 0 0000 30 4641
|
||
4626 0 0001 26 4696
|
||
4627 0 0000 30 4583
|
||
4628 0 0001 26 4695
|
||
4629 0 0000 30 4583
|
||
4630 0 0000 41 4689
|
||
4631 0 2211 40 4694
|
||
4632 0 0001 26 4698
|
||
4633 0 0000 30 4583
|
||
4634 0 0000 49 0001
|
||
4635 0 0000 41 4700
|
||
4636 0 0001 49 0001
|
||
4637 0 0001 40 4700
|
||
4638 0 2201 27 4694
|
||
4639 0 0001 26 4697
|
||
4640 0 0000 30 4583
|
||
4641 0 0000 41 4696
|
||
4642 0 0000 37 4668
|
||
4643 0 0000 42 4695
|
||
4644 0 0000 11 4700
|
||
4645 0 0000 49 0008
|
||
4646 0 9999 20 4648
|
||
4647 0 0000 43 0000
|
||
4648 0 0001 13 4702
|
||
4649 0 1210 27 4702
|
||
4650 0 0000 32 4655
|
||
4651 0 0000 43 0000
|
||
4652 0 0000 30 4655
|
||
4653 0 0000 42 4699
|
||
4654 0 0001 21 4663
|
||
4655 0 0000 44 0090
|
||
4656 0 0000 30 0091
|
||
4657 0 0009 33 4665
|
||
4658 0 0412 40 4656
|
||
4659 0 0000 42 4673
|
||
4660 0 0000 46 4694
|
||
4661 0 0070 29 4694
|
||
4662 0 0000 30 0097
|
||
4663 0 0000 40 4702
|
||
4664 0 0000 30 4659
|
||
4665 0 0000 42 4565
|
||
4666 0 0000 46 4693
|
||
4667 1 0000 30 0000
|
||
4668 0 0000 10 4700
|
||
4669 0 0000 42 4698
|
||
4670 0 9999 20 4673
|
||
4671 0 0000 12 4694
|
||
4672 0 0000 22 4694
|
||
4673 0 0000 42 4695
|
||
4674 0 9999 20 4653
|
||
4675 0 0001 43 0000
|
||
4676 0 0000 30 4653
|
||
4677 0 0810 18 4690
|
||
4678 0 0001 35 4583
|
||
4679 0 0000 41 0102
|
||
4680 0 0011 18 4691
|
||
4681 0 0001 35 4583
|
||
4682 0 0000 42 4693
|
||
4683 1 0001 26 0000
|
||
4684 0 0000 30 4665
|
||
4685 0 0000 42 4697
|
||
4686 0 9999 20 4601
|
||
4687 0 0502 27 4703
|
||
4688 0 0000 30 4613
|
||
4689 0 5800 00 0000
|
||
4690 2 8562 45 5563
|
||
4691 2 4955 45 5300
|
||
4692 0 5110 00 0000
|
||
4693 0 0000 00 0000
|
||
4694 0 0000 00 0000
|
||
4695 0 0000 00 0000
|
||
4696 0 0000 00 0000
|
||
4697 0 0000 00 0000
|
||
4698 0 0000 00 0000
|
||
4699 0 0000 00 0000
|
||
4700 0 0000 00 0000
|
||
4701 0 0000 00 0000
|
||
4702 0 0000 00 0000
|
||
4703 0 0000 00 0000
|
||
4256 0 0000 44 4289
|
||
4257 0 0000 30 4285
|
||
4258 2 5945 62 6453
|
||
4259 2 6300 56 6463
|
||
4260 2 0056 46 0059
|
||
4261 2 4155 47 4500
|
||
4262 2 4955 00 0000
|
||
4263 0 0000 44 4289
|
||
4264 0 0000 30 4285
|
||
4265 2 5945 62 6453
|
||
4266 2 6300 64 5544
|
||
4267 2 4546 49 5545
|
||
4268 2 4400 46 5659
|
||
4269 2 0000 00 0000
|
||
4270 0 0000 44 4289
|
||
4271 0 0000 30 4285
|
||
4272 2 5945 62 6453
|
||
4273 2 6300 49 5353
|
||
4274 2 2044 45 4649
|
||
4275 2 5545 44 0046
|
||
4276 2 5659 00 0000
|
||
4277 0 0000 46 4298
|
||
4278 0 0000 44 4289
|
||
4279 0 0000 30 4287
|
||
4280 0 0000 00 0000
|
||
4281 2 0041 59 4963
|
||
4282 2 4854 45 6349
|
||
4283 2 4300 56 6545
|
||
4284 2 5946 53 5666
|
||
4285 0 0001 40 4298
|
||
4286 0 0003 45 0000
|
||
4287 0 0412 40 4292
|
||
4288 0 0000 42 4309
|
||
4289 0 0050 29 0033
|
||
4290 0 1106 44 4310 REED
|
||
4291 0 4293 30 4310 REED
|
||
4292 0 0000 30 0036
|
||
4293 0 0000 00 0000
|
||
4294 0 0000 00 0000
|
||
4295 0 0000 00 0000
|
||
4296 0 0000 00 0000
|
||
4297 0 0000 00 0000
|
||
4298 0 0000 00 0000
|
||
4299 0 0000 00 0000
|
||
4300 0 0000 00 0000
|
||
4301 0 0000 00 0000
|
||
4302 0 0000 00 0000
|
||
4303 0 0000 00 0000
|
||
4304 0 5000 00 0000
|
||
4305 0 0000 00 0000
|
||
4306 0 0000 00 0000
|
||
4307 0 5110 00 0000
|
||
4308 0 0000 00 0000
|
||
4309 0 0000 00 4293
|
||
9947 0 9950 01 0000
|
||
9948 0 4400 28 9947
|
||
9949 0 0000 31 4277 ERROR
|
||
9950 0 0001 40 4305 ERROR
|
||
9951 0 0000 40 4308 ERROR
|
||
9952 0 0000 42 9947
|
||
9953 0 0000 36 9976
|
||
9954 0 0000 10 4307 ERROR
|
||
9955 1 0000 37 0000
|
||
9956 0 0000 40 4306 ERROR
|
||
9957 0 0000 10 4304 ERROR
|
||
9958 0 0000 14 4305 ERROR
|
||
9959 0 0000 40 4305 ERROR
|
||
9960 0 0000 37 9964
|
||
9961 0 0000 10 4306 ERROR
|
||
9962 0 0000 24 4308 ERROR
|
||
9963 0 0000 40 4306 ERROR
|
||
9964 0 0000 41 4305 ERROR
|
||
9965 0 0000 37 9970
|
||
9966 0 0000 10 4308 ERROR
|
||
9967 0 0000 24 4308 ERROR
|
||
9968 0 0000 40 4308 ERROR
|
||
9969 0 0000 30 9957
|
||
9970 0 0000 31 9980
|
||
9971 0 1210 37 9973
|
||
9972 1 0000 30 0000
|
||
9973 0 0000 10 4307 ERROR
|
||
9974 0 0000 25 4306 ERROR
|
||
9975 1 0000 30 0000
|
||
9976 0 0001 18 4280 ERROR
|
||
9977 1 0000 34 0000
|
||
9978 0 0000 41 9984
|
||
9979 0 0000 30 4263 ERROR
|
||
9980 0 0001 45 0000
|
||
9981 1 1210 37 0000
|
||
9982 0 0000 41 9984
|
||
9983 0 0000 30 4256 ERROR
|
||
9984 2 4653 46 6700
|
||
4253 0 0000 01 0000
|
||
4254 0 1000 60 0117
|
||
4255 0 0000 30 4576 READ
|
||
4310 0 0000 01 0000
|
||
4311 0 0000 42 4310
|
||
4312 1 0000 10 9999
|
||
4313 0 0000 48 0006
|
||
4314 0 0410 40 4320
|
||
4315 1 0000 10 9998
|
||
4316 0 4210 40 4320
|
||
4317 0 0000 48 0008
|
||
4318 0 0000 13 4322
|
||
4319 0 0000 19 4320
|
||
4320 0 2000 61 0000
|
||
4321 1 0000 30 0000
|
||
4322 0 0000 00 0001
|
||
COMPILED PROGRAM ENDS AT 1512
|
||
PROGRAM VARIABLES BEGIN AT 4253
|
||
|
||
|
||
|
||
|
||
|
||
4
|
||
|
||
.90000000, 01 .30000000, 01 .30000000, 01 .30000000, 01 .24000000, 02 .30000000, 01
|
||
.10000000, 02 -.20000000, 01 -.20000000, 01 .17000000, 02 .30000000, 01 -.20000000, 01
|
||
.18000000, 02 .10000000, 02 .45000000, 02 .30000000, 01 -.20000000, 01 .10000000, 02
|
||
.10000000, 02 .29000000, 02
|
||
|
||
4
|
||
|
||
.24000000, 02 .17000000, 02 .45000000, 02 .29000000, 02
|
||
|
||
.90000000, 01 .30000000, 01 .30000000, 01 .30000000, 01 .24000000, 02 .33333333, 00
|
||
.90000010, 01 -.29999999, 01 -.29999999, 01 .90000010, 01 .33333333, 00 -.33333328, 00
|
||
.16000001, 02 .80000010, 01 .40000000, 02 .33333333, 00 -.33333328, 00 .50000003, 00
|
||
.40000000, 01 .39999990, 01
|
||
|
||
4 1 2 3 4
|
||
|
||
.10000003, 01 .19999996, 01 .20000000, 01 .99999975, 00
|
||
|
||
.51840008, 04
|
||
|
||
0
|
||
|
||
.24000000, 02 .90000010, 01 .40000000, 02 .39999990, 01
|
||
|
||
.10000003, 01 .19999996, 01 .20000000, 01 .99999975, 00
|
||
|