Приближенное вычисление определенного интеграла

[Novichok2016]

Задание: Функция f(x)=x*x*arcctg(x/2)/x, отрезок a=1, b=2, точность E=0,005. Использовать квадратурную формулу трапеция, язык программирования С++, многопоточное программирование.

Вопрос: Подскажите,пожалуйста, что я делаю неправильно, почему результат уходит в бесконечность?

#include "stdafx.h"
#include <iostream>
#include <thread>
#include <ctime>
#include <fstream>

using namespace std;

double f(double x)
{
	return (x*x*(atan(-x / 2) + 3.14 / 2) / x);
}

bool func1(double a, double b, double e, double *func)
{
	double f1, f2, fsr, s1, s2, sl;
	double sr = (a + b) / 2;
	double h = b - a;
	f1 = f(a);
	f2 = f(b);
	fsr = f(sr);
	s1 = 0.5*h*(f1 + f2);
	sl = 0.25*h*(f1 + fsr);
	sr = 0.25*h*(fsr + f2);
	s2 = sl + sr;
	*func = s2;
	if (fabs(s1 - s2)<e)
		return true;
	else
		return false;
}

double func2(double a, double b, double e, int T)
{
	if (T == 1)
	{
		double fx;
		if (func1(a, b, e, &fx))
			return fx;
		double m = (a + b) / 2;
		double e1 = e / 2;
		double a1, b1;
		a1 = func2(a, m, e1, T);
		b1 = func2(m, b, e1, T);
		return a1 + b1;
	}
	else
	{
		T = T / 2;
		thread func_thread1(func2, a, b, e, T);
		thread func_thread2(func2, a, b, e, T);
		if (func_thread1.joinable())
			func_thread1.join();
		if (func_thread2.joinable())
			func_thread2.join();
	}
}

int main(int argc, char *argv[])
{
	setlocale(0, "");
	srand(time(0));
	int T, n = 100;
	double a = 1, b = 2, e = 0.005, I = 0.0, h;
	cout << "Введите кол-во потоков: ";
	cin >> T;
	h = (b - a) / n;
	//n--;

	I = func2(a, b, e, T);
	cout << "I= " << I << endl;

	cout << endl << endl;
	system("pause");
	return 0;
}

[type_Oleg]

А что за странная функция f(x)=x*x*arcctg(x/2)/x ?
Почему нельзя сократить x ?

]f(x) = x*x*arcctg(x/2)/x = x*arcctg(x/2)

[Smitt&Wesson]

Можно. Но это-ж такая фигня , получится...

[Novichok2016]

Можно сократить, но ничего не изменится - интеграл не вычисляется...