funkcja_pole

#include <iostream>
#include <cmath>
using namespace std;

double funkcja(double x)
{
    //return (x*x)/3+2;
    return pow(x,3)-(4*x);
}

int main()
{
    double n,a=3,b=4,bp=a,fa=funkcja(a),fb=funkcja(a),wynik=0;
    cin>>n;
    double dlugosc=(b-a)/n;
    for(int i=0;i<n;i++)
    {
        bp+=dlugosc;
        fa=fb;
        fb=funkcja(bp);
        wynik+=(dlugosc*(fa+fb))/2;
    }
    cout <<wynik;
    return 0;
}