Методы вычисления корней не линейных уравнений

//f(x)=x^6-3x^2+x-1

#include <iostream>
#include <iomanip>
#include <vector>
#include <math.h>
#include <fstream>
 
using namespace std;
 
//f(x)=x^6-3x^2+x-1
 
struct MyStruct
{   double x, min, max;//min - нижний предел, max - верхний предел
    double eps; // точность
};
 
double f(double x)//функция 
{
    //  return pow(2.0,x) - 2*pow(x, 2.0) - 1;
    return pow(x, 6.0) - 3 * pow(x, 2.0) + x - 1;
}
 
double df(double x)//производная
{
    return pow(x, 5) * 6 - x * 6 + 1;
}
 
double g_iter(double x)
{
    return x + 0.05*f(x);
}
 
double g_newton(double x)
{
    return x - f(x) / df(x);
}
 
double MtIter(double x, double eps, double min, double max)
{   x = abs(max-min)/2;
 
    //cout << "Enter initial root value   : "; cin >> x;
    for (double iter = 1; eps < fabs(f(x)); iter = iter + 1)
    {
        system("cls");
        cout << "Iteration : " << setprecision(0) << iter << endl;
        cout << "x    = " << x << endl;
        //cout << "f(g(x)) = " << f(g(x)) << endl;
        //cout << "g(x) = " << g(x) << endl;
        //cout << "f(x) = " << f(x) << endl;
 
    if (f(x)*f(g_iter(x)) > 0)
            x = g_iter(x);
        else
            break;
    }
    system("pause");
    return x;
}
double MtDih(double x, double eps, double min, double max)
{   double c; 
    while (abs(max - min) > eps) { 
        system("cls");
        c = (min + max) / 2; 
        if ((f(c) > 0 && f(min) > 0) || (f(c) < 0 && f(min) < 0)) {
            if (min == c)
            {break;}
                min = c;
        } else 
        { if(max == c)
        {break;}
            max = c; 
        } 
    } 
    cout << "корень: " << (min + max) / 2; 
    x = (min + max) / 2;
    return x, min, max; 
} 
 
double MtNewton(double x, double eps, double min,double max)
{
    x = abs(max - min) / 2;
    for (double iter = 1; eps < fabs(f(x)); iter = iter + 1)
    {
        system("cls");
        //*Итераций может быть очень много, поэтому рекомендую забыть
        //о целых а использовать дабл как счётчик, хотя в принципе если 
        //решение не нашли за 10-100 итераций то решения для данного коэффициента
        //при f(x) в g(x) нет и надо его менять
        cout << "Iteration : " << setprecision(0) << iter << endl;
        if (df(x) == 0)//Чёртовский важный момент(!)
        {       break;
    }//ведь если df(x) == 0, то будет деление на ноль x - f(x)/df(x)
        x = g_newton(x);
        
        cout << "x    = " << x << endl;
        //cout << "g(x) = " << g(x) << endl;
        //cout << "df(x)= " << df(x) << endl;
        //cout << "f(x) = " << f(x) << endl;
        
    }
    system("pause");
    return x;
}
 
 
void main()
{
    setlocale(LC_ALL, "");//подключение языков
    vector<MyStruct> date;
 
    MyStruct boofer;
    
    int menu = -1;
    while (menu = -1)
    {
    boofer.x = 0;
    cout << "введите начало и конец отрезка:";
    boofer.min;
    boofer.max;
    cin >> boofer.min >> boofer.max;
    cout << "Введите эпсилон"; cin >> boofer.eps;
    boofer.eps;
 
        cout << "(1) - метод итераций" << endl << "(2) - метод Ньютона" << endl << "(3) - метод деления на 2" << endl << "(0) - выход" << endl;
        cin >> menu;
        switch (menu)
        {
        case 1:
            MtIter(boofer.x, boofer.min, boofer.max, boofer.eps);
            break;
        case 2:
            MtNewton(boofer.x, boofer.min, boofer.max, boofer.eps);
            break;
        case 3:
            MtDih(boofer.x, boofer.min, boofer.max, boofer.eps);
            break;
        default:
            cout << "Такого ответа нет";
            break;
        }
    }
}