Метод деления отрезка пополам C++

[buk_bear]

Вот решение:

#include <iostream>
#include <cmath>
#include <iomanip>

using namespace std;
 
const float eps = 1e-5;
 
float f(const float x)
{
    return x*x*x-8*x*x+11*x+20;
}
 
float Root(const float a, const float b)
{
    if(f(a)*f(b) > 0)
    {
        return log(-1);
    }
    if(fabs(a-b) < eps)
    {
        return (a + b)/2;
    }
    else
    {
        if (f(a)*f((a+b)/2) < 0)
        {
            return Root(a, (a+b)/2);
        }
        else
        {
            return Root((a+b)/2, b);
        }
    }
}
 
int main()
{
	cout<<"Roots"<<endl;
    cout<<setprecision(6)<<Root(-4, 3)<<endl;
    cout<<setprecision(6)<<Root(3, 4.5)<<endl;
    cout<<setprecision(6)<<Root(4.5, 8)<<endl;
}

[digitalis]

   if(f(a)*f(b) > 0)
    {
        return log(-1);
    }

Это чтобы не выдавать сообщение "Корень не найден" воспользовались системным сообщением об ошибке? Однако... Ну и что хорошего? При первой же неудаче прога вылетает, вместо того, чтоб выдать сообщение и продолжать искать других корней.
И даже рекурсия! Могучий источник для сдирания найден. Жаль - препод не впечатлится, потребует объяснить.