迭代法求算术平方根，是牛顿?算出来不对

#include<stdio.h>
int main()
{
double x0,x1,a;
scanf("%lf",&a);
x0=a/2.0;
x1=(x0+a/x0)/2.0;
while((x0-x1)>=1e-5)
{ x0=x1;
x1 = (x0+a/x0)/2.0;
}
printf("%.4lf",x1);
return 0;
}