Рекурсивная функция, которая методом деления отрезка пополам находит корень уравнения

 #include <iostream>
 #include <cmath>
 
 double  f( const double &x )
 {
   return x*x - 81;
 }
 //--------------------------------------------------------------
 // Половинное деление.
 // eps > 0, а < b, и f(x) непрерывная на отрезке [а, b] функция,
 // имеющая на этом отрезке единственный корень.
 //--------------------------------------------------------------
 double dihotomia( const double   &a,
                   const double   &b,
                   const double   &eps )
 {
    double  ksi = ( a + b ) / 2;
       if (    ( fabs( f( a ) - f( b ) ) <= eps )
            || ( fabs( f( ksi) )         <= eps )  )
          return ( a + b ) / 2;
 
       if ( f( a )*f( ksi) <= 0.0 )
          dihotomia( a, ksi, eps );
       else
          dihotomia( ksi, b, eps );
 }
 //--------------------------------------------------------------
 int  main()
 {
   const double  a   = -1;
   const double  b   =  10;
   const double  eps = 0.001;
   std::cout << dihotomia( a, b, eps );
   //std::cin.get();
   return 0;
 }