Метод деления отрезка пополам для решения нелинейных уравнений (метод дихотомии)

#include <conio.h>
#include <math.h>
#include <iostream>
#include <fstream>
using namespace std;
double func(double x)
    {
        //return (x * x * x - x - 1);//уравнение
        return (x * x * x + 2 * x * x + 3 * x + 5);
    }
double dichotomy(double infinum, double supremum, double epsilon) {
        double x;
        while (supremum - infinum > epsilon) {
            x = (infinum + supremum) / 2;
            if (func(supremum) * func(x) < 0)
                infinum = x;
            else
                supremum = x;
        }
        return (infinum + supremum) / 2;
    }
int main () {
        freopen ("out.txt", "w", stdout);
        printf ("Хприближенное = %0.3f", dichotomy(-2,-1,0.001));
    }

//------------------------------------------------------------------------------
//исключения
class bad_bracket_end_points {};
class too_many_iteration {};
//------------------------------------------------------------------------------
long root_iteration = 1000000;
//------------------------------------------------------------------------------
double dichotomic_root
(
    double (* const y)(const double &x),
    double x_min, double x_max,
    const double &dx, const double &dy
)
{
    double y_x_min = y(x_min);
    double y_x_max = y(x_max);
 
    double x_mid, y_x_mid;
 
    if (x_min >= x_max || (y_x_min > 0) == (y_x_max > 0))
        throw bad_bracket_end_points();
 
    long k = 0;
 
    do
    {
        x_mid = (x_min + x_max) / 2;
        y_x_mid = y(x_mid);
 
        if ((y_x_min > 0) == (y_x_mid > 0))
            x_min = x_mid;
        else
            x_max = x_mid;
 
        k++;
        if (k > root_iteration) throw too_many_iteration();
    }
    while (x_max - x_min > dx && fabs(y_x_mid) > dy);
 
    return x_mid;
}