#define ABS(val) (((val)>0)?(val):(-(val)))

double my_sqrt(double a)
{
    double k = 1.0;
    while(ABS(k*k-a)>1e-9) {
        k = (k+a/k)/2;
    }
    return k;
}