Проинтегрировать функцию f(x) = x^2 * arctg(x)

#include 

#include 

#include 

#include 

void main()

{

        using namespace std;

        setlocale(LC_ALL, "Rus");

        system("cls");

        int N, i;

        float a, b, h, res = 0, f1, fx, x, f2, f3 = 0;

        cout << "Проинтегрировать функцию f(x) = x^2 * arctg(x) : \n";

        cout << "Введите a,b и количество шагов N:";

        cin >> a >> b >> N; //вводим а, b и N

        h = (b - a) / N; //Шаг сетки

        x = a;

        for (i = 1; i <= N; i++)

        {

                x = a + (i + 0, 5)*h; //Вычисляем значение в средней точке

                f1 = x * x * atan(x); //подынтегральная функция

                res += f1; //добавляем в сумму

                f2 = h * res;

        }

        cout << "Определенный интеграл по правилу прямоугольников равен " << f2 << endl;

        f1 = x * x * atan(x); //подынтегральная функция

        fx = (x + h) * (x + h) * atan(x + h);

        while (x + h <= b)

        {

                f3 += h * (f1 + fx) / 2.0;

                x += h;

        }

        cout << "Определенный интеграл по правилу трапеции равен " << f3 << endl;