#include <stdio.h>
#include <conio.h>
#include <math.h>

double f(double x)
{
	double y;
	y = pow((1 + x), (-1. / 3.));
	return y;
}

void main()
{
	double ii = 1.03730716743448; //истинное значение интеграла
	double a = 0, b = 1.2, h, n;
	double it1, it2, it0, itn, it2n;

	h = b / 3;
	it0 = ((f(a) + f(b))*h) / 2;
	it1 = 0;
	for (n = h; n <= (b - h); n += h)
	{
		itn = h*f(n);
		it1 += itn;
	}
	it1 += it0;
	it2 = it1;
	while (fabs(it2 - ii) > 0.00001)
	{
		it1 = it2;
		it2n = 0;
		h = h / 2;
		for (n = h; n <= (b - h); n += 2 * h)
		{
			itn = h*f(n);
			it2n += itn;
		}
		it2 = (it1 / 2) + it2n;
	}
	printf("%.8lf", it2);
	_getch();
}