github.com/pkimpel.retro-220
1
0
Fork 0
mirror of https://github.com/pkimpel/retro-220.git synced 2026-03-02 01:40:31 +00:00
Code Issues Releases Wiki Activity
Files
master
pkimpel.retro-220/software/BALGOL/BALGOL-Examples/Reference-Manual/Example-1/Example-1-Code-List.lst
Paul Kimpel 83a0b62faf Update BALGOL Library and Examples 
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.
2018-08-06 17:59:30 -07:00

2905 lines
186 KiB
Plaintext
Raw Permalink Blame History

This file contains invisible Unicode characters
This file contains invisible Unicode characters that are indistinguishable to humans but may be processed differently by a computer. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.
0200 BAC-220 STANDARD VERSION 2/1/62
0200 COMMENT FIRST EXAMPLE PROGRAM FROM BALGOL MANUAL, MARCH 1963.
0200
0200 J.G. HERRIOT, OF STANFORD UNIVERSITY, HAS WRITTEN THE FOLLOWING
0200 PROGRAM TO DETERMINE AN APROXIMATION OF HARMONIC-BOUNDARY VALUES,
0200 USING ORTHONORMAL FUNCTIONS$
0200
0200 COMMENT THIS PROGRAM FIRST CONSTRUCTS A SET OF ORTHONORMAL FUNCTIONS
0200 AND THEN USES THEM TO FIND AN APPROXIMATION TO THE SOLUTION OF A
0200 HARMONIC BOUNDARY-VALUE PROBLEM$
0200
0200 COMMENT WE FIRST CONSTRUCT THE ORTHONORMAL FUNCTIONS$
0200
0200 INTEGER I, J, K, L, M, N, NU, TH$
0200 ARRAY R(29), HFN(29), DSUM(24), HFCN(5), HFCEN(6),
0200 FA(25,25), A(25,25), B(25,25), HA(47), HAA(24)$
0200 0 0000 30
0201 0 0000 30
0202 0 0000 10 2960 CONST
0203 0 0000 40 2959 I
0204 0 0000 44
0205 0 0000 30
0206 0 0000 10 2960 CONST
0207 0 0000 12 2959 I
0208 0 0000 40 2959 I
0205 0209
0209 0 0000 13 2958 CONST
0210 0 1111 31 0212
0211 0 0000 36
0212 0 0001 33
0213 0 0000 30
0212 0214
0211 0214
0214 0 0000 42 0214
0215 0 0002 20 0201
0216 0 0000 42 2959 I
0217 1 0000 40 4970 R
0213 0219
0204 0218
0218 0 0000 30
0219 0 0009 43 0000
0220 0 0000 30 0201
0200 0221
0221 0 0000 30
0222 0 0000 30
0200 INPUT DATA (FOR I=(1,1,29)$ R(I)), DIMEN(N)$
0223 0 0100 42 0223
0224 0 0002 20 0222
0225 0 0000 40 2957 N
0226 0 0009 43 0000
0227 0 0000 30 0222
0221 0228
0228 0 0000 30
0229 0 0000 30
0230 0 0000 10 2960 CONST
0231 0 0000 40 2959 I
0232 0 0000 44
0233 0 0000 30
0234 0 0000 10 2960 CONST
0235 0 0000 12 2959 I
0236 0 0000 40 2959 I
0233 0237
0237 0 0000 13 2957 N
0238 0 1111 31 0240
0239 0 0000 36
0240 0 0001 33
0241 0 0000 30
0240 0242
0239 0242
0242 0 0000 10 2960 CONST
0243 0 0000 40 2956 J
0244 0 0000 44
0245 0 0000 30
0246 0 0000 10 2960 CONST
0247 0 0000 12 2956 J
0248 0 0000 40 2956 J
0245 0249
0249 0 0000 13 2957 N
0250 0 1111 31 0252
0251 0 0000 36
0252 0 0001 33
0253 0 0000 30
0252 0254
0251 0254
0254 0 0000 10 2959 I
0255 0 0000 14 2955 CONST
0256 0 0001 49 0010
0257 0 0000 12 2956 J
0258 0 0000 40 2954 TEMP
0259 0 0000 42 2954 TEMP
0260 1 0000 10 4256 FA
0261 0 0000 42 0261
0262 0 0002 20 0229
0253 0264
0244 0263
0263 0 0000 30
0241 0265
0232 0264
0264 0 0000 30
0228 OUTPUT FRESULTS (FOR I=(1,1,N)$ FOR J=(1,1,N)$ FA(I,J)),
0265 0 0009 43 0000
0266 0 0000 30 0229
0228 0267
0267 0 0000 30
0268 0 0000 30
0269 0 0000 10 2960 CONST
0270 0 0000 40 2959 I
0271 0 0000 44
0272 0 0000 30
0273 0 0000 10 2960 CONST
0274 0 0000 12 2959 I
0275 0 0000 40 2959 I
0272 0276
0276 0 0000 13 2957 N
0277 0 1111 31 0279
0278 0 0000 36
0279 0 0001 33
0280 0 0000 30
0279 0281
0278 0281
0281 0 0000 10 2960 CONST
0282 0 0000 40 2956 J
0283 0 0000 44
0284 0 0000 30
0285 0 0000 10 2960 CONST
0286 0 0000 12 2956 J
0287 0 0000 40 2956 J
0284 0288
0288 0 0000 13 2957 N
0289 0 1111 31 0291
0290 0 0000 36
0291 0 0001 33
0292 0 0000 30
0291 0293
0290 0293
0293 0 0000 10 2959 I
0294 0 0000 14 2955 CONST
0295 0 0001 49 0010
0296 0 0000 12 2956 J
0297 0 0000 40 2954 TEMP
0298 0 0000 42 2954 TEMP
0299 1 0000 10 3631 A
0300 0 0000 42 0300
0301 0 0002 20 0268
0292 0303
0283 0302
0302 0 0000 30
0280 0304
0271 0303
0303 0 0000 30
0267 ARESULTS (FOR I=(1,1,N)$ FOR J=(1,1,N)$ A(I,J)),
0304 0 0009 43 0000
0305 0 0000 30 0268
0267 0306
0306 0 0000 30
0307 0 0000 30
0308 0 0000 10 2960 CONST
0309 0 0000 40 2959 I
0310 0 0000 44
0311 0 0000 30
0312 0 0000 10 2960 CONST
0313 0 0000 12 2959 I
0314 0 0000 40 2959 I
0311 0315
0315 0 0000 13 2957 N
0316 0 1111 31 0318
0317 0 0000 36
0318 0 0001 33
0319 0 0000 30
0318 0320
0317 0320
0320 0 0000 10 2960 CONST
0321 0 0000 40 2956 J
0322 0 0000 44
0323 0 0000 30
0324 0 0000 10 2960 CONST
0325 0 0000 12 2956 J
0326 0 0000 40 2956 J
0323 0327
0327 0 0000 13 2957 N
0328 0 1111 31 0330
0329 0 0000 36
0330 0 0001 33
0331 0 0000 30
0330 0332
0329 0332
0332 0 0000 10 2959 I
0333 0 0000 14 2955 CONST
0334 0 0001 49 0010
0335 0 0000 12 2956 J
0336 0 0000 40 2954 TEMP
0337 0 0000 42 2954 TEMP
0338 1 0000 10 3006 B
0339 0 0000 42 0339
0340 0 0002 20 0307
0331 0342
0322 0341
0341 0 0000 30
0319 0343
0310 0342
0342 0 0000 30
0306 BRESULTS (FOR I=(1,1,N)$ FOR J=(1,1,N)$ B(I,J)),
0343 0 0009 43 0000
0344 0 0000 30 0307
0306 0345
0345 0 0000 30
0346 0 0000 30
0347 0 0000 10 2953 CONST
0348 0 0000 40 2952 NU
0349 0 0000 44
0350 0 0000 30
0351 0 0000 10 2953 CONST
0352 0 0000 12 2952 NU
0353 0 0000 40 2952 NU
0350 0354
0354 0 0000 40 2954 TEMP
0355 0 0000 10 2957 N
0356 0 0000 13 2960 CONST
0357 0 0000 13 2954 TEMP
0358 0 1111 31 0360
0359 0 0000 36
0360 0 0000 33
0361 0 0000 30
0360 0362
0359 0362
0362 0 0000 10 2951 CONST
0363 0 0000 14 2952 NU
0364 0 0001 49 0010
0364 0 0001 40 2954 TEMP
0365 0 0000 42 2954 TEMP
0366 1 0000 10 2983 HA
0367 0 0000 42 0367
0368 0 0002 20 0346
0361 0370
0349 0369
0369 0 0000 30
0345 COEFFS (FOR NU=(4,4,N-1)$ HA(2NU-1)),
0370 0 0009 43 0000
0371 0 0000 30 0346
0345 0372
0372 0 0000 30
0373 0 0000 30
0374 0 0000 10 2960 CONST
0375 0 0000 40 2950 K
0376 0 0000 44
0377 0 0000 30
0378 0 0000 10 2960 CONST
0379 0 0000 12 2950 K
0380 0 0000 40 2950 K
0377 0381
0381 0 0000 13 2958 CONST
0382 0 1111 31 0384
0383 0 0000 36
0384 0 0001 33
0385 0 0000 30
0384 0386
0383 0386
0386 0 0000 42 2950 K
0387 1 0000 10 4941 HFN
0388 0 0000 42 0388
0389 0 0002 20 0373
0385 0391
0376 0390
0390 0 0000 30
0372 HFNRES (FOR K=(1,1,29)$ HFN(K)),
0391 0 0009 43 0000
0392 0 0000 30 0373
0372 0393
0393 0 0000 30
0394 0 0000 30
0395 0 0000 10 2949 CONST
0396 0 0000 42 0396
0397 0 0002 20 0394
0398 0 0009 43 0000
0399 0 0000 30 0394
0393 0400
0400 0 0000 30
0401 0 0000 30
0402 0 0000 10 2948 TH
0403 0 0100 42 0403
0404 0 0002 20 0401
0405 0 0000 10 2960 CONST
0406 0 0000 40 2950 K
0407 0 0000 44
0408 0 0000 30
0409 0 0000 10 2960 CONST
0410 0 0000 12 2950 K
0411 0 0000 40 2950 K
0408 0412
0412 0 0000 13 2947 CONST
0413 0 1111 31 0415
0414 0 0000 36
0415 0 0001 33
0416 0 0000 30
0415 0417
0414 0417
0417 0 0000 42 2950 K
0418 1 0000 10 4912 HFCN
0419 0 0000 42 0419
0420 0 0002 20 0401
0416 0422
0407 0421
0421 0 0000 30
0393 CRES(CONST), HFCNRES (TH, FOR K=(1,1,5)$ HFCN(K)),
0422 0 0009 43 0000
0423 0 0000 30 0401
0400 0424
0424 0 0000 30
0425 0 0000 30
0426 0 0000 10 2948 TH
0427 0 0100 42 0427
0428 0 0002 20 0425
0429 0 0000 10 2960 CONST
0430 0 0000 40 2950 K
0431 0 0000 44
0432 0 0000 30
0433 0 0000 10 2960 CONST
0434 0 0000 12 2950 K
0435 0 0000 40 2950 K
0432 0436
0436 0 0000 13 2946 CONST
0437 0 1111 31 0439
0438 0 0000 36
0439 0 0001 33
0440 0 0000 30
0439 0441
0438 0441
0441 0 0000 42 2950 K
0442 1 0000 10 4906 HFCEN
0443 0 0000 42 0443
0444 0 0002 20 0425
0440 0446
0431 0445
0445 0 0000 30
0424 HFCENRES(TH, FOR K=(1,1,6)$ HFCEN(K))$
0446 0 0009 43 0000
0447 0 0000 30 0425
0424 0448
0448 0 0000 30
0449 0 0004 20 0800
0450 0 0064 60 1608
0451 0 0006 60 0000
0452 1 0000 00 0449
0448 FORMAT VECTOR (B8,6F16.8,W0),
0453 0 0004 20 4800
0454 2 4659 45 6264
0455 2 5363 62 2346
0456 2 4124 49 2351
0457 3 0414 00 0000
0458 0 0006 60 0300
0459 0 0006 60 0200
0460 1 0000 00 0453
0448 FTITLE (B48,*FRESULTS,FA(I,J)*,W3,W2),
0461 0 0004 20 4800
0462 2 4259 45 6264
0463 2 5363 62 2342
0464 2 2449 23 5104
0465 3 1400 00 0000
0466 0 0006 60 0300
0467 0 0006 60 0200
0468 1 0000 00 0461
0448 BTITLE (B48,*BRESULTS,B(I,J)*,W3,W2),
0469 0 0004 20 4800
0470 2 4159 45 6264
0471 2 5363 62 2341
0472 2 2449 23 5104
0473 3 1400 00 0000
0474 0 0006 60 0300
0475 0 0006 60 0200
0476 1 0000 00 0469
0448 ATITLE (B48,*ARESULTS,A(I,J)*,W3,W2),
0477 0 0004 20 3000
0478 2 4841 24 8855
0479 3 6420 81 0414
0480 0 0006 60 0200
0481 1 0000 00 0477
0448 COEFTITLE (B30,*HA(8NU-1)*,W2),
0482 0 0004 20 4200
0483 2 5759 45 5349
0484 2 5449 55 4159
0485 2 6800 42 5664
0486 2 5544 41 5968
0487 2 0065 41 5364
0488 3 4562 14 0000
0489 0 0006 60 0300
0490 0 0006 60 0200
0491 1 0000 00 0482
0448 BDYVALUES (B42,*PRELIMINARY BOUNDARY VALUES*,W3,W2),
0492 0 0004 20 4300
0493 2 4356 59 5945
0494 2 4363 45 4400
0495 2 4256 64 5544
0496 2 4159 68 0065
0497 2 4153 64 4562
0498 3 1400 00 0000
0499 0 0006 60 0200
0500 1 0000 00 0492
0448 CBDYVALUES (B43,*CORRECTED BOUNDARY VALUES*,W2),
0501 0 0004 20 5000
0502 2 4356 55 6263
0503 3 4155 63 1400
0504 0 0006 60 0200
0505 1 0000 00 0501
0448 CONTITLE (B50,*CONSTANT*,W2),
0506 0 0004 20 0800
0507 0 0004 90 0200
0508 0 0004 20 0600
0509 0 0064 60 1608
0510 0 0006 60 0000
0511 1 0000 00 0506
0448 TABLE (B8,I2,B6,6F16.8,W0),
0512 0 0004 20 4000
0513 2 6348 45 0065
0514 2 4153 64 4562
0515 2 0056 46 0048
0516 2 2459 48 5623
0517 2 6348 04 0049
0518 3 5500 42 1400
0519 0 0006 60 0300
0520 0 0006 60 0200
0521 1 0000 00 0512
0448 TABLEHEAD (B40, *THE VALUES OF H(RHO,TH) IN B*, W3,W2),
0522 0 0004 20 1300
0523 3 5948 56 1400
0524 0 0004 20 0600
0525 3 8003 85 1400
0526 0 0004 20 1300
0527 3 8103 80 1400
0528 0 0004 20 1300
0529 3 8103 85 1400
0530 0 0004 20 1300
0448 TABLELINE (B13,*RHO*,B6,*0.5*,B13,*1.0*,B13,*1.5*,B13,
0531 3 8203 80 1400
0532 0 0004 20 1300
0533 3 8203 85 1400
0534 0 0004 20 1300
0535 3 8303 80 1400
0536 0 0006 60 0000
0537 1 0000 00 0522
0448 *2.0*,B13,*2.5*,B13,*3.0*,W0),
0538 0 0004 20 0800
0539 3 6348 14 0000
0540 0 0006 60 0000
0541 1 0000 00 0538
0448 TABLETH (B8,*TH*,W0)$
0448 0542
0542 START..
0542 READ ($$ DATA)$
0542 0 0000 10 0543
0543 0 0000 01 0201 DATA
0544 0 0000 44 2807 READ
0545 0 0000 30 2807
0546 RDIM..
0546 READ ($$ DIMEN)$
0546 0 0000 10 0547
0547 0 0000 01 0222 DIMEN
0548 0 0000 44 2807 READ
0549 0 0000 30 2807
0550 0 0000 10 2960 CONST
0551 0 0000 40 2959 I
0552 0 0000 44
0553 0 0000 30
0554 0 0000 10 2960 CONST
0555 0 0000 12 2959 I
0556 0 0000 40 2959 I
0553 0557
0557 0 0000 13 2957 N
0558 0 1111 31 0560
0559 0 0000 36
0560 0 0001 33
0550 FOR I = (1,1,N)$
0561 0 0000 30
0560 0562
0559 0562
0562 0 0000 10 2959 I
0563 0 0000 40 2956 J
0564 0 0000 44
0565 0 0000 30
0566 0 0000 10 2953 CONST
0567 0 0000 12 2956 J
0568 0 0000 40 2956 J
0565 0569
0569 0 0000 13 2957 N
0570 0 1111 31 0572
0571 0 0000 36
0572 0 0001 33
0561 FOR J = (I,4,N)$
0573 0 0000 30
0572 0574
0571 0574
0573 BEGIN
0574 0 0000 10 2959 I
0575 0 0000 13 2956 J
0576 0 0000 40 2806 L
0577 L = I-J$ K = I+J$
0577 0 0000 10 2956 J
0578 0 0000 12 2959 I
0579 0 0000 40 2950 K
0580 0 0000 41 2950 K
0581 0 0000 10 4971 R
0582 0 0000 44 2768 FL*FX
0583 0 0000 30 2768
0584 0 0000 41 2950 K
0585 0 0000 40 2954 TEMP
0586 0 0000 10 4988 R
0587 0 0000 44 2768 FL*FX
0588 0 0000 30 2768
0589 0 0000 24 2767 CONST
0580 SUM = R(1)*K + 1.5R(18)*K.COS(0.59341195.L)
0590 0 0000 40 2756 TEMP
0591 0 0000 10 2806 L
0592 0 0000 44 2742 FLOAT
0593 0 0000 30 2742
0594 0 0000 24 2741 CONST
0595 0 0000 44 2757 COS
0596 0 0000 30 2757
0597 0 0000 24 2756 TEMP
0598 0 0000 22 2954 TEMP
0599 0 0000 41 2950 K
0600 0 0000 40 2954 TEMP
0601 0 0000 10 4999 R
0602 0 0000 44 2768 FL*FX
0603 0 0000 30 2768
0604 0 0000 24 2740 CONST
0605 0 0000 40 2756 TEMP
0606 0 0000 10 2806 L
0607 0 0000 44 2742 FLOAT
0608 0 0000 30 2742
0609 0 0000 24 2739 CONST
0595 + 0.5.R(29)*K.COS(0.78539816.L)$
0610 0 0000 44 2757 COS
0611 0 0000 30 2757
0612 0 0000 24 2756 TEMP
0613 0 0000 22 2954 TEMP
0614 0 0000 40 2738 SUM
0615 0 0000 10 2951 CONST
0616 0 0000 40 2737 M
0617 0 0000 44
0618 0 0000 30
0619 0 0000 10 2960 CONST
0620 0 0000 12 2737 M
0621 0 0000 40 2737 M
0618 0622
0622 0 0000 13 2736 CONST
0623 0 1111 31 0625
0624 0 0000 36
0625 0 0001 33
0615 FOR M = (2,1,17)$
0626 0 0000 30
0625 0627
0624 0627
0627 0 0000 41 2950 K
0628 0 0000 42 2737 M
0629 1 0000 10 4970 R
0630 0 0000 44 2768 FL*FX
0631 0 0000 30 2768
0632 0 0000 24 2735 CONST
0633 0 0000 40 2954 TEMP
0634 0 0000 10 2737 M
0635 0 0000 13 2960 CONST
0636 0 0000 44 2742 FLOAT
0637 0 0000 30 2742
0638 0 0000 24 2734 CONST
0639 0 0000 40 2756 TEMP
0640 0 0000 10 2806 L
0641 0 0000 44 2742 FLOAT
0642 0 0000 30 2742
0643 0 0000 24 2756 TEMP
0626 SUM = SUM + 2.0.R(M)*K.COS((M-1).0.034906585.L)$
0644 0 0000 44 2757 COS
0645 0 0000 30 2757
0646 0 0000 24 2954 TEMP
0647 0 0000 22 2738 SUM
0648 0 0000 40 2738 SUM
0626 0650
0617 0649
0649 0 0000 30
0650 0 0000 10 2733 CONST
0651 0 0000 40 2737 M
0652 0 0000 44
0653 0 0000 30
0654 0 0000 10 2960 CONST
0655 0 0000 12 2737 M
0656 0 0000 40 2737 M
0653 0657
0657 0 0000 13 2732 CONST
0658 0 1111 31 0660
0659 0 0000 36
0660 0 0001 33
0650 FOR M = (19,1,28)$
0661 0 0000 30
0660 0662
0659 0662
0662 0 0000 41 2950 K
0663 0 0000 42 2737 M
0664 1 0000 10 4970 R
0665 0 0000 44 2768 FL*FX
0666 0 0000 30 2768
0667 0 0000 40 2954 TEMP
0668 0 0000 10 2737 M
0669 0 0000 13 2731 CONST
0670 0 0000 44 2742 FLOAT
0671 0 0000 30 2742
0672 0 0000 24 2730 CONST
0673 0 0000 22 2741 CONST
0674 0 0000 40 2756 TEMP
0675 0 0000 10 2806 L
0676 0 0000 44 2742 FLOAT
0677 0 0000 30 2742
0678 0 0000 24 2756 TEMP
0661 SUM = SUM + R(M)*K.COS((0.59341195 + (M-18).0.017453293).L)$
0679 0 0000 44 2757 COS
0680 0 0000 30 2757
0681 0 0000 24 2954 TEMP
0682 0 0000 22 2738 SUM
0683 0 0000 40 2738 SUM
0661 0685
0652 0684
0684 0 0000 30
0685 0 0000 10 2959 I
0686 0 0000 14 2955 CONST
0687 0 0001 49 0010
0688 0 0000 12 2956 J
0689 0 0000 40 2954 TEMP
0690 0 0000 10 2950 K
0691 0 0000 44 2742 FLOAT
0692 0 0000 30 2742
0693 0 0000 40 2756 TEMP
0694 0 0000 10 2729 CONST
0695 0 0002 45 0000
0696 0 0000 25 2756 TEMP
0697 0 0000 24 2730 CONST
0685 FA(I,J) = (8.0/K).0.017453293.SUM
0698 0 0000 24 2738 SUM
0699 0 0000 42 2954 TEMP
0700 1 0000 40 4256 FA
0701 END$
0573 0702
0564 0701
0701 0 0000 30
0561 0703
0552 0702
0702 0 0000 30
0703 WRITE ($$ FTITLE)$
0703 0 0000 10 0704
0704 0 0000 01 0453 FTITL
0705 0 0000 44 2489 WRITE
0706 0 0000 30 2489
0707 0 0000 10 0708
0708 0 0000 01 0229 FRESU
0709 0 4400 28 2489
0710 1 0000 40 0000
0707 WRITE ($$ FRESULTS, VECTOR)$
0711 0 0000 10 0712
0712 0 0000 01 0449 VECTO
0713 0 0000 44 2489 WRITE
0714 0 0100 30 2489
0715 0 0000 10 2960 CONST
0716 0 0000 40 2956 J
0717 0 0000 44
0718 0 0000 30
0719 0 0000 10 2960 CONST
0720 0 0000 12 2956 J
0721 0 0000 40 2956 J
0718 0722
0722 0 0000 13 2957 N
0723 0 1111 31 0725
0724 0 0000 36
0725 0 0001 33
0715 FOR J = (1,1,N)$
0726 0 0000 30
0725 0727
0724 0727
0726 B(1,J) = FA(1,J)$
0727 0 0000 42 2956 J
0728 1 0000 10 4281 FA
0729 1 0000 40 3031 B
0726 0731
0717 0730
0730 0 0000 30
0731 0 0000 10 2951 CONST
0732 0 0000 40 2959 I
0733 0 0000 44
0734 0 0000 30
0735 0 0000 10 2960 CONST
0736 0 0000 12 2959 I
0737 0 0000 40 2959 I
0734 0738
0738 0 0000 13 2957 N
0739 0 1111 31 0741
0740 0 0000 36
0741 0 0001 33
0731 FOR I = (2,1,N)$
0742 0 0000 30
0741 0743
0740 0743
0742 BEGIN
0743 0 0000 10 2960 CONST
0744 0 0000 40 2956 J
0745 0 0000 44
0746 0 0000 30
0747 0 0000 10 2960 CONST
0748 0 0000 12 2956 J
0749 0 0000 40 2956 J
0746 0750
0750 0 0000 40 2954 TEMP
0751 0 0000 10 2959 I
0752 0 0000 13 2960 CONST
0753 0 0000 13 2954 TEMP
0754 0 1111 31 0756
0755 0 0000 36
0756 0 0000 33
0743 FOR J = (1,1,I-1)$
0757 0 0000 30
0756 0758
0755 0758
0758 0 0000 10 2959 I
0759 0 0000 14 2955 CONST
0760 0 0001 49 0010
0761 0 0000 12 2956 J
0762 0 0000 40 2954 TEMP
0763 0 0000 10 2956 J
0764 0 0000 14 2955 CONST
0765 0 0001 49 0010
0766 0 0000 12 2959 I
0767 0 0000 40 2756 TEMP
0768 0 0000 10 2956 J
0769 0 0000 14 2955 CONST
0770 0 0001 49 0010
0771 0 0000 12 2956 J
0757 B(I,J) = -B(J,I)/B(J,J)$
0772 0 0000 40 2488 TEMP
0773 0 0000 42 2756 TEMP
0774 1 0000 11 3006 B
0775 0 0002 45 0000
0776 0 0000 42 2488 TEMP
0777 1 0000 25 3006 B
0778 0 0000 42 2954 TEMP
0779 1 0000 40 3006 B
0757 0781
0745 0780
0780 0 0000 30
0781 0 0000 10 2959 I
0782 0 0000 40 2956 J
0783 0 0000 44
0784 0 0000 30
0785 0 0000 10 2960 CONST
0786 0 0000 12 2956 J
0787 0 0000 40 2956 J
0784 0788
0788 0 0000 13 2957 N
0789 0 1111 31 0791
0790 0 0000 36
0791 0 0001 33
0781 FOR J = (I,1,N)$
0792 0 0000 30
0791 0793
0790 0793
0792 BEGIN
0793 0 0000 10 2959 I
0794 0 0000 14 2955 CONST
0795 0 0001 49 0010
0796 0 0000 12 2956 J
0797 0 0000 40 2954 TEMP
0798 0 0000 10 2959 I
0799 0 0000 14 2955 CONST
0800 0 0001 49 0010
0801 0 0000 12 2956 J
0793 B(I,J) = FA(I,J)$
0802 0 0000 40 2488 TEMP
0803 0 0000 42 2488 TEMP
0804 1 0000 10 4256 FA
0805 0 0000 42 2954 TEMP
0806 1 0000 40 3006 B
0807 0 0000 10 2960 CONST
0808 0 0000 40 2950 K
0809 0 0000 44
0810 0 0000 30
0811 0 0000 10 2960 CONST
0812 0 0000 12 2950 K
0813 0 0000 40 2950 K
0810 0814
0814 0 0000 40 2954 TEMP
0815 0 0000 10 2959 I
0816 0 0000 13 2960 CONST
0817 0 0000 13 2954 TEMP
0818 0 1111 31 0820
0819 0 0000 36
0820 0 0000 33
0807 FOR K = (1,1,I-1)$
0821 0 0000 30
0820 0822
0819 0822
0822 0 0000 10 2959 I
0823 0 0000 14 2955 CONST
0824 0 0001 49 0010
0825 0 0000 12 2956 J
0826 0 0000 40 2954 TEMP
0827 0 0000 10 2959 I
0828 0 0000 14 2955 CONST
0829 0 0001 49 0010
0830 0 0000 12 2956 J
0831 0 0000 40 2488 TEMP
0832 0 0000 10 2959 I
0833 0 0000 14 2955 CONST
0834 0 0001 49 0010
0835 0 0000 12 2950 K
0836 0 0000 40 2756 TEMP
0837 0 0000 10 2950 K
0838 0 0000 14 2955 CONST
0839 0 0001 49 0010
0821 B(I,J) = B(I,J) + B(I,K).B(K,J)
0840 0 0000 12 2956 J
0841 0 0000 40 2487 TEMP
0842 0 0000 42 2756 TEMP
0843 1 0000 10 3006 B
0844 0 0000 42 2487 TEMP
0845 1 0000 24 3006 B
0846 0 0000 42 2488 TEMP
0847 1 0000 22 3006 B
0848 0 0000 42 2954 TEMP
0849 1 0000 40 3006 B
0821 0851
0809 0850
0850 0 0000 30
0851 END$
0792 0852
0783 0851
0851 0 0000 30
0852 0 0000 10 2960 CONST
0853 0 0000 40 2956 J
0854 0 0000 44
0855 0 0000 30
0856 0 0000 10 2960 CONST
0857 0 0000 12 2956 J
0858 0 0000 40 2956 J
0855 0859
0859 0 0000 40 2954 TEMP
0860 0 0000 10 2959 I
0861 0 0000 13 2960 CONST
0862 0 0000 13 2954 TEMP
0863 0 1111 31 0865
0864 0 0000 36
0865 0 0000 33
0852 FOR J = (1,1,I-1)$
0866 0 0000 30
0865 0867
0864 0867
0867 0 0000 10 2959 I
0868 0 0000 14 2955 CONST
0869 0 0001 49 0010
0870 0 0000 12 2956 J
0871 0 0000 40 2954 TEMP
0872 0 0000 10 2959 I
0873 0 0000 14 2955 CONST
0874 0 0001 49 0010
0875 0 0000 12 2956 J
0876 0 0000 40 2488 TEMP
0877 0 0000 10 2956 J
0878 0 0000 14 2955 CONST
0879 0 0001 49 0010
0880 0 0000 12 2956 J
0881 0 0000 40 2487 TEMP
0882 0 0000 10 2959 I
0883 0 0000 14 2955 CONST
0884 0 0001 49 0010
0885 0 0000 12 2959 I
0866 B(I,J) = B(I,J).SQRT(B(J,J)/B(I,I))
0886 0 0000 40 2756 TEMP
0887 0 0000 42 2487 TEMP
0888 1 0000 10 3006 B
0889 0 0002 45 0000
0890 0 0000 42 2756 TEMP
0891 1 0000 25 3006 B
0892 0 0000 44 2442 SQRT
0893 0 0000 30 2442
0894 0 0000 42 2488 TEMP
0895 1 0000 24 3006 B
0896 0 0000 42 2954 TEMP
0897 1 0000 40 3006 B
0866 0899
0854 0898
0898 0 0000 30
0892 END$
0742 0900
0733 0899
0899 0 0000 30
0900 0 0000 10 2960 CONST
0901 0 0000 40 2959 I
0902 0 0000 44
0903 0 0000 30
0904 0 0000 10 2960 CONST
0905 0 0000 12 2959 I
0906 0 0000 40 2959 I
0903 0907
0907 0 0000 13 2957 N
0908 0 1111 31 0910
0909 0 0000 36
0910 0 0001 33
0900 FOR I = (1,1,N)$
0911 0 0000 30
0910 0912
0909 0912
0912 0 0000 10 2959 I
0913 0 0000 14 2955 CONST
0914 0 0001 49 0010
0915 0 0000 12 2959 I
0916 0 0000 40 2954 TEMP
0917 0 0000 10 2959 I
0918 0 0000 14 2955 CONST
0919 0 0001 49 0010
0920 0 0000 12 2959 I
0921 0 0000 40 2488 TEMP
0922 0 0000 42 2488 TEMP
0923 1 0000 10 3006 B
0924 0 0000 44 2442 SQRT
0925 0 0000 30 2442
0926 0 0000 40 2488 TEMP
0927 0 0000 10 2959 I
0928 0 0000 44 2742 FLOAT
0929 0 0000 30 2742
0930 0 0000 24 2488 TEMP
0911 B(I,I) = 1.0/(SQRT(B(I,I)).I)$
0931 0 0000 40 2488 TEMP
0932 0 0000 10 2441 CONST
0933 0 0002 45 0000
0934 0 0000 25 2488 TEMP
0935 0 0000 42 2954 TEMP
0936 1 0000 40 3006 B
0911 0938
0902 0937
0937 0 0000 30
0938 WRITE ($$ BTITLE)$
0938 0 0000 10 0939
0939 0 0000 01 0461 BTITL
0940 0 0000 44 2489 WRITE
0941 0 0000 30 2489
0942 0 0000 10 0943
0943 0 0000 01 0307 BRESU
0944 0 4400 28 2489
0945 1 0000 40 0000
0942 WRITE ($$ BRESULTS, VECTOR)$
0946 0 0000 10 0947
0947 0 0000 01 0449 VECTO
0948 0 0000 44 2489 WRITE
0949 0 0100 30 2489
0950 0 0000 10 2960 CONST
0951 0 0000 40 2959 I
0952 0 0000 44
0953 0 0000 30
0954 0 0000 10 2960 CONST
0955 0 0000 12 2959 I
0956 0 0000 40 2959 I
0953 0957
0957 0 0000 13 2957 N
0958 0 1111 31 0960
0959 0 0000 36
0960 0 0001 33
0950 FOR I = (1,1,N)$
0961 0 0000 30
0960 0962
0959 0962
0962 0 0000 10 2960 CONST
0963 0 0000 40 2956 J
0964 0 0000 44
0965 0 0000 30
0966 0 0000 10 2960 CONST
0967 0 0000 12 2956 J
0968 0 0000 40 2956 J
0965 0969
0969 0 0000 13 2957 N
0970 0 1111 31 0972
0971 0 0000 36
0972 0 0001 33
0961 FOR J = (1,1,N)$
0973 0 0000 30
0972 0974
0971 0974
0974 0 0000 10 2959 I
0975 0 0000 14 2955 CONST
0976 0 0001 49 0010
0977 0 0000 12 2956 J
0978 0 0000 40 2954 TEMP
0973 A(I,J) = 0$
0979 0 0000 42 2954 TEMP
0980 1 0000 46 3631 A
0973 0982
0964 0981
0981 0 0000 30
0961 0983
0952 0982
0982 0 0000 30
0983 A(1,1) = B(1,1)$
0983 0 0000 10 3032 B
0984 0 0000 40 3657 A
0985 0 0000 10 2951 CONST
0986 0 0000 40 2959 I
0987 0 0000 44
0988 0 0000 30
0989 0 0000 10 2960 CONST
0990 0 0000 12 2959 I
0991 0 0000 40 2959 I
0988 0992
0992 0 0000 13 2957 N
0993 0 1111 31 0995
0994 0 0000 36
0995 0 0001 33
0985 FOR I = (2,1,N)$
0996 0 0000 30
0995 0997
0994 0997
0996 BEGIN
0997 0 0000 10 2960 CONST
0998 0 0000 40 2956 J
0999 0 0000 44
1000 0 0000 30
1001 0 0000 10 2960 CONST
1002 0 0000 12 2956 J
1003 0 0000 40 2956 J
1000 1004
1004 0 0000 40 2954 TEMP
1005 0 0000 10 2959 I
1006 0 0000 13 2960 CONST
1007 0 0000 13 2954 TEMP
1008 0 1111 31 1010
1009 0 0000 36
1010 0 0000 33
0997 FOR J= (1,1,I-1)$
1011 0 0000 30
1010 1012
1009 1012
1011 BEGIN
1012 0 0000 10 2959 I
1013 0 0000 14 2955 CONST
1014 0 0001 49 0010
1015 0 0000 12 2956 J
1016 0 0000 40 2954 TEMP
1012 A(I,J) = 0$
1017 0 0000 42 2954 TEMP
1018 1 0000 46 3631 A
1019 0 0000 10 2956 J
1020 0 0000 40 2950 K
1021 0 0000 44
1022 0 0000 30
1023 0 0000 10 2960 CONST
1024 0 0000 12 2950 K
1025 0 0000 40 2950 K
1022 1026
1026 0 0000 40 2954 TEMP
1027 0 0000 10 2959 I
1028 0 0000 13 2960 CONST
1029 0 0000 13 2954 TEMP
1030 0 1111 31 1032
1031 0 0000 36
1032 0 0000 33
1019 FOR K = (J,1,I-1)$
1033 0 0000 30
1032 1034
1031 1034
1034 0 0000 10 2959 I
1035 0 0000 14 2955 CONST
1036 0 0001 49 0010
1037 0 0000 12 2956 J
1038 0 0000 40 2954 TEMP
1039 0 0000 10 2959 I
1040 0 0000 14 2955 CONST
1041 0 0001 49 0010
1042 0 0000 12 2956 J
1043 0 0000 40 2488 TEMP
1044 0 0000 10 2959 I
1045 0 0000 14 2955 CONST
1046 0 0001 49 0010
1047 0 0000 12 2950 K
1048 0 0000 40 2756 TEMP
1049 0 0000 10 2950 K
1050 0 0000 14 2955 CONST
1051 0 0001 49 0010
1033 A(I,J) = A(I,J) + B(I,K).A(K,J)
1052 0 0000 12 2956 J
1053 0 0000 40 2487 TEMP
1054 0 0000 42 2756 TEMP
1055 1 0000 10 3006 B
1056 0 0000 42 2487 TEMP
1057 1 0000 24 3631 A
1058 0 0000 42 2488 TEMP
1059 1 0000 22 3631 A
1060 0 0000 42 2954 TEMP
1061 1 0000 40 3631 A
1033 1063
1021 1062
1062 0 0000 30
1063 END$
1011 1064
0999 1063
1063 0 0000 30
1064 0 0000 10 2959 I
1065 0 0000 14 2955 CONST
1066 0 0001 49 0010
1067 0 0000 12 2959 I
1068 0 0000 40 2954 TEMP
1069 0 0000 10 2959 I
1070 0 0000 14 2955 CONST
1071 0 0001 49 0010
1064 A(I,I) = B(I,I)
1072 0 0000 12 2959 I
1073 0 0000 40 2488 TEMP
1074 0 0000 42 2488 TEMP
1075 1 0000 10 3006 B
1076 0 0000 42 2954 TEMP
1077 1 0000 40 3631 A
1078 END$
0996 1079
0987 1078
1078 0 0000 30
1079 WRITE ($$ ATITLE)$
1079 0 0000 10 1080
1080 0 0000 01 0469 ATITL
1081 0 0000 44 2489 WRITE
1082 0 0000 30 2489
1083 0 0000 10 1084
1084 0 0000 01 0268 ARESU
1085 0 4400 28 2489
1086 1 0000 40 0000
1083 WRITE ($$ ARESULTS, VECTOR)$
1083
1087 0 0000 10 1088
1088 0 0000 01 0449 VECTO
1089 0 0000 44 2489 WRITE
1090 0 0100 30 2489
1091 COMMENT NOW CONSTRUCT THE APROXIMATION TO THE SOLUTION$
1091
1091 0 0000 10 2953 CONST
1092 0 0000 40 2956 J
1093 0 0000 44
1094 0 0000 30
1095 0 0000 10 2953 CONST
1096 0 0000 12 2956 J
1097 0 0000 40 2956 J
1094 1098
1098 0 0000 40 2954 TEMP
1099 0 0000 10 2957 N
1100 0 0000 13 2960 CONST
1101 0 0000 13 2954 TEMP
1102 0 1111 31 1104
1103 0 0000 36
1104 0 0000 33
1091 FOR J = (4,4,N-1)$
1105 0 0000 30
1104 1106
1103 1106
1105 BEGIN
1105 DSUM(J) = 0$
1106 0 0000 42 2956 J
1107 1 0000 46 4917 DSUM
1108 0 0000 10 2960 CONST
1109 0 0000 40 2737 M
1110 0 0000 44
1111 0 0000 30
1112 0 0000 10 2960 CONST
1113 0 0000 12 2737 M
1114 0 0000 40 2737 M
1111 1115
1115 0 0000 13 2736 CONST
1116 0 1111 31 1118
1117 0 0000 36
1118 0 0001 33
1108 FOR M = (1,1,17)$
1119 0 0000 30
1118 1120
1117 1120
1120 0 0000 42 2737 M
1121 1 0000 10 4970
1122 1 0000 24 4970 R
1123 0 0000 40 2954 TEMP
1124 1 0000 10 4971
1125 1 0000 24 4971 R
1126 0 0000 22 2954 TEMP
1119 DSUM(J) = DSUM(J) + (R(M)*2 + R(M+1)*2).
1127 0 0000 41 2956 J
1128 0 0000 40 2954 TEMP
1129 1 0000 10 4971 R
1130 0 0000 44 2768 FL*FX
1131 0 0000 30 2768
1119 (R(M+1)*J.SIN(M.0.034906585.J)
1132 0 0000 40 2488 TEMP
1133 0 0000 10 2956 J
1134 0 0000 44 2742 FLOAT
1135 0 0000 30 2742
1136 0 0000 24 2734 CONST
1137 0 0000 40 2487 TEMP
1138 0 0000 10 2737 M
1139 0 0000 44 2742 FLOAT
1140 0 0000 30 2742
1141 0 0000 24 2487 TEMP
1142 0 0000 44 2385 SIN
1143 0 0000 30 2385
1144 0 0000 24 2488 TEMP
1145 0 0000 41 2956 J
1146 0 0000 40 2488 TEMP
1147 0 0000 42 2737 M
1148 1 0000 10 4970 R
1149 0 0000 44 2768 FL*FX
1150 0 0000 30 2768
1151 0 0000 40 2487 TEMP
1152 0 0000 10 2737 M
1153 0 0000 13 2960 CONST
1154 0 0000 44 2742 FLOAT
1155 0 0000 30 2742
1156 0 0000 24 2734 CONST
1157 0 0000 40 2756 TEMP
1158 0 0000 10 2956 J
1159 0 0000 44 2742 FLOAT
1160 0 0000 30 2742
1161 0 0000 24 2756 TEMP
1162 0 0000 44 2385 SIN
1163 0 0000 30 2385
1164 0 0000 24 2487 TEMP
1165 0 0000 23 2488 TEMP
1142 - R(M)*J.SIN((M-1).0.034906585.J))$
1166 0 0000 24 2954 TEMP
1167 0 0000 42 2956 J
1168 1 0000 23 4917 DSUM
1169 0 0000 40 2954 TEMP
1170 0 0000 11 2954 TEMP
1171 1 0000 40 4917 DSUM
1119 1173
1110 1172
1172 0 0000 30
1173 0 0000 10 2731 CONST
1174 0 0000 40 2737 M
1175 0 0000 44
1176 0 0000 30
1177 0 0000 10 2960 CONST
1178 0 0000 12 2737 M
1179 0 0000 40 2737 M
1176 1180
1180 0 0000 13 2732 CONST
1181 0 1111 31 1183
1182 0 0000 36
1183 0 0001 33
1173 FOR M = (18,1,28)$
1184 0 0000 30
1183 1185
1182 1185
1185 0 0000 42 2737 M
1186 1 0000 10 4970
1187 1 0000 24 4970 R
1188 0 0000 40 2954 TEMP
1189 1 0000 10 4971
1190 1 0000 24 4971 R
1184 DSUM(J) = DSUM(J) + (R(M)*2 + R(M+1)*2.(R(M+1)*J.
1191 0 0000 41 2956 J
1192 0 0000 40 2488 TEMP
1193 1 0000 10 4971 R
1194 0 0000 44 2768 FL*FX
1195 0 0000 30 2768
1196 0 0000 40 2487 TEMP
1197 0 0000 10 2737 M
1198 0 0000 13 2736 CONST
1199 0 0000 44 2742 FLOAT
1200 0 0000 30 2742
1201 0 0000 24 2730 CONST
1202 0 0000 22 2741 CONST
1199 SIN((0.59341195 + (M-17).0.017453293).J)
1203 0 0000 40 2756 TEMP
1204 0 0000 10 2956 J
1205 0 0000 44 2742 FLOAT
1206 0 0000 30 2742
1207 0 0000 24 2756 TEMP
1208 0 0000 44 2385 SIN
1209 0 0000 30 2385
1210 0 0000 24 2487 TEMP
1211 0 0000 41 2956 J
1212 0 0000 40 2487 TEMP
1213 0 0000 42 2737 M
1214 1 0000 10 4970 R
1215 0 0000 44 2768 FL*FX
1216 0 0000 30 2768
1208 - R(M)*J.SIN((0.59341195
1217 0 0000 40 2756 TEMP
1218 0 0000 10 2737 M
1219 0 0000 13 2731 CONST
1220 0 0000 44 2742 FLOAT
1221 0 0000 30 2742
1222 0 0000 24 2730 CONST
1223 0 0000 22 2741 CONST
1224 0 0000 40 2384 TEMP
1225 0 0000 10 2956 J
1226 0 0000 44 2742 FLOAT
1227 0 0000 30 2742
1228 0 0000 24 2384 TEMP
1229 0 0000 44 2385 SIN
1230 0 0000 30 2385
1231 0 0000 24 2756 TEMP
1232 0 0000 23 2487 TEMP
1220 + (M-18).0.017453293).J)))
1233 0 0000 24 2488 TEMP
1234 0 0000 23 2954 TEMP
1235 0 0000 42 2956 J
1236 1 0000 23 4917 DSUM
1237 0 0000 40 2954 TEMP
1238 0 0000 11 2954 TEMP
1239 1 0000 40 4917 DSUM
1184 1241
1175 1240
1240 0 0000 30
1235 END$
1105 1242
1093 1241
1241 0 0000 30
1242 0 0000 10 2953 CONST
1243 0 0000 40 2952 NU
1244 0 0000 44
1245 0 0000 30
1246 0 0000 10 2953 CONST
1247 0 0000 12 2952 NU
1248 0 0000 40 2952 NU
1245 1249
1249 0 0000 40 2954 TEMP
1250 0 0000 10 2957 N
1251 0 0000 13 2960 CONST
1252 0 0000 13 2954 TEMP
1253 0 1111 31 1255
1254 0 0000 36
1255 0 0000 33
1242 FOR NU = (4,4,N-1)$
1256 0 0000 30
1255 1257
1254 1257
1256 BEGIN
1257 0 0000 10 2951 CONST
1258 0 0000 14 2952 NU
1259 0 0001 49 0010
1259 0 0001 40 2954 TEMP
1256 HA(2NU-1) = 0$
1260 0 0000 42 2954 TEMP
1261 1 0000 46 2983 HA
1262 0 0000 10 2953 CONST
1263 0 0000 40 2956 J
1264 0 0000 44
1265 0 0000 30
1266 0 0000 10 2953 CONST
1267 0 0000 12 2956 J
1268 0 0000 40 2956 J
1265 1269
1269 0 0000 13 2952 NU
1270 0 1111 31 1272
1271 0 0000 36
1272 0 0001 33
1262 FOR J = (4,4,NU)$
1273 0 0000 30
1272 1274
1271 1274
1274 0 0000 10 2951 CONST
1275 0 0000 14 2952 NU
1276 0 0001 49 0010
1276 0 0001 40 2954 TEMP
1277 0 0000 10 2951 CONST
1278 0 0000 14 2952 NU
1279 0 0001 49 0010
1279 0 0001 40 2488 TEMP
1280 0 0000 10 2952 NU
1281 0 0000 14 2955 CONST
1282 0 0001 49 0010
1283 0 0000 12 2956 J
1284 0 0000 40 2487 TEMP
1273 HA(2NU-1) = HA(2NU-1) + A(NU,J).DSUM(J)$
1285 0 0000 42 2487 TEMP
1286 1 0000 10 3631 A
1287 0 0000 42 2956 J
1288 1 0000 24 4917 DSUM
1289 0 0000 42 2488 TEMP
1290 1 0000 22 2983 HA
1291 0 0000 42 2954 TEMP
1292 1 0000 40 2983 HA
1273 1294
1264 1293
1293 0 0000 30
1294 0 0000 10 2951 CONST
1295 0 0000 14 2952 NU
1296 0 0001 49 0010
1296 0 0001 40 2954 TEMP
1297 0 0000 10 2951 CONST
1298 0 0000 14 2952 NU
1299 0 0001 49 0010
1294 HA(2NU-1) = 4.0.HA(2NU-1)
1299 0 0001 40 2488 TEMP
1300 0 0000 10 2383 CONST
1301 0 0000 42 2488 TEMP
1302 1 0000 24 2983 HA
1303 0 0000 42 2954 TEMP
1304 1 0000 40 2983 HA
1305 END$
1256 1306
1244 1305
1305 0 0000 30
1306 WRITE ($$ COEFTITLE)$
1306 0 0000 10 1307
1307 0 0000 01 0477 COEFT
1308 0 0000 44 2489 WRITE
1309 0 0000 30 2489
1310 0 0000 10 1311
1311 0 0000 01 0346 COEFF
1312 0 4400 28 2489
1313 1 0000 40 0000
1310 WRITE ($$ COEFFS, VECTOR)$
1314 0 0000 10 1315
1315 0 0000 01 0449 VECTO
1316 0 0000 44 2489 WRITE
1317 0 0100 30 2489
1318 0 0000 10 2953 CONST
1319 0 0000 40 2956 J
1320 0 0000 44
1321 0 0000 30
1322 0 0000 10 2953 CONST
1323 0 0000 12 2956 J
1324 0 0000 40 2956 J
1321 1325
1325 0 0000 40 2954 TEMP
1326 0 0000 10 2957 N
1327 0 0000 13 2960 CONST
1328 0 0000 13 2954 TEMP
1329 0 1111 31 1331
1330 0 0000 36
1331 0 0000 33
1318 FOR J = (4,4,N-1)$
1332 0 0000 30
1331 1333
1330 1333
1332 BEGIN
1332 HAA(J) = 0$
1333 0 0000 42 2956 J
1334 1 0000 46 2960 HAA
1335 0 0000 10 2956 J
1336 0 0000 40 2952 NU
1337 0 0000 44
1338 0 0000 30
1339 0 0000 10 2953 CONST
1340 0 0000 12 2952 NU
1341 0 0000 40 2952 NU
1338 1342
1342 0 0000 40 2954 TEMP
1343 0 0000 10 2957 N
1344 0 0000 13 2960 CONST
1345 0 0000 13 2954 TEMP
1346 0 1111 31 1348
1347 0 0000 36
1348 0 0000 33
1335 FOR NU = (J,4,N-1)$
1349 0 0000 30
1348 1350
1347 1350
1350 0 0000 10 2951 CONST
1351 0 0000 14 2952 NU
1352 0 0001 49 0010
1352 0 0001 40 2954 TEMP
1353 0 0000 10 2952 NU
1354 0 0000 14 2955 CONST
1355 0 0001 49 0010
1349 HAA(J) = HAA(J) + HA(2NU-1).A(NU,J)
1356 0 0000 12 2956 J
1357 0 0000 40 2488 TEMP
1358 0 0000 42 2954 TEMP
1359 1 0000 10 2983 HA
1360 0 0000 42 2488 TEMP
1361 1 0000 24 3631 A
1362 0 0000 42 2956 J
1363 1 0000 22 2960 HAA
1364 1 0000 40 2960 HAA
1349 1366
1337 1365
1365 0 0000 30
1366 END$
1332 1367
1320 1366
1366 0 0000 30
1367 0 0000 10 2960 CONST
1368 0 0000 40 2737 M
1369 0 0000 44
1370 0 0000 30
1371 0 0000 10 2960 CONST
1372 0 0000 12 2737 M
1373 0 0000 40 2737 M
1370 1374
1374 0 0000 13 2731 CONST
1375 0 1111 31 1377
1376 0 0000 36
1377 0 0001 33
1367 FOR M = (1,1,18)$
1378 0 0000 30
1377 1379
1376 1379
1378 BEGIN
1378 HFN(M) = 0$
1379 0 0000 42 2737 M
1380 1 0000 46 4941 HFN
1381 0 0000 10 2953 CONST
1382 0 0000 40 2956 J
1383 0 0000 44
1384 0 0000 30
1385 0 0000 10 2953 CONST
1386 0 0000 12 2956 J
1387 0 0000 40 2956 J
1384 1388
1388 0 0000 40 2488 TEMP
1389 0 0000 10 2957 N
1390 0 0000 13 2960 CONST
1391 0 0000 13 2488 TEMP
1392 0 1111 31 1394
1393 0 0000 36
1394 0 0000 33
1381 FOR J = (4,4,N-1)$
1395 0 0000 30
1394 1396
1393 1396
1396 0 0000 41 2956 J
1397 0 0000 42 2737 M
1398 1 0000 10 4970 R
1399 0 0000 44 2768 FL*FX
1400 0 0000 30 2768
1401 0 0000 42 2956 J
1402 1 0000 24 2960 HAA
1403 0 0000 40 2488 TEMP
1404 0 0000 10 2737 M
1405 0 0000 13 2960 CONST
1406 0 0000 44 2742 FLOAT
1407 0 0000 30 2742
1408 0 0000 24 2734 CONST
1395 HFN(M) = HFN(M) + HAA(J).R(M)*J.COS((M-1).0.034906585.J)
1409 0 0000 40 2954 TEMP
1410 0 0000 10 2956 J
1411 0 0000 44 2742 FLOAT
1412 0 0000 30 2742
1413 0 0000 24 2954 TEMP
1414 0 0000 44 2757 COS
1415 0 0000 30 2757
1416 0 0000 24 2488 TEMP
1417 0 0000 42 2737 M
1418 1 0000 22 4941 HFN
1419 1 0000 40 4941 HFN
1395 1421
1383 1420
1420 0 0000 30
1414 END$
1378 1422
1369 1421
1421 0 0000 30
1422 0 0000 10 2733 CONST
1423 0 0000 40 2737 M
1424 0 0000 44
1425 0 0000 30
1426 0 0000 10 2960 CONST
1427 0 0000 12 2737 M
1428 0 0000 40 2737 M
1425 1429
1429 0 0000 13 2958 CONST
1430 0 1111 31 1432
1431 0 0000 36
1432 0 0001 33
1422 FOR M = (19,1,29)$
1433 0 0000 30
1432 1434
1431 1434
1433 BEGIN
1433 HFN(M) = 0$
1434 0 0000 42 2737 M
1435 1 0000 46 4941 HFN
1436 0 0000 10 2953 CONST
1437 0 0000 40 2956 J
1438 0 0000 44
1439 0 0000 30
1440 0 0000 10 2953 CONST
1441 0 0000 12 2956 J
1442 0 0000 40 2956 J
1439 1443
1443 0 0000 40 2488 TEMP
1444 0 0000 10 2957 N
1445 0 0000 13 2960 CONST
1446 0 0000 13 2488 TEMP
1447 0 1111 31 1449
1448 0 0000 36
1449 0 0000 33
1436 FOR J = (4,4,N-1)$
1450 0 0000 30
1449 1451
1448 1451
1450 HFN(M) = HFN(M) + HAA(J).R(M)*J.
1451 0 0000 41 2956 J
1452 0 0000 42 2737 M
1453 1 0000 10 4970 R
1454 0 0000 44 2768 FL*FX
1455 0 0000 30 2768
1456 0 0000 42 2956 J
1457 1 0000 24 2960 HAA
1458 0 0000 40 2488 TEMP
1459 0 0000 10 2737 M
1460 0 0000 13 2731 CONST
1461 0 0000 44 2742 FLOAT
1462 0 0000 30 2742
1463 0 0000 24 2730 CONST
1464 0 0000 22 2741 CONST
1461 COS((0.59341195 + (M-18).0.017453293).J)
1465 0 0000 40 2954 TEMP
1466 0 0000 10 2956 J
1467 0 0000 44 2742 FLOAT
1468 0 0000 30 2742
1469 0 0000 24 2954 TEMP
1470 0 0000 44 2757 COS
1471 0 0000 30 2757
1472 0 0000 24 2488 TEMP
1473 0 0000 42 2737 M
1474 1 0000 22 4941 HFN
1475 1 0000 40 4941 HFN
1450 1477
1438 1476
1476 0 0000 30
1470 END$
1433 1478
1424 1477
1477 0 0000 30
1478 WRITE ($$ BDYVALUES)$
1478 0 0000 10 1479
1479 0 0000 01 0482 BDYVA
1480 0 0000 44 2489 WRITE
1481 0 0000 30 2489
1482 0 0000 10 1483
1483 0 0000 01 0373 HFNRE
1484 0 4400 28 2489
1485 1 0000 40 0000
1482 WRITE ($$ HFNRES, VECTOR)$
1486 0 0000 10 1487
1487 0 0000 01 0449 VECTO
1488 0 0000 44 2489 WRITE
1489 0 0100 30 2489
1490 AVT = 0$
1490 0 0000 46 2382 AVT
1491 0 0000 10 2960 CONST
1492 0 0000 40 2737 M
1493 0 0000 44
1494 0 0000 30
1495 0 0000 10 2960 CONST
1496 0 0000 12 2737 M
1497 0 0000 40 2737 M
1494 1498
1498 0 0000 13 2958 CONST
1499 0 1111 31 1501
1500 0 0000 36
1501 0 0001 33
1491 FOR M = (1,1,29)$
1502 0 0000 30
1501 1503
1500 1503
1503 0 0000 42 2737 M
1504 1 0000 10 4970
1505 1 0000 24 4970 R
1506 0 0000 22 2382 AVT
1502 AVT = AVT + R(M)*2 - HFN(M)$
1507 1 0000 23 4941 HFN
1508 0 0000 40 2382 AVT
1502 1510
1493 1509
1509 0 0000 30
1510 CONST = AVT/29.0$
1510 0 0000 10 2382 AVT
1511 0 0002 45 0000
1512 0 0000 25 2381 CONST
1513 0 0000 40 2949 CONST
1514 WRITE ($$ CONTITLE)$
1514 0 0000 10 1515
1515 0 0000 01 0501 CONTI
1516 0 0000 44 2489 WRITE
1517 0 0000 30 2489
1518 0 0000 10 1519
1519 0 0000 01 0394 CRES
1520 0 4400 28 2489
1521 1 0000 40 0000
1518 WRITE ($$ CRES, VECTOR)$
1522 0 0000 10 1523
1523 0 0000 01 0449 VECTO
1524 0 0000 44 2489 WRITE
1525 0 0100 30 2489
1526 0 0000 10 2960 CONST
1527 0 0000 40 2737 M
1528 0 0000 44
1529 0 0000 30
1530 0 0000 10 2960 CONST
1531 0 0000 12 2737 M
1532 0 0000 40 2737 M
1529 1533
1533 0 0000 13 2958 CONST
1534 0 1111 31 1536
1535 0 0000 36
1536 0 0001 33
1526 FOR M = (1,1,29)$
1537 0 0000 30
1536 1538
1535 1538
1537 HFN(M) = CONST + HFN(M)$
1538 0 0000 42 2737 M
1539 1 0000 10 4941 HFN
1540 0 0000 22 2949 CONST
1541 1 0000 40 4941 HFN
1537 1543
1528 1542
1542 0 0000 30
1543 WRITE ($$ CBDYVALUES)$
1543 0 0000 10 1544
1544 0 0000 01 0492 CBDYV
1545 0 0000 44 2489 WRITE
1546 0 0000 30 2489
1547 0 0000 10 1548
1548 0 0000 01 0373 HFNRE
1549 0 4400 28 2489
1550 1 0000 40 0000
1547 WRITE ($$ HFNRES, VECTOR)$
1551 0 0000 10 1552
1552 0 0000 01 0449 VECTO
1553 0 0000 44 2489 WRITE
1554 0 0100 30 2489
1555 0 0000 10 2960 CONST
1556 0 0000 40 2959 I
1557 0 0000 44
1558 0 0000 30
1559 0 0000 10 2960 CONST
1560 0 0000 12 2959 I
1561 0 0000 40 2959 I
1558 1562
1562 0 0000 13 2947 CONST
1563 0 1111 31 1565
1564 0 0000 36
1565 0 0001 33
1555 FOR I = (1,1,5)$
1566 0 0000 30
1565 1567
1564 1567
1566 BEGIN
1567 0 0000 10 2959 I
1568 0 0000 13 2960 CONST
1569 TH = 5.(I-1)$
1569 0 0000 14 2947 CONST
1570 0 0001 49 0010
1570 0 0001 40 2948 TH
1571 0 0000 10 2960 CONST
1572 0 0000 40 2956 J
1573 0 0000 44
1574 0 0000 30
1575 0 0000 10 2960 CONST
1576 0 0000 12 2956 J
1577 0 0000 40 2956 J
1574 1578
1578 0 0000 13 2947 CONST
1579 0 1111 31 1581
1580 0 0000 36
1581 0 0001 33
1571 FOR J = (1,1,5)$
1582 0 0000 30
1581 1583
1580 1583
1582 BEGIN
1582 HFCN(J) = CONST$
1583 0 0000 10 2949 CONST
1584 0 0000 42 2956 J
1585 1 0000 40 4912 HFCN
1586 0 0000 10 2953 CONST
1587 0 0000 40 2737 M
1588 0 0000 44
1589 0 0000 30
1590 0 0000 10 2953 CONST
1591 0 0000 12 2737 M
1592 0 0000 40 2737 M
1589 1593
1593 0 0000 40 2488 TEMP
1594 0 0000 10 2957 N
1595 0 0000 13 2960 CONST
1596 0 0000 13 2488 TEMP
1597 0 1111 31 1599
1598 0 0000 36
1599 0 0000 33
1586 FOR M = (4,4,N-1)$
1600 0 0000 30
1599 1601
1598 1601
1601 0 0000 10 2956 J
1602 0 0000 44 2742 FLOAT
1603 0 0000 30 2742
1604 0 0000 24 2740 CONST
1605 0 0000 41 2737 M
1606 0 0000 44 2768 FL*FX
1607 0 0000 30 2768
1608 0 0000 42 2737 M
1609 1 0000 24 2960 HAA
1610 0 0000 40 2488 TEMP
1611 0 0000 10 2959 I
1612 0 0000 13 2960 CONST
1613 0 0000 44 2742 FLOAT
1614 0 0000 30 2742
1615 0 0000 24 2380 CONST
1600 HFCN(J) = HFCN(J) + HAA(M).(0.5.J)*M.COS((I-1).0.087266463.M)
1616 0 0000 40 2954 TEMP
1617 0 0000 10 2737 M
1618 0 0000 44 2742 FLOAT
1619 0 0000 30 2742
1620 0 0000 24 2954 TEMP
1621 0 0000 44 2757 COS
1622 0 0000 30 2757
1623 0 0000 24 2488 TEMP
1624 0 0000 42 2956 J
1625 1 0000 22 4912 HFCN
1626 1 0000 40 4912 HFCN
1600 1628
1588 1627
1627 0 0000 30
1621 END$
1582 1629
1573 1628
1628 0 0000 30
1629 WRITE ($$ TABLEHEAD)$
1629 0 0000 10 1630
1630 0 0000 01 0512 TABLE
1631 0 0000 44 2489 WRITE
1632 0 0000 30 2489
1633 WRITE ($$ TABLELINE)$
1633 0 0000 10 1634
1634 0 0000 01 0522 TABLE
1635 0 0000 44 2489 WRITE
1636 0 0000 30 2489
1637 WRITE ($$ TABLETH)$
1637 0 0000 10 1638
1638 0 0000 01 0538 TABLE
1639 0 0000 44 2489 WRITE
1640 0 0000 30 2489
1641 0 0000 10 1642
1642 0 0000 01 0401 HFCNR
1643 0 4400 28 2489
1644 1 0000 40 0000
1641 WRITE ($$ HFCNRES, TABLE)
1645 0 0000 10 1646
1646 0 0000 01 0506 TABLE
1647 0 0000 44 2489 WRITE
1648 0 0100 30 2489
1645 END$
1566 1650
1557 1649
1649 0 0000 30
1650 0 0000 10 2946 CONST
1651 0 0000 40 2959 I
1652 0 0000 44
1653 0 0000 30
1654 0 0000 10 2960 CONST
1655 0 0000 12 2959 I
1656 0 0000 40 2959 I
1653 1657
1657 0 0000 13 2379 CONST
1658 0 1111 31 1660
1659 0 0000 36
1660 0 0001 33
1650 FOR I = (6,1,10)$
1661 0 0000 30
1660 1662
1659 1662
1661 BEGIN
1662 0 0000 10 2959 I
1663 0 0000 13 2960 CONST
1664 TH = 5.(I-1)$
1664 0 0000 14 2947 CONST
1665 0 0001 49 0010
1665 0 0001 40 2948 TH
1666 0 0000 10 2960 CONST
1667 0 0000 40 2956 J
1668 0 0000 44
1669 0 0000 30
1670 0 0000 10 2960 CONST
1671 0 0000 12 2956 J
1672 0 0000 40 2956 J
1669 1673
1673 0 0000 13 2946 CONST
1674 0 1111 31 1676
1675 0 0000 36
1676 0 0001 33
1666 FOR J = (1,1,6)$
1677 0 0000 30
1676 1678
1675 1678
1677 BEGIN
1677 HFCEN(J) = CONST$
1678 0 0000 10 2949 CONST
1679 0 0000 42 2956 J
1680 1 0000 40 4906 HFCEN
1681 0 0000 10 2953 CONST
1682 0 0000 40 2737 M
1683 0 0000 44
1684 0 0000 30
1685 0 0000 10 2953 CONST
1686 0 0000 12 2737 M
1687 0 0000 40 2737 M
1684 1688
1688 0 0000 40 2488 TEMP
1689 0 0000 10 2957 N
1690 0 0000 13 2960 CONST
1691 0 0000 13 2488 TEMP
1692 0 1111 31 1694
1693 0 0000 36
1694 0 0000 33
1681 FOR M = (4,4,N-1)$
1695 0 0000 30
1694 1696
1693 1696
1696 0 0000 10 2956 J
1697 0 0000 44 2742 FLOAT
1698 0 0000 30 2742
1699 0 0000 24 2740 CONST
1695 HFCEN(J) = HFCEN(J) + HAA(M).(0.5.J)*M.
1700 0 0000 41 2737 M
1701 0 0000 44 2768 FL*FX
1702 0 0000 30 2768
1703 0 0000 42 2737 M
1704 1 0000 24 2960 HAA
1705 0 0000 40 2488 TEMP
1706 0 0000 10 2959 I
1707 0 0000 13 2960 CONST
1708 0 0000 44 2742 FLOAT
1709 0 0000 30 2742
1710 0 0000 24 2380 CONST
1708 COS((I-1).0.087266463.M)
1711 0 0000 40 2954 TEMP
1712 0 0000 10 2737 M
1713 0 0000 44 2742 FLOAT
1714 0 0000 30 2742
1715 0 0000 24 2954 TEMP
1716 0 0000 44 2757 COS
1717 0 0000 30 2757
1718 0 0000 24 2488 TEMP
1719 0 0000 42 2956 J
1720 1 0000 22 4906 HFCEN
1721 1 0000 40 4906 HFCEN
1695 1723
1683 1722
1722 0 0000 30
1716 END$
1677 1724
1668 1723
1723 0 0000 30
1724 0 0000 10 1725
1725 0 0000 01 0401 HFCNR
1726 0 4400 28 2489
1727 1 0000 40 0000
1724 WRITE ($$ HFCNRES, TABLE)
1728 0 0000 10 1729
1729 0 0000 01 0506 TABLE
1730 0 0000 44 2489 WRITE
1731 0 0100 30 2489
1728 END$
1661 1733
1652 1732
1732 0 0000 30
1733 STOP 1234$
1733 0 0000 10 2378 CONST
1734 0 0137 00 7310
1735 GO TO RDIM$
1735 0 0000 30 0546 RDIM
1736 FINISH$
1736 0 9669 00 9669
1737 0 1000 60 0000
2960 0 0000 00 0001 POOL
2951 0 0000 00 0002
2958 0 0000 00 0029
2947 0 0000 00 0005
2946 0 0000 00 0006
2955 0 0000 00 0025
2953 0 0000 00 0004
2731 0 0000 00 0018
2736 0 0000 00 0017
2733 0 0000 00 0019
2732 0 0000 00 0028
2379 0 0000 00 0010
2378 0 0000 00 1234
2441 0 5110 00 0000
2735 0 5120 00 0000
2767 0 5115 00 0000
2741 0 5059 34 1195
2740 0 5050 00 0000
2739 0 5078 53 9816
2734 0 4934 90 6585
2730 0 4917 45 3293
2729 0 5180 00 0000
2383 0 5140 00 0000
2381 0 5229 00 0000
2380 0 4987 26 6463
2742 0 0006 45 0000
2743 0 2200 36 2746
2744 0 0000 48 0001
2745 0 0001 20 2743
2746 0 0001 48 0008
2747 1 0000 10 2752
2748 0 0001 49 0008
2749 0 0000 22 2755
2750 0 0000 42 2742
2751 1 0000 30 0000
2752 0 0000 00 0058
2753 0 0000 00 0059
2754 0 0000 00 0060
2755 0 5800 00 0000
2489 0 0100 01 0000
2490 0 0410 40 2525
2491 0 0000 42 2493
2492 0 0010 29 2725
2493 0 0300 29 0101
2494 0 0000 42 2489
2495 0 0412 40 2574
2496 1 0000 41 9999
2497 0 0009 43 7557
2498 0 2200 37 2502
2499 0 0000 42 0100
2500 1 0000 44 0000
2501 1 2201 37 0001
2502 0 0000 40 0132
2503 0 0412 40 2556
2504 0 0000 41 0130
2505 0 0000 37 2523
2506 0 0000 10 0134
2507 0 0000 46 0133
2508 0 5246 37 2685
2509 0 5267 37 2661
2510 0 5249 37 2599
2511 0 5262 37 2712
2512 0 5241 37 2514
2513 0 0000 30 2700
2514 0 0000 41 2532
2515 0 0001 27 0129
2516 0 0000 44 2657
2517 0 0000 30 2648
2518 0 0000 10 0129
2519 0 0000 36 2539
2520 0 2299 37 2553
2521 0 0001 49 0010
2522 0 0000 30 2515
2523 0 0000 42 2490
2524 0 0000 31 2332 ERROR
2525 0 0000 10 9999
2526 0 0401 26 2525
2527 0 0000 33 2541
2528 0 0000 41 2532
2529 0 0001 33 2557
2530 0 2214 36 2525
2531 0 0000 44 2657
2532 0 9900 30 2648
2533 0 2299 37 2525
2534 0 0001 49 0010
2535 0 0000 30 2530
2536 0 0000 10 0129
2537 0 0000 19 0125
2538 0 0000 19 0125
2539 0 0000 10 0130
2540 0 3300 36 2523
2541 0 3300 36 2543
2542 0 0000 13 2726
2543 0 0000 40 0130
2544 0 0210 40 0131
2545 0 0001 48 0002
2546 0 0310 40 0129
2547 0 7242 36 2536
2548 0 7266 36 2577
2549 0 7257 36 2585
2550 0 7263 36 2587
2551 0 7243 36 2576
2552 0 0300 36 2525
2553 0 0000 10 0132
2554 0 0000 40 0134
2555 0 0009 33 2536
2556 0 0000 30 9999
2557 0 0000 42 2525
2558 0 6600 36 2573
2559 0 0001 49 0003
2560 0 0300 37 2567
2561 0 0000 13 2726
2562 0 3300 36 2569
2563 1 6301 27 9999
2564 1 0000 41 9999
2565 0 0411 40 2525
2566 0 0000 30 2525
2567 1 3310 40 9999
2568 0 0000 30 2561
2569 1 0000 10 9999
2570 0 0001 48 0003
2571 1 6610 40 9999
2572 0 0000 30 2525
2573 0 0000 10 0132
2574 0 0009 33 9999
2575 0 0000 30 2564
2576 0 1008 61 0124
2577 0 0001 48 0003
2578 0 3111 40 2579
2579 0 2408 44 2365 REED
2580 0 0101 30 2365 REED
2581 0 0000 42 2583
2582 0 0010 29 2725
2583 0 0240 29 0101
2584 0 0000 30 2539
2585 0 1008 61 0124
2586 0 0000 30 2581
2587 0 0300 36 2592
2588 0 0010 09 2591
2589 0 0001 12 2569
2590 0 0000 30 2587
2591 2 1602 02 0202
2592 0 0000 10 0125
2593 0 0000 36 2539
2594 0 0000 12 2662
2595 0 0001 48 0004
2596 0 3211 40 2597
2597 0 0990 09 0101
2598 0 0000 30 2581
2599 0 0000 42 2609
2600 0 0000 36 2658
2601 0 1100 36 2659
2602 0 0000 40 0134
2603 0 0003 45 0080
2604 0 0002 40 0128
2605 0 0000 46 0126
2606 0 0410 40 0126
2607 0 0001 40 0127
2608 0 0000 10 0134
2609 0 0002 48 0010
2610 0 0001 43 0023
2611 0 0000 12 0129
2612 0 0000 13 0128
2613 0 0001 33 2700
2614 0 0000 19 0125
2615 0 5000 19 0125
2616 0 0000 10 0134
2617 0 0000 41 2653
2618 0 0000 44 2657
2619 0 0001 33 2647
2620 0 0001 27 0128
2621 0 0000 32 2633
2622 0 0000 41 0127
2623 0 0000 37 2539
2624 0 0000 41 2610
2625 0 0000 44 2657
2626 0 0000 30 2647
2627 0 0000 10 0127
2628 0 0000 13 2615
2629 0 0000 42 2545
2630 0 0000 46 0129
2631 0 0003 26 0129
2632 0 0000 30 2602
2633 0 0001 27 0126
2634 0 0000 41 0126
2635 0 0000 37 2641
2636 0 0000 41 0133
2637 0 0000 37 2643
2638 0 0001 27 0133
2639 0 0000 41 2603
2640 0 0000 30 2647
2641 0 0000 41 2570
2642 0 0000 30 2647
2643 0 0000 10 0134
2644 0 0000 41 2662
2645 0 0001 49 0001
2646 0 0000 40 0134
2647 0 0001 49 0018
2648 0 0001 48 0008
2649 0 0001 42 0125
2650 1 0000 49 0008
2651 0 9400 28 0125
2652 0 0024 21 2657
2653 0 0000 43 0020
2654 1 0000 12 0125
2655 1 0010 40 0125
2656 0 0002 26 0125
2657 0 0000 30 9999
2658 0 0009 21 2603
2659 0 0000 49 0051
2660 0 0001 21 2601
2661 0 0002 45 7557
2662 0 0001 48 0008
2663 0 0410 18 2727
2664 0 0001 40 0134
2665 0 0001 34 2675
2666 0 0410 40 2673
2667 0 0449 27 2673
2668 0 0000 42 2673
2669 0 0000 10 0130
2670 0 0001 48 0004
2671 0 4211 40 2674
2672 0 0000 10 2673
2673 0 0002 45 9999
2674 0 0000 20 2604
2675 0 0000 43 7557
2676 0 0000 13 2727
2677 0 0210 40 0133
2678 0 0211 18 0130
2679 0 0000 35 2682
2680 0 0000 10 2727
2681 0 0000 30 2666
2682 0 0000 10 2659
2683 0 1110 40 0134
2684 0 0000 30 2666
2685 0 0304 27 0129
2686 0 0000 32 2689
2687 0 0000 30 2699
2688 0 0000 10 2728
2689 0 0000 36 2688
2690 0 0002 45 7557
2691 0 2210 40 2614
2692 0 0001 48 0008
2693 0 0001 40 0134
2694 0 0002 48 0010
2695 0 0000 12 0130
2696 0 0000 17 2659
2697 0 0310 18 0129
2698 0 0001 34 2708
2699 0 0004 26 0129
2700 0 0000 10 0129
2701 0 0000 19 0125
2702 0 0000 19 0125
2703 0 0002 27 0125
2704 0 0000 41 2715
2705 0 0000 44 2657
2706 0 0000 30 2647
2707 0 0000 30 2539
2708 0 0000 41 2614
2709 0 0000 10 2645
2710 0 0000 42 0131
2711 0 0001 20 2604
2712 0 0002 45 7557
2713 0 0001 48 0008
2714 0 0001 40 0134
2715 0 0000 43 0014
2716 0 0000 13 2727
2717 0 0002 45 7557
2718 0 0001 33 2723
2719 0 0210 18 0130
2720 0 0000 34 2700
2721 0 0000 12 2645
2722 0 0000 30 2710
2723 0 0010 40 0133
2724 0 0000 30 2709
2725 2 0000 00 0000
2726 0 0010 00 0000
2727 0 0000 00 0050
2728 0 5000 00 0000
2807 0 2935 01 0000
2808 0 0410 40 2897
2809 0 0410 40 2898
2810 0 0401 26 2898
2811 0 0000 44 2904
2812 0 0000 30 2897
2813 0 0000 42 2935
2814 1 0000 46 0000
2815 0 0000 46 2943
2816 0 0000 44 2308 RITE
2817 0 0117 30 2308 RITE
2818 0 0000 46 2945
2819 0 0502 26 2945
2820 0 0000 42 2935
2821 0 9999 20 2825
2822 0 6400 28 2943
2823 0 0000 10 0101
2824 0 9999 20 2919
2825 0 9400 28 2945
2826 1 0000 10 0101
2827 0 0000 42 2945
2828 1 0000 49 0000
2829 0 0502 26 2945
2830 0 0000 41 2943
2831 0 0000 37 2845
2832 0 0160 21 2816
2833 0 4400 28 2943
2834 0 2213 36 2860
2835 0 0000 48 0008
2836 1 0000 49 0000
2837 0 0000 19 2942
2838 0 4102 27 2943
2839 0 0000 32 2825
2840 0 0000 10 2942
2841 0 0000 44 2904
2842 0 0000 30 2897
2843 0 5212 27 2943
2844 0 0000 30 2825
2845 0 0162 21 2816
2846 0 1108 36 2876
2847 0 2203 36 2872
2848 0 2220 36 2870
2849 0 2234 36 2870
2850 0 2223 36 2865
2851 0 2213 36 2927
2852 0 2214 36 2858
2853 0 0000 42 2939
2854 0 9999 20 2825
2855 0 0000 44 2904
2856 0 0000 30 2883
2857 0 0000 30 2825
2858 0 9416 26 2945
2859 0 0000 30 2853
2860 0 0000 46 2943
2861 0 0000 10 2942
2862 0 0000 44 2904
2863 0 9992 20 2897
2864 0 0000 30 2825
2865 0 0001 26 2941
2866 0 0000 44 2904
2867 0 0000 30 2883
2868 0 0001 26 2938
2869 0 0000 30 2825
2870 0 0001 26 2937
2871 0 0000 30 2825
2872 0 0000 41 2931
2873 0 2211 40 2936
2874 0 0001 26 2940
2875 0 0000 30 2825
2876 0 0000 49 0001
2877 0 0000 41 2942
2878 0 0001 49 0001
2879 0 0001 40 2942
2880 0 2201 27 2936
2881 0 0001 26 2939
2882 0 0000 30 2825
2883 0 0000 41 2938
2884 0 0000 37 2910
2885 0 0000 42 2937
2886 0 0000 11 2942
2887 0 0000 49 0008
2888 0 9999 20 2890
2889 0 0000 43 0000
2890 0 0001 13 2944
2891 0 1210 27 2944
2892 0 0000 32 2897
2893 0 0000 43 0000
2894 0 0000 30 2897
2895 0 0000 42 2941
2896 0 0001 21 2905
2897 0 0000 44 0090
2898 0 0000 30 0091
2899 0 0009 33 2907
2900 0 0412 40 2898
2901 0 0000 42 2915
2902 0 0000 46 2936
2903 0 0070 29 2936
2904 0 0000 30 0097
2905 0 0000 40 2944
2906 0 0000 30 2901
2907 0 0000 42 2807
2908 0 0000 46 2935
2909 1 0000 30 0000
2910 0 0000 10 2942
2911 0 0000 42 2940
2912 0 9999 20 2915
2913 0 0000 12 2936
2914 0 0000 22 2936
2915 0 0000 42 2937
2916 0 9999 20 2895
2917 0 0001 43 0000
2918 0 0000 30 2895
2919 0 0810 18 2932
2920 0 0001 35 2825
2921 0 0000 41 0102
2922 0 0011 18 2933
2923 0 0001 35 2825
2924 0 0000 42 2935
2925 1 0001 26 0000
2926 0 0000 30 2907
2927 0 0000 42 2939
2928 0 9999 20 2843
2929 0 0502 27 2945
2930 0 0000 30 2855
2931 0 5800 00 0000
2932 2 8562 45 5563
2933 2 4955 45 5300
2934 0 5110 00 0000
2935 0 0000 00 0000
2936 0 0000 00 0000
2937 0 0000 00 0000
2938 0 0000 00 0000
2939 0 0000 00 0000
2940 0 0000 00 0000
2941 0 0000 00 0000
2942 0 0000 00 0000
2943 0 0000 00 0000
2944 0 0000 00 0000
2945 0 0000 00 0000
2311 0 0000 44 2344
2312 0 0000 30 2340
2313 2 5945 62 6453
2314 2 6300 56 6463
2315 2 0056 46 0059
2316 2 4155 47 4500
2317 2 4955 00 0000
2318 0 0000 44 2344
2319 0 0000 30 2340
2320 2 5945 62 6453
2321 2 6300 64 5544
2322 2 4546 49 5545
2323 2 4400 46 5659
2324 2 0000 00 0000
2325 0 0000 44 2344
2326 0 0000 30 2340
2327 2 5945 62 6453
2328 2 6300 49 5353
2329 2 2044 45 4649
2330 2 5545 44 0046
2331 2 5659 00 0000
2332 0 0000 46 2353
2333 0 0000 44 2344
2334 0 0000 30 2342
2335 0 0000 00 0000
2336 2 0041 59 4963
2337 2 4854 45 6349
2338 2 4300 56 6545
2339 2 5946 53 5666
2340 0 0001 40 2353
2341 0 0003 45 0000
2342 0 0412 40 2347
2343 0 0000 42 2364
2344 0 0050 29 0033
2345 0 1106 44 2365 REED
2346 0 2348 30 2365 REED
2347 0 0000 30 0036
2348 0 0000 00 0000
2349 0 0000 00 0000
2350 0 0000 00 0000
2351 0 0000 00 0000
2352 0 0000 00 0000
2353 0 0000 00 0000
2354 0 0000 00 0000
2355 0 0000 00 0000
2356 0 0000 00 0000
2357 0 0000 00 0000
2358 0 0000 00 0000
2359 0 5000 00 0000
2360 0 0000 00 0000
2361 0 0000 00 0000
2362 0 5110 00 0000
2363 0 0000 00 0000
2364 0 0000 00 2348
2442 0 0004 45 0000
2443 0 0000 42 2442
2444 1 0000 36 0000
2445 0 0001 33 2476
2446 0 0810 40 2478
2447 0 2310 40 2479
2448 0 0000 10 2359 ERROR
2449 0 0000 14 2479
2450 0 0000 13 2480
2451 0 2310 40 2479
2452 0 3110 18 2479
2453 0 0000 35 2456
2454 0 0001 33 2456
2455 0 2201 26 2479
2456 0 0000 10 2481
2457 0 0000 24 2478
2458 0 0000 22 2482
2459 0 0000 24 2478
2460 0 0000 22 2483
2461 0 0000 40 2360 ERROR
2462 0 0000 10 2478
2463 0 0000 25 2360 ERROR
2464 0 0000 22 2360 ERROR
2465 0 0000 24 2484
2466 0 0000 40 2360 ERROR
2467 0 0000 10 2478
2468 0 0000 25 2360 ERROR
2469 0 0000 22 2360 ERROR
2470 0 0000 12 2479
2471 0 0000 35 2474
2472 0 0000 24 2485
2473 1 0000 30 0000
2474 0 0000 24 2484
2475 1 0000 30 0000
2476 0 0000 41 2486
2477 0 0000 30 2318 ERROR
2478 0 5110 00 0000
2479 0 0000 00 0000
2480 0 2550 00 0000
2481 1 4916 45 0338
2482 0 5041 11 7101
2483 0 5062 69 7923
2484 0 5050 00 0000
2485 0 5015 81 1388
2486 2 6258 59 6300
2768 0 2771 01 0000
2769 0 4400 28 2768
2770 0 0000 31 2332 ERROR
2771 0 0001 40 2360 ERROR
2772 0 0000 40 2363 ERROR
2773 0 0000 42 2768
2774 0 0000 36 2797
2775 0 0000 10 2362 ERROR
2776 1 0000 37 0000
2777 0 0000 40 2361 ERROR
2778 0 0000 10 2359 ERROR
2779 0 0000 14 2360 ERROR
2780 0 0000 40 2360 ERROR
2781 0 0000 37 2785
2782 0 0000 10 2361 ERROR
2783 0 0000 24 2363 ERROR
2784 0 0000 40 2361 ERROR
2785 0 0000 41 2360 ERROR
2786 0 0000 37 2791
2787 0 0000 10 2363 ERROR
2788 0 0000 24 2363 ERROR
2789 0 0000 40 2363 ERROR
2790 0 0000 30 2778
2791 0 0000 31 2801
2792 0 1210 37 2794
2793 1 0000 30 0000
2794 0 0000 10 2362 ERROR
2795 0 0000 25 2361 ERROR
2796 1 0000 30 0000
2797 0 0001 18 2335 ERROR
2798 1 0000 34 0000
2799 0 0000 41 2805
2800 0 0000 30 2318 ERROR
2801 0 0001 45 0000
2802 1 1210 37 0000
2803 0 0000 41 2805
2804 0 0000 30 2311 ERROR
2805 2 4653 46 6700
2385 0 0006 45 0000
2386 0 0000 41 2431
2387 0 0411 40 2429
2388 0 0000 46 2363 ERROR
2389 0 0000 24 2434
2390 0 0001 49 0002
2391 0 0211 40 2433
2392 0 0000 42 2433
2393 0 9949 20 2399
2394 0 0007 21 2428
2395 1 0001 49 0018
2396 0 0000 40 2363 ERROR
2397 0 0000 10 2398
2398 0 0001 49 0050
2399 0 0001 49 0018
2400 0 0000 40 2361 ERROR
2401 0 0000 22 2361 ERROR
2402 0 0010 18 2362 ERROR
2403 0 0001 34 2407
2404 0 0000 43 0000
2405 0 0000 23 2435
2406 0 0001 49 0000
2407 0 0000 40 2361 ERROR
2408 0 0000 24 2361 ERROR
2409 0 0000 40 2360 ERROR
2410 0 0000 24 2436
2411 0 0000 22 2437
2412 0 0000 24 2360 ERROR
2413 0 0000 22 2438
2414 0 0000 24 2360 ERROR
2415 0 0000 22 2439
2416 0 0000 24 2360 ERROR
2417 0 0000 22 2440
2418 0 0000 24 2361 ERROR
2419 0 0000 22 2361 ERROR
2420 0 0000 40 2361 ERROR
2421 0 0000 10 2359 ERROR
2422 0 0000 14 2363 ERROR
2423 0 0000 10 2361 ERROR
2424 0 0000 42 2385
2425 1 0000 37 0000
2426 0 0000 11 2361 ERROR
2427 1 0000 30 0000
2428 0 0000 42 2385
2429 0 0000 41 2432
2430 0 0000 30 2325 ERROR
2431 0 0000 00 2432
2432 2 6249 55 0000
2433 0 0000 00 0000
2434 0 5031 83 0989
2435 0 5120 00 0000
2436 0 4715 14 8419
2437 1 4846 73 7656
2438 0 4979 68 9679
2439 1 5064 59 6371
2440 0 5057 07 9632
2757 0 0006 45 0000
2758 0 0000 41 2757
2759 0 0411 40 2385 SIN
2760 0 0000 41 2764
2761 0 0000 22 2766
2762 0 0000 22 2766
2763 0 0000 30 2387 SIN
2764 0 0000 00 2765
2765 2 4356 62 0000
2766 0 5078 53 9816
2308 0 0000 01 0000
2309 0 1000 60 0117
2310 0 0000 30 2818 READ
2365 0 0000 01 0000
2366 0 0000 42 2365
2367 1 0000 10 9999
2368 0 0000 48 0006
2369 0 0410 40 2375
2370 1 0000 10 9998
2371 0 4210 40 2375
2372 0 0000 48 0008
2373 0 0000 13 2377
2374 0 0000 19 2375
2375 0 2000 61 0000
2376 1 0000 30 0000
2377 0 0000 00 0001
COMPILED PROGRAM ENDS AT 1737
PROGRAM VARIABLES BEGIN AT 2308
FRESULTS,FA(I,J)
.46363974, 00 .00000000, 00 .00000000, 00 .00000000, 00 .23764104,-01 .00000000, 00
.00000000, 00 .00000000, 00 -.71568347,-02 .00000000, 00 .00000000, 00 .10859770, 00
.00000000, 00 .00000000, 00 .00000000, 00 .11793551,-01 .00000000, 00 .00000000, 00
.00000000, 00 -.46675578,-02 .00000000, 00 .00000000, 00 .45596588,-01 .00000000, 00
.00000000, 00 .00000000, 00 .66280898,-02 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .23645548,-01 .00000000, 00 .00000000, 00
.00000000, 00 .40349226,-02 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .13723587,-01 .00000000, 00 .00000000, 00 .00000000, 00
.25944878,-02 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .85504586,-02 .00000000, 00 .00000000, 00 .00000000, 00 .17352776,-02
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.55925943,-02 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .37885967,-02
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .26349488,-02 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .18701734,-02
BRESULTS,B(I,J)
.14686208, 01 .00000000, 00 .00000000, 00 .00000000, 00 .23764104,-01 .00000000, 00
.00000000, 00 .00000000, 00 -.71568347,-02 .00000000, 00 .00000000, 00 .15172589, 01
.00000000, 00 .00000000, 00 .00000000, 00 .11793551,-01 .00000000, 00 .00000000, 00
.00000000, 00 -.46675578,-02 .00000000, 00 .00000000, 00 .15610347, 01 .00000000, 00
.00000000, 00 .00000000, 00 .66280898,-02 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .16257932, 01 .00000000, 00 .00000000, 00
.00000000, 00 .40349226,-02 .00000000, 00 .00000000, 00 -.31208993, 00 .00000000, 00
.00000000, 00 .00000000, 00 .17884578, 01 .00000000, 00 .00000000, 00 .00000000, 00
.29613152,-02 .00000000, 00 .00000000, 00 -.41973587, 00 .00000000, 00 .00000000, 00
.00000000, 00 .19547473, 01 .00000000, 00 .00000000, 00 .00000000, 00 .22421675,-02
.00000000, 00 .00000000, 00 -.45621904, 00 .00000000, 00 .00000000, 00 .00000000, 00
.20996803, 01 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 -.47127538, 00 .00000000, 00 .00000000, 00 .00000000, 00 .22450411, 01
.00000000, 00 .00000000, 00 .24615538, 00 .00000000, 00 .00000000, 00 .00000000, 00
-.62017180, 00 .00000000, 00 .00000000, 00 .00000000, 00 .26021732, 01 .00000000, 00
.00000000, 00 .45290409, 00 .00000000, 00 .00000000, 00 .00000000, 00 -.84088537, 00
.00000000, 00 .00000000, 00 .00000000, 00 .31976202, 01
ARESULTS,A(I,J)
.14686208, 01 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .15172589, 01
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .15610347, 01 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .16257932, 01 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 -.45834176, 00 .00000000, 00
.00000000, 00 .00000000, 00 .17884578, 01 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 .00000000, 00 -.63684798, 00 .00000000, 00 .00000000, 00
.00000000, 00 .19547473, 01 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 .00000000, 00 -.71217375, 00 .00000000, 00 .00000000, 00 .00000000, 00
.20996803, 01 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00 .00000000, 00
.00000000, 00 -.76619630, 00 .00000000, 00 .00000000, 00 .00000000, 00 .22450411, 01
.00000000, 00 .00000000, 00 .64575954, 00 .00000000, 00 .00000000, 00 .00000000, 00
-.11091510, 01 .00000000, 00 .00000000, 00 .00000000, 00 .26021732, 01 .00000000, 00
.00000000, 00 .12226889, 01 .00000000, 00 .00000000, 00 .00000000, 00 -.16437184, 01
.00000000, 00 .00000000, 00 .00000000, 00 .31976202, 01
HA(8NU-1)
.50545980, 01 .42274664, 01
PRELIMINARY BOUNDARY VALUES
.49796109,-03 .79116912,-02 .39293030,-01 .12059672, 00 .28013604, 00 .52195968, 00
.74267286, 00 .54385490, 00 -.10833408, 01 .15377240,-03 .12638557,-03 .35622571,-04
-.14939181,-03 -.46393693,-03 -.94592558,-03 -.16334454,-02 -.25616621,-02 -.37591919,-02
-.49674820,-02 -.64370066,-02 -.81955319,-02 -.10268529,-01 -.12678090,-01 -.15441817,-01
-.18571706,-01 -.22073060,-01 -.25943490,-01 -.30172023,-01 -.34738354,-01
CONSTANT
.93175934,-01
CORRECTED BOUNDARY VALUES
.93673895,-01 .10108762, 00 .13246896, 00 .21377265, 00 .37331197, 00 .61513561, 00
.83584879, 00 .63703083, 00 -.99016490, 00 .93329706,-01 .93302319,-01 .93211556,-01
.93026543,-01 .92711998,-01 .92230009,-01 .91542489,-01 .90614272,-01 .89416743,-01
.88208452,-01 .86738928,-01 .84980403,-01 .82907405,-01 .80497844,-01 .77734117,-01
.74604228,-01 .71102874,-01 .67232444,-01 .63003911,-01 .58437580,-01
THE VALUES OF H(RHO,TH) IN B
RHO 0.5 1.0 1.5 2.0 2.5 3.0
TH
0 .44141586, 00 .14562673, 02 .26853738, 03 .25094056, 04 .14676438, 05
THE VALUES OF H(RHO,TH) IN B
RHO 0.5 1.0 1.5 2.0 2.5 3.0
TH
5 .41397669, 00 .12041991, 02 .21011010, 03 .19361709, 04 .11276597, 05
THE VALUES OF H(RHO,TH) IN B
RHO 0.5 1.0 1.5 2.0 2.5 3.0
TH
10 .33798099, 00 .55551219, 01 .61639132, 02 .48302088, 03 .26638265, 04
THE VALUES OF H(RHO,TH) IN B
RHO 0.5 1.0 1.5 2.0 2.5 3.0
TH
15 .23022240, 00 -.21629082, 01 -.10892441, 03 -.11749041, 04 -.71435976, 04
THE VALUES OF H(RHO,TH) IN B
RHO 0.5 1.0 1.5 2.0 2.5 3.0
TH
20 .11237168, 00 -.79607559, 01 -.22410068, 03 -.22692020, 04 -.13574638, 05
25 .11237168, 00 -.79607559, 01 -.22410068, 03 -.22692020, 04 -.13574638, 05
30 .11237168, 00 -.79607559, 01 -.22410068, 03 -.22692020, 04 -.13574638, 05
35 .11237168, 00 -.79607559, 01 -.22410068, 03 -.22692020, 04 -.13574638, 05
40 .11237168, 00 -.79607559, 01 -.22410068, 03 -.22692020, 04 -.13574638, 05
45 .11237168, 00 -.79607559, 01 -.22410068, 03 -.22692020, 04 -.13574638, 05
Reference in New Issue View Git Blame Copy Permalink