牛顿迭代法求开根号。 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:
");
	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);
}