求开方的函数

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);

}