利用牛顿方法计算正数的平方根

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int main()
{
  double x,y=1,y1;
  scanf("%lf",&x);
  y1=x/y;
  for(y=1;fabs(y1-y)>=y1*1e-6;)
{  y=y1;
  y1=(y+x/y)/2;
}
printf("%.6lf",y);
    return 0;
}
/*

Hint

是求正数的平方根，因此输入应该定义为double类型。
注意：当y的新旧值之差的绝对值小于0.000001和y的乘积时程序终止，输出y的旧值。*/