Решение определенного интеграла методом прямоугольников на C++

Реализуем метод прям угольников на C++. Для этого возьмем функцию

#include "pch.h"
#include "iostream"
#include "cstdlib"
using namespace std;
double metodpryamoygol( double a , double b , double n);
double h = (b — a) / (n); //Шаг интегрирования.
cout << "H=" << h << endl;
for ( double i = 0; i < n; i++)
{
x = a + i * h;
fx = f(x);
summa += fx;
}
return (summa * h);
}
int main() {
setlocale(LC_ALL, «Russian»);
double n = 5; double I = 0; double a, b, eps; //Задаем нижний и верхний предел интегрирования
cout << "Введите a = ";
cin >> a;
cout << "\nВведите b = ";
cin >> b;
cout << "\nВведите точность eps = ";
cin >> eps;
I = metodpryamoygol(a, b, n);
cout << "I= " << I;
cin.get();
}

Теперь введем значения a ,b и eps и посмотрим результат программы