24 lines
520 B
Text
24 lines
520 B
Text
def MaxIt = 13, MaxItJ = 10;
|
|
real A, A1, A2, D, D1, X, Y;
|
|
int I, J, K;
|
|
[A1:= 1.; A2:= 0.; D1:= 3.2;
|
|
Text(0, " i d^m^j");
|
|
for I:= 2 to MaxIt do
|
|
[A:= A1 + (A1-A2)/D1;
|
|
for J:= 1 to MaxItJ do
|
|
[X:= 0.; Y:= 0.;
|
|
for K:= 1 to 1<<I do
|
|
[Y:= 1. - 2.*Y*X;
|
|
X:= A - X*X;
|
|
];
|
|
A:= A - X/Y;
|
|
];
|
|
D:= (A1-A2) / (A-A1);
|
|
Format(2, 0); RlOut(0, float(I));
|
|
Format(5, 8); RlOut(0, D);
|
|
CrLf(0);
|
|
D1:= D;
|
|
A2:= A1;
|
|
A1:= A;
|
|
];
|
|
]
|