Найти значение интеграла используя метод средних прямоугольников

#include <iostream>
#include <stdio.h>
#include <math.h>
 
using namespace std;
 
void IRect(double A, double B, double E, double &I)
//численное интегрирование функций методом средних прямоугольников
//A,B - интервал интегрирования
//E - заданная точность
//Y -  значение интеграла
{
    double Res, //значение интеграла
           ts, //погрешность
           h, //шаг
           z, //значение интеграла на предыдущем шаге
           s, //интегральная сумма
           x; //значение аргумента
    int j, //параметр цикла
        N; //число отрезков разбиения
    N = 100; //начальное число отрезков разбиения
    Res = 1E10; //начальное значение интеграла
    do //цикл вычисления интеграла с заданной точностью
    {
        N = N * 2; //удваиваем число отрезков разбиения
        z = Res; //запомнили предыдущее значение интеграла
        h = (B - A) / N; //шаг
        s = 0; //значение интегральной суммы
        for (j = 1; j <= N; j++) //вычисление интеграла по формуле средних
        {                        //прямоугольников для заданного числа отрезков
            
            x = A + j * h; //вычислили значение аргумента
            s = s + (5 * x * sin(x)); //добавили в интегральную сумму
        }
        Res = s * h; //текущее значение интеграла
        ts = fabs(Res - z)/24; //вычислили погрешность
    }
    while (ts >= E); //пока погрешность не станет меньше точности
    I = Res - ts;
}
 
int main(int argc, char** argv)
{
    double a, b;//пределы интегрирования
    double eps;//точность интегрирования
    double I;//значение интеграла
    setlocale(LC_ALL,"Russian");
    cout << "Введите интервал пределы интегрирования a, b: ";
    cin >> a >> b;//ввод
    cout << "Введите точность eps: ";
    cin >> eps;//ввод
    IRect(a, b, eps, I);
    cout << "Значение интеграла: " << I << endl;//вывод интеграла
    system("pause");//задержка экрана до нажатия любой клавиши
    return 0;
}