【项目3:Bessel函数】
Bessel函数Jn(X)有以下的递推关系:
编写程序,利用递推关系,由任意的n和x≠0求J n(X)。
Bessel函数Jn(X)有以下的递推关系:
编写程序,利用递推关系,由任意的n和x≠0求J n(X)。
[参考解答]
#include <stdio.h>
#include <math.h>
int main( )
{
double jn, j0, j1, x;
int n, i;
scanf("%d %lf", &n, &x);
j0=sin(x)/x;
j1=sin(x)/(x*x) - cos(x)/x;
for (i=2; i<=n; i++)
{
jn=(2*i-1) * j1 / x - j0; //见后注
j0 = j1;
j1 = jn;
}
printf("Jn = %.4f\n", jn);
return 0;
}
注:变量jn代表的实际是公式中的j
n+1,推导一下,2*i-1,而不是2*i+1,这是对的