用迭代法求x=根号a

公式 x2=0.5*(x1+a/x1);

#include<stdio.h>
#include<math.h>
int main(void){

    double x1,x2;
    double a=0;
    printf("input a number:");
    scanf("%lf",&a);


        x2=a/2;
    do
    {
        x1=x2;
        x2=0.5*(x1+a/x1);
    }while(fabs(x2-x1)>=1e-5);

        printf("result is %lf.",x2);

   }