import std.cstream;
import std.math;
void main() {
int N = 10;
real[] a = new real[N], b = new real[N], t = new real[N], pi = new real[N];
a[0] = 1.0;
b[0] = 1.0 / sqrt(2.0);
t[0] = 1.0 / 4.0;
for(int i=0; i<N-1; i++) {
a[i+1] = (a[i] + b[i]) / 2.0;
b[i+1] = sqrt(a[i] * b[i]);
t[i+1] = t[i] - pow(2.0, cast(real)i) * pow(a[i] - a[i+1], 2.0);
}
for(int i=0; i<N; i++) {
dout.writefln("%40.37f", pow(a[i] + b[i], 2.0) / (4.0 * t[i]));
}
dout.writefln("%40.37f", real.epsilon);
}