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

Commit transcription corrections #1 to DIRICHLET example.

Browse Source
This commit is contained in:
Paul Kimpel
2018-02-24 07:13:57 -08:00
parent 714993f651
commit 68b293813d
1 changed files with 186 additions and 189 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -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