#include <stdio.h>
double f(double x)
{
return 4.0 / (1.0 + x * x);
}
int main()
{
int k, n, i;
double d, x, h, S;
for (k = 1; k <= 10; k++) {
n = 1 << k;
d = 0.5 * f(0.0);
for (i = 1; i < n; i++) {
x = (double)i / n;
d += f(x);
}
d += 0.5 * f(1.0);
h = 1.0 / n;
S = h * d;
}
return 0;
}
I2luY2x1ZGUgPHN0ZGlvLmg+Cgpkb3VibGUgZihkb3VibGUgeCkKewoJcmV0dXJuIDQuMCAvICgxLjAgKyB4ICogeCk7Cn0KCmludCBtYWluKCkKewoJaW50IGssIG4sIGk7Cglkb3VibGUgZCwgeCwgaCwgUzsKCglmb3IgKGsgPSAxOyBrIDw9IDEwOyBrKyspIHsKCQluID0gMSA8PCBrOwoJCWQgPSAwLjUgKiBmKDAuMCk7CgkJZm9yIChpID0gMTsgaSA8IG47IGkrKykgewoJCQl4ID0gKGRvdWJsZSlpIC8gbjsKCQkJZCArPSBmKHgpOwoJCX0KCQlkICs9IDAuNSAqIGYoMS4wKTsKCQloID0gMS4wIC8gbjsKCQlTID0gaCAqIGQ7CgkJcHJpbnRmKCJrPSVkIFM9JWZcbiIsIGssIFMpOwoJfQoJcmV0dXJuIDA7Cn0K