# include <stdio.h>

double JD = 1.0E-13;

double mysqrt(double x)
{
    if(x < JD)    //x<=0时, 输出
        return -1;

    double p, q, s;

    p = 0;
    q = x;
    s = x / 2;

    while((s*s-x)>JD || (s*s-x)<-JD)
    {
        if(s*s-x > JD)
        {
            q = s;
            s = p + (s-p) / 2;
        }
        else
        {
            p = s;
            s = s + (q-p) / 2;
        }
    }

    return s;
}

int main(void)
{
    int a;

    printf("your num: ");
    scanf("%d", &a);

    printf("sqrt(%d) = %lf\n", a, mysqrt(a));
   
    return 0;
}

表示10的负13次方，很小的一个数。代码里就是让x*x接近s

定义了一个非常小的变量

双精度的全局变量 JD（是不是精度的缩写），赋值10的负13次方

double类型的JD  10的负13次方
代码中用于判断