double sqrt(double n) {
    double x,newx;
    x = n;
    for (unsigned i = 0; i < 10; i++) {
        newx = 0.5*(x + n/x);
        printf("%u\t%g\t%g\n",i,x,newx);
        x = newx;
    }
    return x;
}