牛顿迭代法求开根号。 a^1/2_______Xn+1=1/2*(Xn+a/Xn)

#include <stdio.h>
#include <math.h>

int main(void)
{
double a,x1=1.0,x2;
printf("please input a number:\n");
scanf("%lf",&a);
x2=x1;
x1=0.5*(x1+a/x1);
for(;fabs(x1-x2)>=1e-5;)
{
x2=x1;
x1=0.5*(x1+a/x1);
}
printf("a^1/2=%.10lf",x1);
}