T015 迭代求根

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

int main()
{
    double a, x0, x1;
    scanf("%lf", &a);
    if(a<0) printf("Error\n");
    else
    {
        x0=a/2;//设初值
        x1=(x0+a/x0)/2;
        do
        {
            x0=x1;
            x1=(x0+a/x0)/2;
        }   while(fabs(x0-x1)>1e-3);
    }
    printf("%.5lf\n", x1);
    return 0;
}