00001 00029 #ifndef BESSEL_INTERNAL_H 00030 #define BESSEL_INTERNAL_H 00031 00033 00034 #include <cmath> 00035 00036 double chbevl(double x, double array[], int n); 00037 double hyperg(double a, double b, double x); 00038 int airy(double x, double *ai, double *aip, double *bi, double *bip); 00039 double polevl(double x, double coef[], int N); 00040 double p1evl(double x, double coef[], int N); 00041 00042 double i0(double x); 00043 double i0e(double x); 00044 double i1(double x); 00045 double i1e(double x); 00046 00047 double k0(double x); 00048 double k0e(double x); 00049 double k1(double x); 00050 double k1e(double x); 00051 00052 double iv(double nu, double x); 00053 double jv(double nu, double x); 00054 double yv(double nu, double x); 00055 double kn(int n, double x); 00056 00057 double gam(double x); 00058 double lgam(double x); 00059 extern int sgngam; 00060 00062 00063 #endif // #ifndef BESSEL_INTERNAL_H