29 lines
648 B
D
29 lines
648 B
D
import std.stdio;
|
|
|
|
void main() {
|
|
int max_it = 13;
|
|
int max_it_j = 10;
|
|
double a1 = 1.0;
|
|
double a2 = 0.0;
|
|
double d1 = 3.2;
|
|
double a;
|
|
|
|
writeln(" i d");
|
|
for (int i=2; i<=max_it; i++) {
|
|
a = a1 + (a1 - a2) / d1;
|
|
for (int j=1; j<=max_it_j; j++) {
|
|
double x = 0.0;
|
|
double y = 0.0;
|
|
for (int k=1; k <= 1<<i; k++) {
|
|
y = 1.0 - 2.0 * y * x;
|
|
x = a - x * x;
|
|
}
|
|
a -= x / y;
|
|
}
|
|
double d = (a1 - a2) / (a - a1);
|
|
writefln("%2d %.8f", i, d);
|
|
d1 = d;
|
|
a2 = a1;
|
|
a1 = a;
|
|
}
|
|
}
|