Определить корень третьей степени по итерационной формуле

 #include <iostream>
 #include <cmath>
 //--------------------------------------------------------------
 double w( const double &z, const double &ww )
 {
    return  ww + ( z / ww / ww - ww ) / 3;
 }
 //--------------------------------------------------------------
 double cube_root( const double &z, const double &eps )
 {
    double  w_n = 0;
        if ( z >= 0 )
          w_n = z / 3.0;
        else
          w_n = 3*z;
    double  w_n1 = w( z, w_n );
        while ( fabs( w_n - w_n1 ) > eps )
        {
           w_n  = w_n1;
           w_n1 = w( z, w_n );
        }
   return w_n1;
 }
 //--------------------------------------------------------------
 int main()
 {
   double  z   = 1.09;
   double  eps = 0.01;
   std::cout << "approx: " <<  cube_root( z, eps ) << "\n"
             << "exact:  " <<  pow( z, 1.0/3 )     << "\n"
             << "error:  " <<  eps;
   return  0;
 }