【项目3:Bessel函数】
Bessel函数Jn(X)有以下的递推关系:
编写程序,利用递推关系,由任意的n和x≠0求Jn(X)。
Bessel函数Jn(X)有以下的递推关系:
编写程序,利用递推关系,由任意的n和x≠0求Jn(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代表的实际是公式中的jn+1,推导一下,2*i-1,而不是2*i+1,这是对的