double squ(double x, double bb)
{
    return 2 * x * func(bb - pow(x , 2));
}

void midsqu(double aa, double bb, double& s, int n)
{
	int it,tnm;
	double del,ddel,x,sum;
	double b,a;
    b = sqrt(bb - aa);
    a = 0.0;
    if (n == 1)
	{
        s = (b - a) * squ(0.5 * (a + b), bb);
        it = 1;
	}
    else
	{
        it = (int)pow(3, n - 2);
        tnm = it;
        del = (b - a) / (3.0 * tnm);
        ddel = del + del;
        x = a + 0.5 * del;
        sum = 0.0;
        for (int j = 1; j<=it; j++)
		{
            sum = sum + squ(x, bb);
            x = x + ddel;
            sum = sum + squ(x, bb);
            x = x + del;
        }
        s = (s + (b - a) * sum / tnm) / 3.0;
    }
}