c函数返回nan而不是double

Question

我们在我的编程课中学习数值方法，向我们介绍的第一个算法是求根的二分法。这是我使用递归实现它的尝试：

#include <math.h>
#include <stdio.h>

#define tolerance 0.00001

double my_function(double z){
    double answer = 5*pow(z,2) + 5*z - 2;

    return answer;
} 

double bisection(double (*fxn)(double),double a, double b){

    double m = ((a+b)/2);

    if (fabs(b-a) < tolerance){
        double root = a;
        printf("value of a is %lf\n",root);
        return a;
    }



    else if (fxn(m) > 0){
        b = m;
    }

    else if (fxn(m) < 0){
        a = m;  
    } 

    bisection(my_function, a, b);
}

int main(){

    double g = 0.01;
    double z = 1;

    double x = bisection(my_function,g,z);
    printf("root is %lf\n",x);  

return 0;
}

这是输出：

value of a is 0.306225
root is nan

根是正确的（稍微偏离，但在容差范围内），但在返回值和打印它之间的某个地方，它以某种方式变成了 NaN。我难住了。我做错了什么？

Answer 1

我的第一个猜测：

在部分

if (fabs(b-a) < tolerance){
    double root = a;
    printf("value of a is %lf\n",root);
    return a;
}

您返回 a 而不是 root。尝试返回 root 看看是否有帮助。

Answer 2

您没有从递归调用中返回。将bisection中的最后一条语句改为

return bisection(my_function, a, b);