Что не так с моей логикой для нахождения квадратного корня по методу Ньютона?

Question

Я написал следующий код для нахождения квадратного корня методом Ньютона с помощью последовательных приближений, но он не дает мне правильного ответа. Может кто-нибудь объяснить это?

#include<stdio.h>
#include<stdlib.h>
#define square(x) x*x
double rootByNewtonApprox(int n);
double improve(double n);
double average(double a,double b);
int goodEnough(double guess);
double guess(int n);
int number;
int main(void)
{

    double root;
    printf("\nEnter the number you want square root of: ");
    scanf("%d",&number);
    if(number<0)
            number = -1* number;
    root = rootByNewtonApprox(number);
    printf("\nThe square root of %d is %lf\n",number,root);
    return 0;
}
double guess(int n)
{
    return n/2;
}
double rootByNewtonApprox(int n)
{
    if(goodEnough(guess(n)))
            return guess(n);
    else
            rootByNewtonApprox(improve(guess(n)));
}

double improve(double guess)
{
    return average(guess,(number/guess));
}
double average(double a,double b)
{
    return ((a+b)/2);
}
int goodEnough(double guess)
{
    if(abs(square(guess) - number) <= 0.001)
            return 1;
    else
            return 0;
}

Теперь, когда я даю n = 2, он дает выходной сигнал nan, а когда я даю n = 9, он говорит segmentation Fault.

Answer 1

double guess(int n)
{
    return  n / (double) 2;
}

Answer 2

Вы забыли один return

double rootByNewtonApprox(int n)
{
if(goodEnough(guess(n)))
    return guess(n);
else
    return rootByNewtonApprox(improve(guess(n)));
    ^
}