diff --git a/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card b/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card index b1653be..7b881f6 100644 --- a/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card +++ b/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card @@ -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