牛頓疊代法求開根号。 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);

}