Разработать алгоритм и программу приближенного вычисления интеграла методом прямоугольников

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
double f(double x)
{
    return (2*pow(x, 2) - 7*x + 4);
}
int main()
{
    double a, b, S, h, x;
    int n, i;
    
    cout << "Vvedite a: "; cin >> a; cout << "\n";
    cout << "Vvedite b: "; cin >> b; cout << "\n";
    cout << "Vvedite h: "; cin >> h; cout << "\n";
    /*
    a = 4;
    b = 5;
    h = 0.001;
    */
 
    S = (f(a)+f(b))/2;
    //h = (b-a)/n;
    n = (b-a)/h;
 
    for (int i=1; i<=n-1; i++)
    {
        S += f(a+i*h);
    }
    double I = h*S;
    cout << setprecision(10) << "Int = " << I << endl;
    
    return 0;
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
double f(double x)
{
    return (2*pow(x, 2) - 7*x + 4);
}
int main()
{
    double a, b, S, h, x;
    int n, i;
    
    cout << "Vvedite a: "; cin >> a; cout << "\n";
    cout << "Vvedite b: "; cin >> b; cout << "\n";
    cout << "Vvedite h: "; cin >> h; cout << "\n";
    /*
    a = 4;
    b = 5;
    h = 0.001;
    */
        
    n = (b-a)/h;    
    S = 0;
    for(i = 1; i <= n; i++)
    {
        x = a + i*h;
        S += f(x);
    }
    I = h*S;
    cout << "\nI = " << I;
 
    return 0;
}