Программа по поиску минимально и максимального значения функции

#include <iostream>
#include <cmath>
double f(double x);
 
using namespace std;
int main()
{
 
    double x0, h, x, fx;
    int N;
    cout << "Enter x begin > ";     cin >> N;
    cout << "Enter x end > ";     cin >> x0;
    cout << "Enter Dx > ";        cin >> h;
    double x = N + 1;
    double fx = N + 1;
    int min = 0, max = 0;
    for (int i = 0; i <= N; i++)
    {
        x = x0 + i*h;
        fx = f(x);
    }
    for (int i = 1; i <= N; i++)
    {
        if (fx > max)
            max = i;
        if (fx < min)
            min = i;
    }
    cout << "max = " << max << endl;
    cout << "min = " << min << endl;
    cout << "raznost = " << max - min << endl;
    system("pause");
    return 0;
}
double f(double x)
{
    return 4 * x *x*x - 3 * x *x + 2 * x - 1;
}

#include <iostream>
#include <cmath>
double f(double x);
 
using namespace std;
int main()
{
    double xbeg, xend, h, x, fx;
    cout << "Enter x begin > ";     cin >> xbeg;
    cout << "Enter x end > ";     cin >> xend;
    cout << "Enter Dx > ";        cin >> h;
    double fmin,  fmax;
    fmax = fmin = f(xbeg);
    for (x = xbeg+h; x <= xend; x+=h)
    {
        fx = f(x);
        if (fx > fmax) fmax = fx;
        if (fx < fmin) fmin = fx;
    }
    cout << "max = " << fmax << endl;
    cout << "min = " << fmin << endl;
    cout << "raznost = " << fmax - fmin << endl;
    system("pause");
    return 0;
}
double f(double x)
{
    return 4 * x *x*x - 3 * x *x + 2 * x - 1;
}