求开方的函数

double   _sqrt(double   a)
{
    double   x1   =   0.0;
    double   x2   =   a/2;
   
    while(x1!=   x2)
    {
          x1   =   x2;
          x2   =   (x1+a/x1)/2;
    }
 cout<<x1<<endl;
    return   x1;
 
}

void main()
{
 _sqrt(9);

}