使用Zeta函数计算π的C代码输出无穷大，求修复方案

先看你代码里导致无穷大的核心问题，以及对应的修复方案：

问题分析 #

1. 循环起始值错误：Zeta函数ζ(2)的求和项是从n=1开始的，你代码里循环从i=0启动，pow(0,2)等于0，执行1 / pow(0,2)会触发除零错误，直接让结果变成无穷大。
2. 全局变量滥用：result是全局变量，每次调用zeta都会累积之前的计算值，加上main里两次调用pi(term)，会重复累加，进一步放大错误。
3. 重复调用冗余：main里先调用pi(term)却没保存结果，接着又调用一次，导致计算过程重复执行，结果异常。
4. 浮点除法不明确：1是整数，虽然pow(i,2)返回double会触发隐式转换，但写成1.0能更明确地表达浮点除法意图，避免潜在隐患。

修复后的代码 #

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

double zeta(int terms) {
    double result = 0.0;
    // 从i=1开始循环，符合Zeta函数的数学定义
    for (int i = 1; i <= terms; i++) {
        // 用1.0确保执行浮点除法
        result += 1.0 / pow(i, 2);
    }
    return result;
}

double pi(int term) {
    double stage1 = zeta(term) * 6;
    double stage2 = sqrt(stage1);
    return stage2;
}

int main() {
    printf("the terms?: ");
    int term;
    scanf("%d", &term);
    // 仅调用一次pi函数，保存结果后输出
    double pi_value = pi(term);
    printf("%f\n", pi_value);
    return 0;
}

修复说明 #

• 将result移至zeta函数内部，每次调用都会重新初始化，彻底避免全局变量的累积污染。
• 循环起始值改为i=1，解决除零问题，完全匹配Zeta函数的数学逻辑。
• main函数中仅调用一次pi(term)并保存结果，避免重复计算导致的错误。
• 用1.0替代1，明确指定浮点除法，代码逻辑更清晰。

内容的提问来源于stack exchange，提问作者Yazan Ahmad