第35讲项目3——Bessel函数

任务及要求

/*
*Copyright (c)2016,csdn学院
*All rights reserved.
*文件名称：lzz.c
*作    者：兰泽泽
*完成日期：2016年3月16日
*版 本 号：v1.0
*问题描述：Bessel函数Jn(X)
*程序输入：
*程序输出：
*/
#include <stdio.h>
#include <math.h>
int main()
{
    double jn, j0, j1, x;
    int n, i;
    printf("please input n&x :");
    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;
}

运行结果