github.com/pkimpel.retro-220
1
0
Fork 0
mirror of https://github.com/pkimpel/retro-220.git synced 2026-04-10 14:49:05 +00:00
Code Issues Releases Wiki Activity

Commit transcription corrections #3 to DIRICHLET example.

Browse Source
This commit is contained in:
Paul Kimpel
2018-02-24 07:15:35 -08:00
parent fc08dd6bed
commit 3ce65f8c86
1 changed files with 57 additions and 55 deletions
Download Patch File Download Diff File Expand all files Collapse all files

112
software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card
View File

@@ -42,20 +42,20 @@
2 Y(43),
2 W(43),
2 CBV(43),
2 BV(54),
2 Z(11, 54),
2 ORTHV(11, 54),
2 BV(43+11),
2 Z(11, 43+11),
2 ORTHV(11, 43+11),
2 DEV(1, 43),
2 COF(1, 11),
2 STD(1),
2 CV(12, 11),
2 VCV(1, 12, 11),
2 Q(1, 12),
2 CV(11+1, 11),
2 VCV(1, 11+1, 11),
2 Q(1, 11+1),
2 Q2(1, 11),
2 E(1, 11),
2 EP(1, 11),
2 A(11, 11),
2 GF(12),
2 GF(11+1),
2 ENF(11);
2
2
@@ -64,13 +64,14 @@
2 Q(,), Q2(,), E(,), EP(,), A(,), GF(), ENF());
2BEGIN
2REAL FN, GMDT;
2INTEGER N, M, P, R, AI, AUI, ZEI, MUI;
2REAL W, Y, Z, X, DEV, COF, STD, CV, VCV, Q, Q2, E, EP, A, GF, ENF;
2INTEGER N, M, P, R, AI, AUI, ZEI, MUI;
2COMMENT
2 ORTHO IS TAKEN FROM ACM ALGORITHM 127 (COMM. ACM, VOL.5,
2 OCTOBER 1962, P. 511, AUTHOR PHILIP J. WALSH);
2INTEGER NPP, NPM, M1, N2, M2, R1, RBAR, P2, BEI, RHI, I18, GAI, SII, I,
2 J, DEI, NUI, E1Z2, E1Z1, K, THI, ALI, OMI, NII;
2ARRAY PK(N+P), XP (N+P), QK(M+1);
2INTEGER NPP, NPM, M1, N2, M2, R1, RBAR, P2, BEI, RHI, I18, GAI, SII,
2 I, J, DEI, NUI, E1Z2, E1Z1, K, THI, ALI, OMI, NII;
2ARRAY PK(99), XP (99), QK(99);
2REAL DENOM, SUM, DK2, DK, FI, SS, SSQ;
2
2 NPP = N+P; NPM = N+M; M1 = M-1; N2 = N+1; M2 = M+1;
@@ -114,7 +115,7 @@
2 DK = SQRT(DK2);
2 GF(I18) = DK; I18 = I18 + 1;
2 FOR I = (1, 1, NPP);
2 X(K,I) = X(K,I)/DK;;
2 X(K,I) = X(K,I)/DK;
2 OMI = 1; GO TO BOX6;
2BOX8.. SWITCH DEI, (DE1, DE2);
2 DE1.. E1Z1 = -E1Z1; EITHER IF E1Z1 LSS 0; GO TO BOX8B;
@@ -198,6 +199,7 @@
2
2
2PROCEDURE G (I, X, Y; V);
2 BEGIN
2 REAL X, Y, V; INTEGER I;
2 SWITCH I, (G1, G2, G3, G4, G5, G6, G7, G8, G9, G10,
2 G11, G12, G13, G14, G15, G16, G17);
@@ -343,47 +345,47 @@
2DMMP ("E ", E(1,), 0, M);
2DMMP ("ENF ", ENF, 0, M);
2FINISH;
5 0.000, 0.110, 0.01414,
5 -0.050, 0.108, 0.01427,
5 -0.100, 0.115, 0.01963,
5 -0.160, 0.150, 0.02300,
5 -0.220, 0.205, 0.03897,
5 -0.320, 0.300, 0.02792,
5 -0.400, 0.358, 0.03324,
5 -0.500, 0.420, 0.01483,
5 -0.550, 0.436, 0.01423,
5 -0.600, 0.430, 0.01505,
5 -0.644, 0.400, 0.01483,
5 -0.660, 0.350, 0.01420,
5 -0.655, 0.300, 0.02881,
5 -0.635, 0.200, 0.03043,
5 -0.595, 0.100, 0.03076,
5 -0.552, 0.000, 0.03311,
5 -0.500, -0.105, 0.03175,
5 -0.440, -0.200, 0.01809,
5 -0.400, -0.250, 0.01998,
5 -0.350, -0.300, 0.01882,
5 -0.300, -0.344, 0.03140,
5 -0.204, -0.400, 0.03450,
5 -0.100, -0.436, 0.02846,
5 0.000, -0.448, 0.02831,
5 0.100, -0.442, 0.03860,
5 0.230, -0.400, 0.02431,
5 0.300, -0.350, 0.02059,
5 0.353, -0.300, 0.03566,
5 0.430, -0.200, 0.03122,
5 0.477, -0.100, 0.02975,
5 0.510, 0.000, 0.02846,
5 0.522, 0.100, 0.01696,
5 0.520, 0.160, 0.02330,
5 0.500, 0.240, 0.02102,
5 0.456, 0.300, 0.01795,
5 0.400, 0.330, 0.01147,
5 0.360, 0.337, 0.01762,
5 0.300, 0.320, 0.01648,
5 0.250, 0.290, 0.01901,
5 0.300, 0.245, 0.01901,
5 0.150, 0.200, 0.01809,
5 0.100, 0.160, 0.01677,
5 0.050, 0.128, 0.01501,
5 0.000 0.110 0.01414
5 -0.050 0.108 0.01427
5 -0.100 0.115 0.01963
5 -0.160 0.150 0.02300
5 -0.220 0.205 0.03897
5 -0.320 0.300 0.02792
5 -0.400 0.358 0.03324
5 -0.500 0.420 0.01483
5 -0.550 0.436 0.01423
5 -0.600 0.430 0.01505
5 -0.644 0.400 0.01483
5 -0.660 0.350 0.01420
5 -0.655 0.300 0.02881
5 -0.635 0.200 0.03043
5 -0.595 0.100 0.03076
5 -0.552 0.000 0.03311
5 -0.500 -0.105 0.03175
5 -0.440 -0.200 0.01809
5 -0.400 -0.250 0.01998
5 -0.350 -0.300 0.01882
5 -0.300 -0.344 0.03140
5 -0.204 -0.400 0.03450
5 -0.100 -0.436 0.02846
5 0.000 -0.448 0.02831
5 0.100 -0.442 0.03860
5 0.230 -0.400 0.02431
5 0.300 -0.350 0.02059
5 0.353 -0.300 0.03566
5 0.430 -0.200 0.03122
5 0.477 -0.100 0.02975
5 0.510 0.000 0.02846
5 0.522 0.100 0.01696
5 0.520 0.160 0.02330
5 0.500 0.240 0.02102
5 0.456 0.300 0.01795
5 0.400 0.330 0.01147
5 0.360 0.337 0.01762
5 0.300 0.320 0.01648
5 0.250 0.290 0.01901
5 0.300 0.245 0.01901
5 0.150 0.200 0.01809
5 0.100 0.160 0.01677
5 0.050 0.128 0.01501
5 SENTINEL