如何在纯C代码中调用C++贝塞尔函数库及处理头文件

方案1：直接使用C标准库的贝塞尔函数 #

C99及之后的标准库在<math.h>中提供了cyl_bessel_i函数，和C的std::cyl_bessel_i功能完全一致，直接调用即可。输出用printf替代C的std::cout：

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

int main()
{
    // 验证nu=0时的贝塞尔函数值
    const double x = 1.2345;
    printf("I_0(%lf) = %lf\n", x, cyl_bessel_i(0, x));

    // I_0的级数展开计算
    double fct = 1;
    double sum = 0;
    for (int k = 0; k < 5; fct *= ++k)
    {
        double term = pow(x / 2, 2 * k) / pow(fct, 2);
        sum += term;
        printf("sum = %lf\n", sum);
    }

    return 0;
}

编译时需要链接数学库（数学库独立于标准C库），以GCC为例：

gcc your_file.c -o your_program -lm

方案2：手动实现级数展开（无标准库贝塞尔函数依赖） #

如果你的环境不支持C99的cyl_bessel_i，可以直接手动实现原C++代码中的级数展开逻辑：

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

int main()
{
    const double x = 1.2345;
    // 手动计算I_0的级数展开
    double fct = 1;
    double sum = 0;
    for (int k = 0; k < 5; fct *= ++k)
    {
        double half_x = x / 2;
        double numerator = pow(half_x, 2 * k);
        double denominator = pow(fct, 2);
        sum += numerator / denominator;
        printf("sum = %lf\n", sum);
    }
    printf("I_0(%lf) = %lf\n", x, sum);

    return 0;
}

关键说明 #

• C的<math.h>包含了绝大多数C++<cmath>中的数学函数，无需特殊“引入”操作，直接包含头文件后使用无std::前缀的函数名即可。
• 无论哪种方案，编译时都必须加上-lm参数完成数学库的链接。

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