python 贝塞尔函数_第二类整数阶贝塞尔函数(诺伊曼函数) | 学步园

#include "stdio.h"

#include "math.h"

/************************************************************

第一类贝塞尔函数:贝塞尔函数

*************************************************************/

double Bessel(int n,double x)

{

int i,m;

double t,y,z,p,q,s,b0,b1;

static double a[6]={ 57568490574.0,-13362590354.0,

651619640.7,-11214424.18,77392.33017,

-184.9052456};

static double b[6]={ 57568490411.0,1029532985.0,

9494680.718,59272.64853,267.8532712,1.0};

static double c[6]={ 72362614232.0,-7895059235.0,

242396853.1,-2972611.439,15704.4826,

-30.16036606};

static double d[6]={ 144725228443.0,2300535178.0,

18583304.74,99447.43394,376.9991397,1.0};

static double e[5]={ 1.0,-0.1098628627e-02,

0.2734510407e-04,-0.2073370639e-05,

0.2093887211e-06};

static double f[5]={ -0.1562499995e-01,

0.1430488765e-03,-0.6911147651e-05,

0.7621095161e-06,-0.934935152e-07};

static double g[5]={ 1.0,0.183105e-02,

-0.3516396496e-04,0.2457520174e-05,

-0.240337019e-06};

static double h[5]={ 0.4687499995e-01,

-0.2002690873e-03,0.8449199096e-05,

-0.88228987e-06,0.105787412e-06};

t=fabs(x);

if (n<0) n=-n;

if (n!=1)