#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;
		printf("k=%d S=%f\n", k, S);
	}
	return 0;
}
