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

main() {
  double C, xold = 1.0, xnew;
  double epsilon = 0.0001;
  int converging;
  printf("dai un numero positivo ");
  scanf("%lf",&C);
  converging = 0;
  while (!converging) {
    xnew = 0.5*(xold + C/xold);
    if(fabs(xnew-xold) < epsilon) {converging = 1;}
    xold = xnew;
  }
  printf("la radice di %lf e' %lf\n",C,xnew);
}