贝塞尔函数 c语言程序,贝塞尔函数Bessel（原理与实现）

程序实现：

// crt_bessel1.c

#include

#include

int main( void )

{

double x = 2.387;

int n = 3, c;

printf( "Bessel functions for x = %f:\n", x );

printf( " Kind Order Function Result\n\n" );

printf( " First 0 _j0( x ) %f\n", _j0( x ) );

printf( " First 1 _j1( x ) %f\n", _j1( x ) );

for( c = 2; c < 5; c++ )

printf( " First %d _jn( %d, x ) %f\n", c, c, _jn( c, x ) );

printf( " Second 0 _y0( x ) %f\n", _y0( x ) );

printf( " Second 1 _y1( x ) %f\n", _y1( x ) );

for( c = 2; c < 5; c++ )

printf( " Second %d _yn( %d, x ) %f\n", c, c, _yn( c, x ) );

} 输出结果：

Bessel functions for x = 2.387000: Kind Order Function Result First 0 _j0( x ) 0.009288 First 1 _j1( x ) 0.5229