C GSL 贝塞尔函数示例

Question

我尝试了以下代码:

# include <stdio.h> 
# include <gsl/gsl_sf_bessel.h> 
int main (void)
{
double x = 5.0; 
double y = gsl_sf_bessel_j0 (x); 
printf ("JO(%g) =  %.18e\n", x , y);
return(0);
}

我得到的结果是:

JO(5) =  -1.917848549326277019e-01

GSL手册给出的结果是

-1.775967713143382920e-01 

知道哪里出了问题吗？

Answer 1

在我看来这确实像是一个错误。我似乎还没有根据 changelist 修复它.但是，有两个可能的原因:

• 文档中的错字。而不是 gsl_sf_bessel_j0 你应该使用 gsl_sf_bessel_J0 (非常奇怪的命名)。这给你 JO(5) = -1.775967713143382642e-01 这足够接近所需的答案(虽然我不确定精度)。
• 文档是正确的，我们需要修复函数本身。

无论如何，由于没有人报告我已经提交了错误报告。我确实计划在获得与此问题相关的任何信息后立即更新此答案。