C语言牛顿迭代法求平方根

double mysqrt(double a)
{
double x,y;
x=1.0;
    
while(x*x-a<-0.00001||x*x-a>0.00001)
{
 y=(x+a/x)/2.0;
 x=y;
}
return x;

}

0的计算不太准确，其他正常