Поиск экстремумов, методы золотого сечения, дихотомии, кв. интерполяции-экстраполяции

#include <iostream>
#include <stdio.h>
#include <conio.h>
#include <math.h>
#define eps 0.000000001
#define h 0.005
#define tay 1.618
void menu()
{
std::cout<<"Choose one:\n1) Dichotomy method\n2) Golden method\n3) I-E method\n4) Exit\nYou choose: ";
}
 
double f(double x)
{
    return (0.1*pow(x,3)-2*pow(x,2)+10*x);
}
 
double dichotomy(double infinum, double supremum) 
{
    double x, a=infinum, b=supremum;
    while (fabs(a-b)>eps)
    {
        x=(a + b)/2;
        double f1=f(x-eps);
        double f2=f(x+eps);
        if (f1<f2)
        {
            b=x;
        }
        else if (f1>f2)
        {
            a=x;
        }
          std::cout<<x<<std::endl;
    }
return f(a-b);
}
 
void golden(double infinum, double supremum)
{
    const double fi=1.618;
    double a=infinum, b=supremum, x1, x2, f1, f2, min, max;
    while ((b-a)/2 < eps) 
    {
        x1=b-(b-a)/fi;
        x2=a+(b-a)/fi;
        f1=f(x1);
        f2=f(x2);
        if (f1>=f2)
        {
            a=x1;
            x1=x2;
            x2=a+(b-a)/fi;
        }
        else
            if (f1<=f2)
            {
                b=x2;
                x2=x1;
                x1=b-(b-a)/fi;
            }
    }
min=(a+b)/2;
    while ((b-a)/2 < eps) 
    {
        x1=b-(b-a)/fi;
        x2=a+(b-a)/fi;
        f1=f(x1);
        f2=f(x2);
        if (f1<=f2)
        {
            a=x1;
            x1=x2;
            x2=a+(b-a)/fi;
        }
        else
            if (f1>=f2)
            {
                b=x2;
                x2=x1;
                x1=b-(b-a)/fi;
            }
    }
max=(a+b)/2;
std::cout<<min<<std::endl<<max<<std::endl;
}
 
double ie(double infinum, double supremum)
{
    double a=infinum, b=supremum, c, d;
    bool flag;
    while (flag!=1)
    {
    a=b-h;
    c=b+h;
    d=((f(a)*(pow(b,2)-pow(c,2))+f(b)*(pow(c,2)-pow(a,2))+f(c)*(pow(a,2)-pow(b,2)))/(f(a)*(b-c)+f(b)*(c-a)+f(c)*(a-b)))/2;
    if ((abs(d-b))/abs(b) && (abs(f(d)-f(b))/abs(f(b))))
        flag=1;
    }
return (b+d)/2;
}
 
void main()
{
    double inf=-5, sup=15, extremum;
    int answ;
    while (1==1)
    {
    menu();
    std::cin>>answ;
    switch (answ)
    {
        case 1:
        {
            extremum=dichotomy(inf, sup);
            std::cout<<"Extremum: "<<extremum<<std::endl;
            break;
        }
    case 2:
        {
            golden(inf, sup);
            break;
        }
    case 3:
        {
            extremum=ie(inf, sup);
            std::cout<<"Extremum: "<<extremum<<std::endl;
            break;
        }
    case 4:
        {
        exit(1);
        break;
        }
    default:
        menu();
    }
    }
_getch();
}