diff --git a/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card b/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card index edd4c31..24c0d3b 100644 --- a/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card +++ b/software/BALGOL/BALGOL-Examples/DIRICHLET/DIRICHLET.card @@ -2,41 +2,40 @@ % DIRICHLET PROBLEM FOR A BEAN-SHAPED REGION. FROM P J DAVIS, % "ORTHONORMALIZING CODES IN NUMERICAL ANALYSIS" IN J TODD, --SURVEY OF % NUMERICAL ANLAYSIS--, MCGRAW-HILL, 1962, P.347. P H KIMPEL 8/15/70 - % MODIFICATION LOG: + % MODIFICATION LOG.. % 92/03/17 P.KIMPEL, PARADIGM CORP, SAN DIEGO, CA 92117. % CONVERT FOR UNISYS A-SERIES MCP 3.8.4. % 2014-11-15 P.KIMPEL % RETRO-CONVERT FROM UNISYS MCP ALGOL BACK TO BURROUGHS B5500 XALGOL. % 2018-02-20 P.KIMPEL - % RETRO-CONVERT FROM B5500 TO BURROUGHS 220 BALGOL. + % RETRO-CONVERT FROM B5500 TO BURROUGHS 220 BALGOL OF 1962. ; - FORMAT F1 (B30,I3,2X10.3,4X10.5), - F2 (B40,*X*,B9,*Y*,B9,*W*,B8,*BV*,B7,*CBV,*,B7,*DEV*), - FT (W4,B40,*PROCESSOR TIME: *, F6.2,* SEC*); - INTEGER I, J; + FORMAT F1 (B30,I3,2X10.3,4X10.5,W), + F2 (B40,*X*,B9,*Y*,B9,*W*,B8,*BV*,B7,*CBV,*,B7,*DEV*,W3,W); + INTEGER I, J, N, M, P; REAL SUM, GMDT; COMMENT ARRAY DIMENSIONS WERE ORIGINALLY DEFINED AS N=43, M=11, P=11.. X, % ABSCISSA VALUES. Y, % ORDINATE VALUES. W, % WEIGHTS. - CBV[0:N1], % BOUNDARY VALUES CALCULATED FROM ORTHO COEFS. - BV[0:0, 0:N1+P], % GIVEN BOUNDARY VALUES. - Z[0:M1, 0:N1+P], % APPROXIMATING VECTORS. - ORTHV[0:M, 0:N1+P], % ORTHONORMAL VECTORS RETURNED BY "ORTHO" - DEV[0:0, 0:N1], % DEVIATIONS. - COF[0:0, 0:P1], % COEFFICIENTS. - STD[0:0], % STANDARD DEVIATION. - CV[0:P, 0:P1], % COVARIANCE MATRIX. - VCV[0:0, 0:P, 0:P1], % VARIANCE/COVARIANCE MATRIX. - Q[0:0, 0:M], % FOURIER COEFFICIENTS. + CBV(0..N-1), % BOUNDARY VALUES CALCULATED FROM ORTHO COEFS + BV(0..0, 0..N-1+P), % GIVEN BOUNDARY VALUES. + Z(0..M-1, 0..N-1+P), % APPROXIMATING VECTORS. + ORTHV(0..M, 0..N-1+P), % ORTHONORMAL VECTORS RETURNED BY "ORTHO" + DEV(0..0, 0..N-1), % DEVIATIONS. + COF(0..0, 0..P-1), % COEFFICIENTS. + STD(0..0), % STANDARD DEVIATION. + CV(0..P, 0..P-1), % COVARIANCE MATRIX. + VCV(0..0, 0..P, 0..P-1),% VARIANCE/COVARIANCE MATRIX. + Q(0..0, 0..M), % FOURIER COEFFICIENTS. Q2, % SQUARED FOURIER COEFFICIENTS. E, % SUM OF SQUARED RESIDUALS. - EP[0:0, 0:M1], % RESIDUALS. - A[0:M1, 0:P1], % LOWER TRIANGULAR MATRIX USED TO CALC CV. - GF[0:M], % GRAM FACTORS. - ENF[0:M1]; % NORMS OF THE APPROXIMATING VECTORS. + EP(0..0, 0..M-1), % RESIDUALS. + A(0..M-1, 0..P-1), % LOWER TRIANGULAR MATRIX USED TO CALC CV. + GF(0..M), % GRAM FACTORS. + ENF(0..M-1); % NORMS OF THE APPROXIMATING VECTORS. ARRAY X(43), @@ -67,102 +66,104 @@ REAL FN, GMDT; INTEGER N, M, P, R, AI, AUI, ZEI, MUI; COMMENT - ORTHO IS TAKEN FROM ACM ALGORITHM 127 [COMM. ACM, VOL.5, - OCTOBER 1962, P. 511, AUTHOR PHILIP J. WALSH]; + ORTHO IS TAKEN FROM ACM ALGORITHM 127 (COMM. ACM, VOL.5, + OCTOBER 1962, P. 511, AUTHOR PHILIP J. WALSH); INTEGER NPP, NPM, M1, N2, M2, R1, RBAR, P2, BEI, RHI, I18, GAI, SII, I, J, DEI, NUI, E1Z2, E1Z1, K, THI, ALI, OMI, NII; - ARRAY PK, XP [N+P], QK[M+1]; + ARRAY PK, XP (N+P), QK(M+1); REAL DENOM, SUM, DK2, DK, FI, SS, SSQ; NPP = N+P; NPM = N+M; M1 = M-1; N2 = N+1; M2 = M+1; R1 = 0; RBAR = R; P2 = P+1; - EITHER IF N EQL M; DENOM = 1.0; OTHERWISE; DENOM = SQRT(N-M); + EITHER IF N EQL M; DENOM = 1.0; OTHERWISE; DENOM = SQRT(N-M); BEI = RHI = I18 = 1; - EITHER IF (P NEQ 0); GAI = SII = 2 OTHERWISE; GAI = SII = 1; + EITHER IF (P NEQ 0); GAI = SII = 2; OTHERWISE; GAI = SII = 1; BOX1.. SWITCH AI, (AT1, AT2); AT1.. FOR J = (1, 1, N); BEGIN - X[2,J] = Z[1,J]; X[1,J] = 1.0 END; + X(2,J) = Z(1,J); X(1,J) = 1.0 END; FOR I = (2, 1, M1); BEGIN FOR J = (1, 1, N); - X[I+1,J] = X[I,J] . X[2,J] END; GO TO BOX2; + X(I+1,J) = X(I,J) . X(2,J) END; GO TO BOX2; AT2.. FOR I = (1, 1, M); BEGIN FOR J = (1, 1, N); - X[I,J] = Z[I,J] END; - BOX2.. EITHER IF P EQL 0; GO TO BOX3 OTHERWISE; SWITCH AUI, (AU1, AU2); + X(I,J) = Z(I,J) END; + BOX2.. EITHER IF P EQL 0; GO TO BOX3; + OTHERWISE; SWITCH AUI, (AU1, AU2); AU1.. FOR I = (1, 1, M); BEGIN FOR J = (N2, 1, NPP); - X[I,J] = 0.0; X[I,N+I] = 1.0 END; GO TO BOX3; + X(I,J) = 0.0; X(I,N+I) = 1.0 END; GO TO BOX3; AU2.. FOR I = (1, 1, M); BEGIN FOR J = (N2, 1, NPP); - X[I,J] = Z[I,J] END; + X(I,J) = Z(I,J) END; BOX3.. DEI = NUI = E1Z1 = E1Z2 = K = 1; BOX4.. THI = 1; - BOX5.. ALI = OMI = 1; EITHER IF P EQL 0; GO TO BOX6 OTHERWISE; - FOR J = (1, 1, P); PK[N+J] = 0.0; + BOX5.. ALI = OMI = 1; EITHER IF P EQL 0; GO TO BOX6; + OTHERWISE; FOR J = (1, 1, P); PK(N+J) = 0.0; BOX6.. GO TO SWITCH MUI, (MU1, MU2); - MU1.. FOR I = (1, 1, N); PK[I] = X[K,I]; + MU1.. FOR I = (1, 1, N); PK(I) = X(K,I); GO TO BOX7; MU2.. FOR I = (1, 1, N); - PK[I] = X[K,I] . W[I]; GO TO BOX7; + PK(I) = X(K,I) . W(I); GO TO BOX7; BOX7.. SWITCH OMI, (OM1, OM2); OM1.. FOR I = (1, 1, K); BEGIN SUM = 0.0; FOR J = (1, 1, NPP); - SUM = SUM + PK[J] . X[I,J]; QK[I] = SUM END; + SUM = SUM + PK(J) . X(I,J); QK(I) = SUM END; GO TO BOX8; OM2.. DK2 = 0.0; FOR I = (1, 1, NPP); - DK2 = DK2 + PK[I] . X[K,I]; + DK2 = DK2 + PK(I) . X(K,I); DK = SQRT(DK2); - GF[I18] = DK; I18 = I18 + 1; + GF(I18) = DK; I18 = I18 + 1; FOR I = (1, 1, NPP); - X[K,I] = X[K,I]/DK;; + X(K,I) = X(K,I)/DK;; OMI = 1; GO TO BOX6; BOX8.. SWITCH DEI, (DE1, DE2); - DE1.. E1Z1 = -E1Z1; EITHER IF E1Z1 < 0; GO TO BOX8B OTHERWISE; - GO TO BOX8A; + DE1.. E1Z1 = -E1Z1; EITHER IF E1Z1 < 0; GO TO BOX8B; + OTHERWISE; GO TO BOX8A; BOX8A.. FOR I = (1, 1, K-1); - QK[I] = -QK[I]; QK[K] = 1.0; + QK(I) = -QK(I); QK(K) = 1.0; FOR I = (1, 1, NPP); BEGIN SUM = 0.0; FOR J = (1, 1, K); - SUM = SUM + X[J,I] . QK[J]; - XP[I] = SUM END; GO TO BOX9; - BOX8B.. ENF[I18] = SQRT(QK[K]); GO TO BOX8A; - DE2.. E1Z2 = -E1Z2; EITHER IF E1Z2 < 0; GO TO BOX8C OTHERWISE; - GO TO BOX8A; + SUM = SUM + X(J,I) . QK(J); + XP(I) = SUM END; GO TO BOX9; + BOX8B.. ENF(I18) = SQRT(QK(K)); GO TO BOX8A; + DE2.. E1Z2 = -E1Z2; EITHER IF E1Z2 < 0; GO TO BOX8C; + OTHERWISE; GO TO BOX8A; BOX8C.. FOR I = (1, 1, M); BEGIN - Q[R1,I] = QK[I]; Q2[R1,I] = QK[I] . QK[I] END; - Q[R1,M2] = QK[M2]; E[R1,1] = Q[R1,M2] - Q2[R1,1]; + Q(R1,I) = QK(I); Q2(R1,I) = QK(I) . QK(I) END; + Q(R1,M2) = QK(M2); E(R1,1) = Q(R1,M2) - Q2(R1,1); FOR J = (2, 1, M); - E[R1,J] = E[R1,J-1] - Q2[R1,J]; + E(R1,J) = E(R1,J-1) - Q2(R1,J); FI = 1.0; FOR I = (1, 1, M); BEGIN - EITHER IF (FN - FI) > 0.0; BEGIN EITHER IF E[R1,I] < 0.0; - BEGIN EP[R1,I] = -SQRT(ABS(E[R1,I])/(FN - FI)); + EITHER IF (FN - FI) > 0.0; BEGIN + EITHER IF E(R1,I) < 0.0; + BEGIN EP(R1,I) = -SQRT(ABS(E(R1,I))/(FN - FI)); GO TO BOX8D; END - OTHERWISE; EP[R1,I] = SQRT(E[R1,I]/(FN - FI)); - GO TO BOX8D; END OTHERWISE; E[R1,I] = -1.0; + OTHERWISE; EP(R1,I) = SQRT(E(R1,I)/(FN - FI)); + GO TO BOX8D; END; OTHERWISE; E(R1,I) = -1.0; BOX8D.. FI = FI + 1.0; END; GO TO BOX8A; BOX9.. SWITCH THI, (TH1, TH2, TH3); TH1.. FOR I = (1, 1, NPP); - X[K,I] = XP[I]; GO TO BOX10; + X(K,I) = XP(I); GO TO BOX10; TH2.. FOR I = (1, 1, N); - DEV[R1,I] = XP[I]; + DEV(R1,I) = XP(I); FOR I = (1, 1, P); - COF[R1,I] = -XP[N+I]; THI = 3; GO TO TH1; + COF(R1,I) = -XP(N+I); THI = 3; GO TO TH1; TH3.. GO TO BOX11; BOX10.. SWITCH ALI, (AL1, AL2); AL1.. OMI = ALI = 2; GO TO BOX6; - AL2.. EITHER IF K < M; BEGIN K = K + 1; GO TO BOX4; END + AL2.. EITHER IF K < M; BEGIN K = K + 1; GO TO BOX4; END; OTHERWISE; GO TO BOX12; BOX11.. SWITCH NUI, (NU1, NU2); NU1.. NUI = 2; GO TO BOX14; NU2.. SS = DK/DENOM; SSQ = SS . SS; - STD[R1] = SS; GO TO BOX14; + STD(R1) = SS; GO TO BOX14; BOX12.. SWITCH BEI, (BE1, BE2); BE1.. FOR I = (1, 1, M); BEGIN FOR J = (1, 1, P); - A[I,J] = X[I,N+J] END; + A(I,J) = X(I,N+J) END; GMDT = 1.0; FOR I = (1, 1, M); - GMDT = GMDT . (GF[I]/ENF[I]); + GMDT = GMDT . (GF(I)/ENF(I)); GMDT = GMDT . GMDT; DEI = BEI = THI = 2; K = K + 1; GO TO BOX13; BE2.. GO TO BOX11; @@ -172,27 +173,28 @@ FOR J = (I, 1, P); BEGIN SUM = 0.0; FOR NII = (1, 1, M); - SUM = SUM + A[NII,I] . A[NII,J]; - CV[I,J] = SUM END END; + SUM = SUM + A(NII,I) . A(NII,J); + CV(I,J) = SUM END END; FOR I = (1, 1, P); - CV[P2,I] = SQRT(CV[I,I]); GAI = 1; GO TO BOX11; + CV(P2,I) = SQRT(CV(I,I)); GAI = 1; GO TO BOX11; BOX14.. SWITCH RHI, (RH1, RH2); - RH1.. EITHER IF RBAR EQL 0; GO TO FINAL OTHERWISE; RBAR = RBAR - 1; + RH1.. EITHER IF RBAR EQL 0; GO TO FINAL; + OTHERWISE; RBAR = RBAR - 1; R1 = R1 + 1; THI = RHI = 2; SWITCH ZEI, (ZE1, ZE2); ZE1.. FOR I = (1, 1, N); - X[M2,I] = Y[R1,I]; + X(M2,I) = Y(R1,I); FOR I = (1, 1, P); - X[M2,N+I] = 0.0; GO TO BOX5; + X(M2,N+I) = 0.0; GO TO BOX5; ZE2.. FOR I = (1, 1, NPP); - X[M2,I] = Y[R1,I]; GO TO BOX5; + X(M2,I) = Y(R1,I); GO TO BOX5; RH2.. SWITCH SII, (SI1, SI2); SI1.. GO TO RH1; SI2.. FOR I = (1, 1, P); BEGIN FOR J = (I, 1, P); - VCV[R1,I,J] = SSQ . CV[I,J] END; + VCV(R1,I,J) = SSQ . CV(I,J) END; FOR I = (1, 1, P); - VCV[R1, P2, I] = SS . CV[P2,I]; GO TO RH1; - FINAL.. RETURN END ORTHO ; + VCV(R1, P2, I) = SS . CV(P2,I); GO TO RH1; + FINAL.. RETURN END ORTHO(); PROCEDURE G (I, X, Y); REAL G; @@ -256,134 +258,129 @@ G = 8.0 . X*7 . Y - 56.0 . X*5 . Y*3 + 56.0 . X*3 . Y*5 - 8.0 . X . Y*7; RETURN; - END G; + END G(); - FOR I = (0, 1, N1); + COMMENT INITIALIZATION; + N = 43; + M = 11; + P = 11; + FOR I = (1, 1, N); BEGIN - READ (CDS, /, X[I], Y[I], W[I]); - BV[0,I] = EXP(X[I]) . COS(Y[I]) + LOG((1 - Y[I])*2 + X[I]*2); - FOR J = (0, 1, M1); - Z[J,I] = G(J+1, X[I], Y[I]); + READ (;; XYWIN); + INPUT XYWIN (X(I), Y(I), W(I)); + BV(1,I) = EXP(X(I)) . COS(Y(I)) + LOG((1 - Y(I))*2 + X(I)*2); + FOR J = (1, 1, M1); + Z(J,I) = G(J+1, X(I), Y(I)); END; - CLOSE (CDS); - ORTHO (W, BV, Z, N , N , M , P, 1, 2, 1, 2, 1, ORTHV, DEV, COF, STD, + + ORTHO (W, BV, Z, N , N , M , P, 1, 2, 1, 2, 1; ORTHV, DEV, COF, STD, CV, VCV, GMDT, Q, Q2, E, EP, A, GF, ENF); - FOR I = (0, 1, N1); + + FOR I = (1, 1, N); BEGIN SUM = 0; - FOR J = (0, 1, M1); SUM = SUM + COF[0,J].G(J+1, X[I], Y[I] - ); - CBV[I] = SUM; + FOR J = (1, 1, M); SUM = SUM + COF(1,J).G(J+1, X(I), Y(I)); + CBV(I) = SUM; END; - WRITE (PR[DBL], F2); - WRITE (PR, F1, FOR I = (0, 1, N1); [I,X[I],Y[I],W[I], BV[0,I], - CBV[I], (CBV[I]-BV[0,I])]); - WRITE (PR, FT, TIME(2)/60); - QUIT.. + WRITE (;; F2); + WRITE (;; RESULTS, F1); + OUTPUT RESULTS (FOR I = (1, 1, N); (I,X(I),Y(I),W(I), BV(1,I), + CBV(I), (CBV(I)-BV(1,I)))); - PROCEDURE DMMP (NAME, ROW, SZ); - VALUE NAME, SZ; - ALPHA NAME, SZ; - ARRAY ROW[0]; + PROCEDURE DMMP (NAME, ROW(), RN, SZ); BEGIN - REAL - I, - UB; - FORMAT - F (A6," = ",/*(6E20.11,/)); - UB = SZ-1; - WRITE (PR[DBL]); - WRITE (PR[DBL]); - WRITE (PR, F, NAME, (SZ+5)DIV 6, - FOR I = (0, 1, UB); ROW[I]); - END DMMP; + INTEGER + NAME, SZ; + INTEGER + I; + FORMAT + F1 (W4,A5,* = *,W), + F2 (W4,A5,*(*,I2,*) = *,W), + F3 (6F20.11,W); + OUTPUT + D1OUT (NAME), + S2OUT (NAME, RN), + ROWOUT (FOR I = (1, 1, SZ); ROW(I)); - ARRAY NR[0:99]; - FOR I = (0, 1, 99); - REPLACE POINTER(NR[I])+6 BY I FOR 2 DIGITS; + EITHER IF SZ GTR 0; + WRITE (;; D2OUT, F2); + OTHERWISE; + WRITE (;; D1OUT, F1); - DMMP ("X ", X, N1+1); - DMMP ("Y ", Y, N1+1); - DMMP ("W ", W, N1+1); - DMMP ("CBV ", CBV, N1+1); - DMMP ("BV ", BV[0,*], N1+P+1); - FOR I = (0, 1, M1); - BEGIN - SUM = "Z " & NR[I] [11:12]; - DMMP (SUM, Z[I,*], N1+P+1); - END; - FOR I = (0, 1, M); - BEGIN - SUM = "ORTHV " & NR[I] [11:12]; - DMMP (SUM, ORTHV[I,*], N1+P+1); - END; - DMMP ("DEV ", DEV[0,*], N1+1); - DMMP ("COF ", COF[0,*], P1+1); - DMMP ("STD ", STD, 1); - FOR I = (0, 1, P); - BEGIN - SUM = "CV " & NR[I] [11:12]; - DMMP (SUM, CV[I,*], P1+1); - END; - FOR I = (0, 1, P); - BEGIN - SUM = "VCV " & NR[I] [11:12]; - DMMP (SUM, VCV[0,I,*], P1+1); - END; - DMMP ("EP ", EP[0,*], M1+1); - FOR I = (0, 1, M1); - BEGIN - SUM = "A " & NR[I] [11:12]; - DMMP (SUM, A[I,*], P1+1); - END; - DMMP ("GF ", GF, M+1); - DMMP ("Q ", Q[0,*], M+1); - DMMP ("Q2 ", Q2[0,*], M1+1); - DMMP ("E ", E[0,*], M1+1); - DMMP ("ENF ", ENF, M1+1); + WRITE (;; ROWOUT, F3); + END DMMP(); + DMMP ("X ", X, 0, N); + DMMP ("Y ", Y, 0, N); + DMMP ("W ", W, 0, N); + DMMP ("CBV ", CBV, 0, N); + DMMP ("BV ", BV(1,), 0, N+P); + FOR I = (1, 1, M); + DMMP ("Z ", Z(I,), I, N+P); + + FOR I = (1, 1, M); + DMMP ("ORTHV", ORTHV(I,), I, N+P); + + DMMP ("DEV ", DEV(1,), 0, N); + DMMP ("COF ", COF(1,), 0, P); + DMMP ("STD ", STD, 0, 1); + FOR I = (1, 1, P); + DMMP ("CV ", CV(I,), I, P); + + FOR I = (1, 1, P); + DMMP ("VCV ", VCV(1,I,), I, P); + + DMMP ("EP ", EP(1,), M); + FOR I = (1, 1, M); + DMMP ("A ", A(I,), I, P); + + DMMP ("GF ", GF, 0, M+1); + DMMP ("Q ", Q(1,), 0, M+1); + DMMP ("Q2 ", Q2(1,), 0, M); + DMMP ("E ", E(1,), 0, M); + DMMP ("ENF ", ENF, 0, M); FINISH; -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 +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