自己动手实现C标准库中sqrt()函数

今天在C和指针P63看到了这个计算公式，a_i+1 = ( a_i + number/a_i ) / 2。之前也有在网上看到过，不过出处忘了，记得好像是关于组合数学还是数论里面的一个内容。

 

公式很简单，理论这个我就不懂了，还希望各位给指导一下。下面贴出自己的代码，作为记忆保留。

float my_sqrt(float number) {
    float new_guess;
    float last_guess;
    
    if (number < 0) {
        printf("Cannot compute the square root of a negative number!\n");
        return -1;
    }
    
    new_guess = 1;
    do {
        last_guess = new_guess;
        new_guess = (last_guess + number / last_guess) / 2;
        printf("%.15e\n", new_guess);
    } while (new_guess != last_guess);
    
    return new_guess;
}