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