Численное интегрирование

#include <iostream>
#include <cmath>
# include <locale.h>
using namespace std;
 
double func(const double &x)
{
    return (2 * x*x*x - 7 * x + 4);
}
 
double integr_lefttriangle(const int &n, const double &a, const double &b)
{
    double sum = 0;
    double h = (b - a) / n;
    for (double i = 0; i < n; i++)
        sum += func(a + i * h);
    return sum * h;
}
double integr_righttriangle(const int &n, const double &a, const double &b)
{
    double sum = 0;
    double h = (b - a) / n;
    double x = a + h; 
    
    for (int i = 0; i < n; i++, x += h) 
    {
        sum += func(x);
    }
    return h * sum;
}
double integr_middle(const int &n, const double &a, const double &b)
{
    double sum = (func(a) + func(b)) / 2;
    double h = (b - a) / n;
    for (int i = 1; i <= n - 1; i++)
    {
        sum += func(a + i * h);
    }
    return h * sum;
}
 
double I(double a, double b, int n, double y) { return ((b - a) / (2 * n)*y); }
double trapeze(const int &n, const double &a, const double &b)
{
    
        double h, sum, x;
        sum = 0;
        h = (b - a) / n;
        sum += func(a);
        for (int i = 1; i < n; i++)
        {
            x = a + i * h;
            sum += 2 * func(x);
        }
        x = a + n * h;
        sum += func(x);
        return(0.5*h*sum);
}
double simpson(const int &n, const double &a, const double &b)
{
    double h = (b - a) / n;
    double sum = 0;
 
    double x0 = a;
    double x1 = a + h;
 
    for (int i = 0; i <= n - 1; i++) {
        sum += func(x0) + 4 * func(x0 + h / 2) + func(x1);
 
        x0 += h;
        x1 += h;
    }
 
    return (h / 6)*sum;
 
}
 
int main()
{
    setlocale(LC_ALL, "rus");
    double a, b;
    int n = 50;
    cout << "Введите нижний предел интегрирования" << endl;
    cin >> a;
    cout << "Введите верхний предел интегрирования" << endl;
    cin >> b;
 
    cout << "Левый прям: " << integr_lefttriangle(n, a, b) << endl;
    cout << "Правый прям: " << integr_righttriangle(n,a,b) << endl;
    cout << "Средние прям: " << integr_middle(n,a,b) << endl;
    cout << "Трапеция: " << trapeze(n, a, b) << endl;
    cout << "Симпсон: " << simpson(n, a, b) << endl;
    system("pause");
    return 0;
}