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

#include <iostream>  // Библиотека ввода ввывода                                                                                                
#include <cmath>  
#include <iomanip>  // библиотека для контроля цифр после запятой 
 
const double e = 2.7;  // число epx
const double pogr = 0.0001; // погрешность 
const double a = 0; // точка а 
const double b = 2; // точка b
 
using namespace std;
 
float vinteg(float way)                // Функция для вычисления интеграла 
{
    float Save = 0.0;                  // Save - сумма всех y
    for (float x = a; x < b; x += way)
        Save += (pow(x,3)*pow(e,3));
    Save *= way;               // умножаем сумму на шаг(way)используя метод прямоугольника.
    return Save;
}
 
float razpervoob()    // функция для вычисления первообразных 
{
    float raz,x;
    raz = fabs(pow(e,3)*(pow(x,3)-3*pow(x,2)+6*x-6));  // вычисляем F(b) - F(a)
    return raz;
}
 
int main()                    //главная функция 
{
    setlocale(0, "");               
    float colway, way;        // colway - количество шагов. way - шаг.
    float otvin2 = razpervoob();
    colway = 2;
    bool yep = false;
 
    do            // Цикл для вычисления кол-во шагов 
    {
        yep = true;
        way = (b - a) / colway;
        if (fabs(fabs(vinteg(way)) - otvin2) >= pogr) { yep = false; colway *= 2.0; }
    } while (yep != true);
 
    float otvin = fabs(vinteg(way));
    cout << "\n\n F(b) - F(a) = " << setprecision(15) << otvin2;      // Вывод F(b) - F(a)
    cout << "\n\n Integral = " << setprecision(15) << otvin;          // Вывод значения интеграла 
    cout << "\n\n Ïîãðåøíîñòü âû÷èñëåíèÿ = " << fabs(otvin - otvin2);        // вывод погрешности 
    cout << "\n\n Êîë-âî ïðîäåëàííûõ øàãîâ = " << setprecision(25) << colway << "\n\n";     // Вывод количества шагов
    system("pause");   
    return 0;
}

raz = fabs((pow(e, 3)*(pow(b, 3) - 3 * pow(b, 2) + 6 * b - 6)) - (pow(e, 3)*(pow(a, 3) - 3 * pow(a, 2) + 6 * a - 6)));

#include <cstdlib> // для system
#include <iostream> 
#include <string>
 
using std::cout;
using std::cin;
using std::endl;
        int n,i;
        double a,b,h,x,s,e;
        double fun(double x);
        double fun_a (double a,double b);
int main()
{    
        setlocale(0, "");
        e=exp(1.0);
        cout<<e<<endl;
        /*
        cout<<"Введите начало интервала a=";
        cin>>a;
        cout<<endl;
        cout<<"Введите конец интервала b=";
        cin>>b;
        cout<<endl;
        cout<<"Введите количество разбиений интервала n=";
        cin>>n;
        cout<<endl;
        */
        n=100;
        a=0;
        b=2;
        h=(b-a)/n;
        x=a;
        s=0;
        for (i = 0; i<n; i++) 
        {
            s=s+0.5*(fun(x)+fun(x+0.5*h))*h;
            x=x+h;
        }
        cout << "Интеграл="<< s <<" "<< fun_a( a,b)<<endl;
        system("pause"); // Только для тех, у кого MS Visual Studio
}
 
double fun (double x)
{
    return pow(x,3)*pow(e,x);
}
double fun_a (double a,double b)
{
    return (pow(e,a)*(pow(a,3)-3*pow(a,2)+6*a-6))-(pow(e,b)*(pow(b,3)-3*pow(b,2)+6*b-6));
}

s=s+fun(x+0.5*h)*h;

fun_a( a,b)

#include <iostream> //ïîäêëþ÷åíèå ñ++
#include <math.h> //áèáëèîòåêà äëÿ âû÷èñëåíèé
#include <cmath>  
using std::cout;
int n, i; 
double a, b, h,s, x, e; //
//ôóíêöèÿ âû÷èñëåíèÿ èíòåãðàëà
float integ(double x)
{
    return pow(x, 3) * pow(e, x);
}
//Ôóíêöèÿ âû÷èñëåíèÿ ïåðâîîáðàçíîé 
float integ_a(double a, double b)
{
    return (pow(e, a) * (pow(a, 3) - 3 * pow(a, 2) + 6 * a - 6)) - (pow(e, b) * (pow(b, 3) - 3 * pow(b, 2) + 6 * b - 6));
}
//ãëàâíàÿ ôóíêöèÿ ïðîãðàììû
int main()
{
    //îáúÿâëåíèå ïåðåìåííûõ ïåðåìåííûõ(êîíñòàíò)
    setlocale(0, ""); 
    e =exp(1.0);
    n = 100;
    a = 0;
    b = 2;
    
    h = (b - a) / n;
    x = a;
    s = 0;
    for (i = 0; i < n; i++)
    {
        s += integ(x + 0.5 * h) * h;
        x  += h;
    }
    cout << "\n\Интеграл: " << s << integ;//âûâîä èíòåãðàëà
     cout <<"\n"  ;
       cout << "\nf(b)-f(a): " << fabs(integ_a(a,b));//âûâîä F(b)-F(a)
       cout <<"\n";
    system("pause");//ïàóçà
}

#include <cstdlib> // для system
#include <iostream> 
#include <string>
 
using std::cout;
using std::cin;
using std::endl;
        int n,i;
        double a,b,h,x,s,e,eps,int_check;
        double fun(double x);
        double fun_a (double a,double b);
int main()
{    
        setlocale(0, "");
        e=exp(1.0);
        /*
        cout<<"Введите начало интервала a=";
        cin>>a;
        cout<<endl;
        cout<<"Введите конец интервала b=";
        cin>>b;
        cout<<endl;
        */
        n=0;
        a=0;
        b=2;
        int_check=fun_a(a,b);
        do{
            n=n+10;
            h=(b-a)/n;
            s=0;
            x=a;
        for (i = 0; i<n; i++) 
        {
            s=s+fun(x+0.5*h)*h;
            x=x+h;
        }
        eps=fabs(s-int_check);
        //cout << "Интеграл="<< s <<" "<< int_check<<" "<<eps<<endl;
        }while (eps>0.001);
        cout << "F(b)-F(a)="<< int_check <<endl;
        cout << "Интеграл="<< s <<endl;
        cout << "Точность="<< eps <<endl;
        cout << "Число шагов разбиения="<< n <<endl;
 
        system("pause"); // Только для тех, у кого MS Visual Studio
}
 
double fun (double x)
{
    return pow(x,3)*pow(e,x);
}
 
double fun_a (double a,double b)
{
    return (pow(e,b)*(pow(b,3)-3*pow(b,2)+6*b-6))-(pow(e,a)*(pow(a,3)-3*pow(a,2)+6*a-6));
}