自己動手實現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;
}