double newton( double arg ) {       // computes square root
    double r=arg, x=1, eps=0.0001;
    while( fabs(r - x) > eps ) {
       r = x;
       x = r - (r * r - arg) / (2 * r);
       }
    return r;
    }