使用泰勒级数计算cosh(x)的CPP代码问题：大x值计算偏差

首先，咱们来拆解你代码里的几个核心问题，然后一步步优化解决：

1. 阶乘计算的溢出隐患 #

你的fact函数里用int类型存储阶乘结果，哪怕最后返回long double，中间计算过程早就溢出了。比如13的阶乘是6227020800，已经超过了32位int的最大值（2147483647），溢出后会变成负数或者乱码，直接导致后续除法结果完全错误。

2. 低效且不稳定的泰勒项计算方式 #

直接调用pow(a,c)和fact(c)来计算每一项，不仅会重复计算大量冗余值，还会因为大指数和大阶乘的数值极端化，导致精度损失甚至溢出。其实泰勒级数的每一项都可以通过前一项递推得到：

cosh(x) = 1 + x²/2! + x⁴/4! + x⁶/6! + ...
第n项（从0开始，n为偶数）：termₙ = termₙ₋₂ * x² / (n*(n-1))

这样完全不需要单独计算阶乘和幂，数值稳定性会好很多。

3. 多余的绝对值操作 #

cosh(x)是偶函数，所有项都是非负的，sum不可能为负，这个判断完全没必要，反而可能在数值异常时引入错误。

修正后的代码 #

下面是优化后的实现，解决了溢出问题，同时提升了计算精度和效率：

#include <iostream>
#include <cmath> // 用std::fabs处理x的绝对值，利用cosh的偶函数特性
using namespace std;

int main() {
    long double x;
    cout << "Enter value of x: " << endl;
    cin >> x;

    // cosh是偶函数，先取绝对值减少计算量
    x = fabs(x);

    long double sum = 1.0; // 泰勒级数的第一项是1
    long double term = 1.0; // 初始项对应x⁰/0! = 1
    int n = 2; // 从第二项（x²/2!）开始递推

    // 循环直到项足够小，保证精度，而非固定循环次数
    while (fabs(term) > 1e-15) {
        term *= (x * x) / (n * (n - 1));
        sum += term;
        n += 2;
    }

    cout << "The answer is: " << sum << endl;
    system("pause");
    return 0;
}

关键优化点说明 #

• 递推计算每一项：彻底避免了大阶乘和大幂的直接计算，从根源上解决溢出问题，同时减少了计算量。
• 动态循环终止条件：不再固定循环12次，而是当项的绝对值小于1e-15时停止，既保证精度，又不会在x较小时做多余计算。
• 利用偶函数特性：先对x取绝对值，因为cosh(-x)=cosh(x)，简化计算逻辑。
• 类型一致性：所有变量都用long double，避免中间类型转换带来的精度损失。

这样修改后，即使x取较大的值（比如x=20），也能得到准确的计算结果，不会再出现阶乘溢出导致的错误。

内容的提问来源于stack exchange，提问作者Mohammad Hassaan