23 #if (AE_COMPILER==AE_MSVC) 24 #pragma warning(disable:4100) 25 #pragma warning(disable:4127) 26 #pragma warning(disable:4702) 27 #pragma warning(disable:4996) 68 return *(
reinterpret_cast<double*
>(&result));
116 return *(
reinterpret_cast<double*
>(&result));
158 return *(
reinterpret_cast<double*
>(&result));
201 return *(
reinterpret_cast<double*
>(&result));
247 return *(
reinterpret_cast<double*
>(&result));
269 return *(
reinterpret_cast<double*
>(&result));
310 return *(
reinterpret_cast<double*
>(&result));
355 return *(
reinterpret_cast<double*
>(&result));
404 return *(
reinterpret_cast<double*
>(&result));
456 return *(
reinterpret_cast<double*
>(&result));
496 void airy(
const double x,
double &ai,
double &aip,
double &bi,
double &bip)
550 return *(
reinterpret_cast<double*
>(&result));
586 return *(
reinterpret_cast<double*
>(&result));
629 return *(
reinterpret_cast<double*
>(&result));
673 return *(
reinterpret_cast<double*
>(&result));
710 return *(
reinterpret_cast<double*
>(&result));
748 return *(
reinterpret_cast<double*
>(&result));
785 return *(
reinterpret_cast<double*
>(&result));
822 return *(
reinterpret_cast<double*
>(&result));
859 return *(
reinterpret_cast<double*
>(&result));
894 return *(
reinterpret_cast<double*
>(&result));
932 return *(
reinterpret_cast<double*
>(&result));
962 double beta(
const double a,
const double b)
970 return *(
reinterpret_cast<double*
>(&result));
1027 return *(
reinterpret_cast<double*
>(&result));
1071 return *(
reinterpret_cast<double*
>(&result));
1118 return *(
reinterpret_cast<double*
>(&result));
1167 return *(
reinterpret_cast<double*
>(&result));
1211 return *(
reinterpret_cast<double*
>(&result));
1239 return *(
reinterpret_cast<double*
>(&result));
1272 return *(
reinterpret_cast<double*
>(&result));
1375 return *(
reinterpret_cast<double*
>(&result));
1422 return *(
reinterpret_cast<double*
>(&result));
1458 return *(
reinterpret_cast<double*
>(&result));
1499 return *(
reinterpret_cast<double*
>(&result));
1545 return *(
reinterpret_cast<double*
>(&result));
1596 return *(
reinterpret_cast<double*
>(&result));
1646 return *(
reinterpret_cast<double*
>(&result));
1689 return *(
reinterpret_cast<double*
>(&result));
1735 return *(
reinterpret_cast<double*
>(&result));
1776 return *(
reinterpret_cast<double*
>(&result));
1823 return *(
reinterpret_cast<double*
>(&result));
1871 return *(
reinterpret_cast<double*
>(&result));
1923 return *(
reinterpret_cast<double*
>(&result));
1974 return *(
reinterpret_cast<double*
>(&result));
2058 return *(
reinterpret_cast<double*
>(&result));
2087 return *(
reinterpret_cast<double*
>(&result));
2194 return *(
reinterpret_cast<double*
>(&result));
2222 return *(
reinterpret_cast<double*
>(&result));
2273 return *(
reinterpret_cast<double*
>(&result));
2302 return *(
reinterpret_cast<double*
>(&result));
2368 return *(
reinterpret_cast<double*
>(&result));
2410 return *(
reinterpret_cast<double*
>(&result));
2445 return *(
reinterpret_cast<double*
>(&result));
2500 return *(
reinterpret_cast<double*
>(&result));
2557 return *(
reinterpret_cast<double*
>(&result));
2590 return *(
reinterpret_cast<double*
>(&result));
2717 static double gammafunc_gammastirf(
double x,
ae_state *_state);
2726 static void bessel_besselmfirstcheb(
double c,
2731 static void bessel_besselmnextcheb(
double x,
2737 static void bessel_besselm1firstcheb(
double c,
2742 static void bessel_besselm1nextcheb(
double x,
2748 static void bessel_besselasympt0(
double x,
2752 static void bessel_besselasympt1(
double x,
2760 static double ibetaf_incompletebetafe(
double a,
2766 static double ibetaf_incompletebetafe2(
double a,
2772 static double ibetaf_incompletebetaps(
double a,
2809 static void trigintegrals_chebiterationshichi(
double x,
2842 #ifndef ALGLIB_INTERCEPTS_SPECFUNCS 2873 z =
ae_pi/(z*gammafunc_gammastirf(q, _state));
2877 z = gammafunc_gammastirf(x, _state);
2892 result = z/((1+0.5772156649015329*
x)*x);
2902 result = z/((1+0.5772156649015329*
x)*x);
2914 pp = 1.60119522476751861407E-4;
2915 pp = 1.19135147006586384913E-3+x*
pp;
2916 pp = 1.04213797561761569935E-2+x*
pp;
2917 pp = 4.76367800457137231464E-2+x*
pp;
2918 pp = 2.07448227648435975150E-1+x*
pp;
2919 pp = 4.94214826801497100753E-1+x*
pp;
2920 pp = 9.99999999999999996796E-1+x*
pp;
2921 qq = -2.31581873324120129819E-5;
2922 qq = 5.39605580493303397842E-4+x*qq;
2923 qq = -4.45641913851797240494E-3+x*qq;
2924 qq = 1.18139785222060435552E-2+x*qq;
2925 qq = 3.58236398605498653373E-2+x*qq;
2926 qq = -2.34591795718243348568E-1+x*qq;
2927 qq = 7.14304917030273074085E-2+x*qq;
2928 qq = 1.00000000000000000320+x*qq;
2932 return _ialglib_i_gammafunction(x);
2971 #ifndef ALGLIB_INTERCEPTS_SPECFUNCS 2989 logpi = 1.14472988584940017414;
2990 ls2pi = 0.91893853320467274178;
3012 result = logpi-
ae_log(z, _state)-
w;
3043 result =
ae_log(z, _state);
3048 b = -1378.25152569120859100;
3049 b = -38801.6315134637840924+x*
b;
3050 b = -331612.992738871184744+x*
b;
3051 b = -1162370.97492762307383+x*
b;
3052 b = -1721737.00820839662146+x*
b;
3053 b = -853555.664245765465627+x*
b;
3055 c = -351.815701436523470549+x*
c;
3056 c = -17064.2106651881159223+x*
c;
3057 c = -220528.590553854454839+x*
c;
3058 c = -1139334.44367982507207+x*
c;
3059 c = -2532523.07177582951285+x*
c;
3060 c = -2018891.41433532773231+x*
c;
3062 result =
ae_log(z, _state)+p;
3065 q = (x-0.5)*
ae_log(x, _state)-x+ls2pi;
3074 q = q+((7.9365079365079365079365*0.0001*p-2.7777777777777777777778*0.001)*p+0.0833333333333333333333)/
x;
3078 a = 8.11614167470508450300*0.0001;
3079 a = -5.95061904284301438324*0.0001+p*
a;
3080 a = 7.93650340457716943945*0.0001+p*
a;
3081 a = -2.77777777730099687205*0.001+p*
a;
3082 a = 8.33333333333331927722*0.01+p*
a;
3088 return _ialglib_i_lngamma(x, sgngam);
3093 static double gammafunc_gammastirf(
double x,
ae_state *_state)
3103 stir = 7.87311395793093628397E-4;
3104 stir = -2.29549961613378126380E-4+w*stir;
3105 stir = -2.68132617805781232825E-3+w*stir;
3106 stir = 3.47222221605458667310E-3+w*stir;
3107 stir = 8.33333333333482257126E-2+w*stir;
3112 v =
ae_pow(x, 0.5*x-0.25, _state);
3117 y =
ae_pow(x, x-0.5, _state)/
y;
3119 result = 2.50662827463100050242*y*
w;
3166 p = 0.007547728033418631287834;
3167 p = -0.288805137207594084924010+xsq*p;
3168 p = 14.3383842191748205576712+xsq*p;
3169 p = 38.0140318123903008244444+xsq*p;
3170 p = 3017.82788536507577809226+xsq*p;
3171 p = 7404.07142710151470082064+xsq*p;
3172 p = 80437.3630960840172832162+xsq*p;
3174 q = 1.00000000000000000000000+xsq*q;
3175 q = 38.0190713951939403753468+xsq*q;
3176 q = 658.070155459240506326937+xsq*q;
3177 q = 6379.60017324428279487120+xsq*q;
3178 q = 34216.5257924628539769006+xsq*q;
3179 q = 80437.3630960840172826266+xsq*q;
3180 result = s*1.1283791670955125738961589031*x*p/q;
3243 p = 0.5641877825507397413087057563+x*p;
3244 p = 9.675807882987265400604202961+x*p;
3245 p = 77.08161730368428609781633646+x*p;
3246 p = 368.5196154710010637133875746+x*p;
3247 p = 1143.262070703886173606073338+x*p;
3248 p = 2320.439590251635247384768711+x*p;
3249 p = 2898.0293292167655611275846+x*p;
3250 p = 1826.3348842295112592168999+x*p;
3252 q = 17.14980943627607849376131193+x*q;
3253 q = 137.1255960500622202878443578+x*q;
3254 q = 661.7361207107653469211984771+x*q;
3255 q = 2094.384367789539593790281779+x*q;
3256 q = 4429.612803883682726711528526+x*q;
3257 q = 6089.5424232724435504633068+x*q;
3258 q = 4958.82756472114071495438422+x*q;
3259 q = 1826.3348842295112595576438+x*q;
3300 result = 0.5*(
errorfunction(x/1.41421356237309504880, _state)+1);
3366 expm2 = 0.13533528323661269189;
3367 s2pi = 2.50662827463100050242;
3389 p0 = -59.9633501014107895267;
3390 p0 = 98.0010754185999661536+y2*p0;
3391 p0 = -56.6762857469070293439+y2*p0;
3392 p0 = 13.9312609387279679503+y2*p0;
3393 p0 = -1.23916583867381258016+y2*p0;
3395 q0 = 1.95448858338141759834+y2*
q0;
3396 q0 = 4.67627912898881538453+y2*
q0;
3397 q0 = 86.3602421390890590575+y2*
q0;
3398 q0 = -225.462687854119370527+y2*
q0;
3399 q0 = 200.260212380060660359+y2*
q0;
3400 q0 = -82.0372256168333339912+y2*
q0;
3401 q0 = 15.9056225126211695515+y2*
q0;
3402 q0 = -1.18331621121330003142+y2*
q0;
3413 p1 = 4.05544892305962419923;
3414 p1 = 31.5251094599893866154+z*p1;
3415 p1 = 57.1628192246421288162+z*p1;
3416 p1 = 44.0805073893200834700+z*p1;
3417 p1 = 14.6849561928858024014+z*p1;
3418 p1 = 2.18663306850790267539+z*p1;
3419 p1 = -1.40256079171354495875*0.1+z*p1;
3420 p1 = -3.50424626827848203418*0.01+z*p1;
3421 p1 = -8.57456785154685413611*0.0001+z*p1;
3423 q1 = 15.7799883256466749731+z*
q1;
3424 q1 = 45.3907635128879210584+z*
q1;
3425 q1 = 41.3172038254672030440+z*
q1;
3426 q1 = 15.0425385692907503408+z*
q1;
3427 q1 = 2.50464946208309415979+z*
q1;
3428 q1 = -1.42182922854787788574*0.1+z*
q1;
3429 q1 = -3.80806407691578277194*0.01+z*
q1;
3430 q1 = -9.33259480895457427372*0.0001+z*
q1;
3435 p2 = 3.23774891776946035970;
3436 p2 = 6.91522889068984211695+z*p2;
3437 p2 = 3.93881025292474443415+z*p2;
3438 p2 = 1.33303460815807542389+z*p2;
3439 p2 = 2.01485389549179081538*0.1+z*p2;
3440 p2 = 1.23716634817820021358*0.01+z*p2;
3441 p2 = 3.01581553508235416007*0.0001+z*p2;
3442 p2 = 2.65806974686737550832*0.000001+z*p2;
3443 p2 = 6.23974539184983293730*0.000000001+z*p2;
3445 q2 = 6.02427039364742014255+z*
q2;
3446 q2 = 3.67983563856160859403+z*
q2;
3447 q2 = 1.37702099489081330271+z*
q2;
3448 q2 = 2.16236993594496635890*0.1+z*
q2;
3449 q2 = 1.34204006088543189037*0.01+z*
q2;
3450 q2 = 3.28014464682127739104*0.0001+z*
q2;
3451 q2 = 2.89247864745380683936*0.000001+z*
q2;
3452 q2 = 6.79019408009981274425*0.000000001+z*
q2;
3498 double igammaepsilon;
3507 igammaepsilon = 0.000000000000001;
3575 double igammaepsilon;
3576 double igammabignumber;
3577 double igammabignumberinv;
3596 igammaepsilon = 0.000000000000001;
3597 igammabignumber = 4503599627370496.0;
3598 igammabignumberinv = 2.22044604925031308085*0.0000000000000001;
3630 pk = pkm1*z-pkm2*
yc;
3631 qk = qkm1*z-qkm2*
yc;
3635 t =
ae_fabs((ans-r)/r, _state);
3648 pkm2 = pkm2*igammabignumberinv;
3649 pkm1 = pkm1*igammabignumberinv;
3650 qkm2 = qkm2*igammabignumberinv;
3651 qkm1 = qkm1*igammabignumberinv;
3698 double igammaepsilon;
3699 double iinvgammabignumber;
3715 igammaepsilon = 0.000000000000001;
3716 iinvgammabignumber = 4503599627370496.0;
3717 x0 = iinvgammabignumber;
3721 dithresh = 5*igammaepsilon;
3725 lgm =
lngamma(a, &tmp, _state);
3750 d = (a-1)*
ae_log(x, _state)-x-lgm;
3766 if(
ae_fp_eq(x0,iinvgammabignumber) )
3772 while(
ae_fp_eq(x0,iinvgammabignumber))
3792 lgm = (x0-x1)/(x1+x0);
3823 d = (y0-yl)/(yh-yl);
3845 d = (y0-yl)/(yh-yl);
3935 sqpii = 5.64189583547756286948E-1;
3936 c1 = 0.35502805388781723926;
3937 c2 = 0.258819403792806798405;
3938 sqrt3 = 1.732050807568877293527;
3952 zeta = -2.0*x*t/3.0;
3957 afn = -1.31696323418331795333E-1;
3958 afn = afn*zz-6.26456544431912369773E-1;
3959 afn = afn*zz-6.93158036036933542233E-1;
3960 afn = afn*zz-2.79779981545119124951E-1;
3961 afn = afn*zz-4.91900132609500318020E-2;
3962 afn = afn*zz-4.06265923594885404393E-3;
3963 afn = afn*zz-1.59276496239262096340E-4;
3964 afn = afn*zz-2.77649108155232920844E-6;
3965 afn = afn*zz-1.67787698489114633780E-8;
3966 afd = 1.00000000000000000000E0;
3967 afd = afd*zz+1.33560420706553243746E1;
3968 afd = afd*zz+3.26825032795224613948E1;
3969 afd = afd*zz+2.67367040941499554804E1;
3970 afd = afd*zz+9.18707402907259625840E0;
3971 afd = afd*zz+1.47529146771666414581E0;
3972 afd = afd*zz+1.15687173795188044134E-1;
3973 afd = afd*zz+4.40291641615211203805E-3;
3974 afd = afd*zz+7.54720348287414296618E-5;
3975 afd = afd*zz+4.51850092970580378464E-7;
3976 uf = 1.0+zz*afn/afd;
3977 agn = 1.97339932091685679179E-2;
3978 agn = agn*zz+3.91103029615688277255E-1;
3979 agn = agn*zz+1.06579897599595591108E0;
3980 agn = agn*zz+9.39169229816650230044E-1;
3981 agn = agn*zz+3.51465656105547619242E-1;
3982 agn = agn*zz+6.33888919628925490927E-2;
3983 agn = agn*zz+5.85804113048388458567E-3;
3984 agn = agn*zz+2.82851600836737019778E-4;
3985 agn = agn*zz+6.98793669997260967291E-6;
3986 agn = agn*zz+8.11789239554389293311E-8;
3987 agn = agn*zz+3.41551784765923618484E-10;
3988 agd = 1.00000000000000000000E0;
3989 agd = agd*zz+9.30892908077441974853E0;
3990 agd = agd*zz+1.98352928718312140417E1;
3991 agd = agd*zz+1.55646628932864612953E1;
3992 agd = agd*zz+5.47686069422975497931E0;
3993 agd = agd*zz+9.54293611618961883998E-1;
3994 agd = agd*zz+8.64580826352392193095E-2;
3995 agd = agd*zz+4.12656523824222607191E-3;
3996 agd = agd*zz+1.01259085116509135510E-4;
3997 agd = agd*zz+1.17166733214413521882E-6;
3998 agd = agd*zz+4.91834570062930015649E-9;
4000 theta = zeta+0.25*
ae_pi;
4001 f =
ae_sin(theta, _state);
4002 g =
ae_cos(theta, _state);
4003 *ai = k*(f*uf-g*ug);
4004 *bi = k*(g*uf+f*ug);
4005 apfn = 1.85365624022535566142E-1;
4006 apfn = apfn*zz+8.86712188052584095637E-1;
4007 apfn = apfn*zz+9.87391981747398547272E-1;
4008 apfn = apfn*zz+4.01241082318003734092E-1;
4009 apfn = apfn*zz+7.10304926289631174579E-2;
4010 apfn = apfn*zz+5.90618657995661810071E-3;
4011 apfn = apfn*zz+2.33051409401776799569E-4;
4012 apfn = apfn*zz+4.08718778289035454598E-6;
4013 apfn = apfn*zz+2.48379932900442457853E-8;
4014 apfd = 1.00000000000000000000E0;
4015 apfd = apfd*zz+1.47345854687502542552E1;
4016 apfd = apfd*zz+3.75423933435489594466E1;
4017 apfd = apfd*zz+3.14657751203046424330E1;
4018 apfd = apfd*zz+1.09969125207298778536E1;
4019 apfd = apfd*zz+1.78885054766999417817E0;
4020 apfd = apfd*zz+1.41733275753662636873E-1;
4021 apfd = apfd*zz+5.44066067017226003627E-3;
4022 apfd = apfd*zz+9.39421290654511171663E-5;
4023 apfd = apfd*zz+5.65978713036027009243E-7;
4024 uf = 1.0+zz*apfn/apfd;
4025 apgn = -3.55615429033082288335E-2;
4026 apgn = apgn*zz-6.37311518129435504426E-1;
4027 apgn = apgn*zz-1.70856738884312371053E0;
4028 apgn = apgn*zz-1.50221872117316635393E0;
4029 apgn = apgn*zz-5.63606665822102676611E-1;
4030 apgn = apgn*zz-1.02101031120216891789E-1;
4031 apgn = apgn*zz-9.48396695961445269093E-3;
4032 apgn = apgn*zz-4.60325307486780994357E-4;
4033 apgn = apgn*zz-1.14300836484517375919E-5;
4034 apgn = apgn*zz-1.33415518685547420648E-7;
4035 apgn = apgn*zz-5.63803833958893494476E-10;
4036 apgd = 1.00000000000000000000E0;
4037 apgd = apgd*zz+9.85865801696130355144E0;
4038 apgd = apgd*zz+2.16401867356585941885E1;
4039 apgd = apgd*zz+1.73130776389749389525E1;
4040 apgd = apgd*zz+6.17872175280828766327E0;
4041 apgd = apgd*zz+1.08848694396321495475E0;
4042 apgd = apgd*zz+9.95005543440888479402E-2;
4043 apgd = apgd*zz+4.78468199683886610842E-3;
4044 apgd = apgd*zz+1.18159633322838625562E-4;
4045 apgd = apgd*zz+1.37480673554219441465E-6;
4046 apgd = apgd*zz+5.79912514929147598821E-9;
4049 *aip = -k*(g*uf+f*ug);
4050 *bip = k*(f*uf-g*ug);
4058 g =
ae_exp(zeta, _state);
4062 an = 3.46538101525629032477E-1;
4063 an = an*z+1.20075952739645805542E1;
4064 an = an*z+7.62796053615234516538E1;
4065 an = an*z+1.68089224934630576269E2;
4066 an = an*z+1.59756391350164413639E2;
4067 an = an*z+7.05360906840444183113E1;
4068 an = an*z+1.40264691163389668864E1;
4069 an = an*z+9.99999999999999995305E-1;
4070 ad = 5.67594532638770212846E-1;
4071 ad = ad*z+1.47562562584847203173E1;
4072 ad = ad*z+8.45138970141474626562E1;
4073 ad = ad*z+1.77318088145400459522E2;
4074 ad = ad*z+1.64234692871529701831E2;
4075 ad = ad*z+7.14778400825575695274E1;
4076 ad = ad*z+1.40959135607834029598E1;
4077 ad = ad*z+1.00000000000000000470E0;
4081 apn = 6.13759184814035759225E-1;
4082 apn = apn*z+1.47454670787755323881E1;
4083 apn = apn*z+8.20584123476060982430E1;
4084 apn = apn*z+1.71184781360976385540E2;
4085 apn = apn*z+1.59317847137141783523E2;
4086 apn = apn*z+6.99778599330103016170E1;
4087 apn = apn*z+1.39470856980481566958E1;
4088 apn = apn*z+1.00000000000000000550E0;
4089 apd = 3.34203677749736953049E-1;
4090 apd = apd*z+1.11810297306158156705E1;
4091 apd = apd*z+7.11727352147859965283E1;
4092 apd = apd*z+1.58778084372838313640E2;
4093 apd = apd*z+1.53206427475809220834E2;
4094 apd = apd*z+6.86752304592780337944E1;
4095 apd = apd*z+1.38498634758259442477E1;
4096 apd = apd*z+9.99999999999999994502E-1;
4101 bn16 = -2.53240795869364152689E-1;
4102 bn16 = bn16*z+5.75285167332467384228E-1;
4103 bn16 = bn16*z-3.29907036873225371650E-1;
4104 bn16 = bn16*z+6.44404068948199951727E-2;
4105 bn16 = bn16*z-3.82519546641336734394E-3;
4106 bd16 = 1.00000000000000000000E0;
4107 bd16 = bd16*z-7.15685095054035237902E0;
4108 bd16 = bd16*z+1.06039580715664694291E1;
4109 bd16 = bd16*z-5.23246636471251500874E0;
4110 bd16 = bd16*z+9.57395864378383833152E-1;
4111 bd16 = bd16*z-5.50828147163549611107E-2;
4115 bppn = 4.65461162774651610328E-1;
4116 bppn = bppn*z-1.08992173800493920734E0;
4117 bppn = bppn*z+6.38800117371827987759E-1;
4118 bppn = bppn*z-1.26844349553102907034E-1;
4119 bppn = bppn*z+7.62487844342109852105E-3;
4120 bppd = 1.00000000000000000000E0;
4121 bppd = bppd*z-8.70622787633159124240E0;
4122 bppd = bppd*z+1.38993162704553213172E1;
4123 bppd = bppd*z-7.14116144616431159572E0;
4124 bppd = bppd*z+1.34008595960680518666E0;
4125 bppd = bppd*z-7.84273211323341930448E-2;
4161 *bi = sqrt3*(uf+ug);
4193 *bip = sqrt3*(uf+ug);
4247 bessel_besselasympt0(x, &pzero, &qzero, _state);
4253 p1 = 26857.86856980014981415848441;
4254 p1 = -40504123.71833132706360663322+xsq*p1;
4255 p1 = 25071582855.36881945555156435+xsq*p1;
4256 p1 = -8085222034853.793871199468171+xsq*p1;
4257 p1 = 1434354939140344.111664316553+xsq*p1;
4258 p1 = -136762035308817138.6865416609+xsq*p1;
4259 p1 = 6382059341072356562.289432465+xsq*p1;
4260 p1 = -117915762910761053603.8440800+xsq*p1;
4261 p1 = 493378725179413356181.6813446+xsq*p1;
4263 q1 = 1363.063652328970604442810507+xsq*
q1;
4264 q1 = 1114636.098462985378182402543+xsq*
q1;
4265 q1 = 669998767.2982239671814028660+xsq*
q1;
4266 q1 = 312304311494.1213172572469442+xsq*
q1;
4267 q1 = 112775673967979.8507056031594+xsq*
q1;
4268 q1 = 30246356167094626.98627330784+xsq*
q1;
4269 q1 = 5428918384092285160.200195092+xsq*
q1;
4270 q1 = 493378725179413356211.3278438+xsq*
q1;
4315 bessel_besselasympt1(x, &pzero, &qzero, _state);
4325 p1 = 2701.122710892323414856790990;
4326 p1 = -4695753.530642995859767162166+xsq*p1;
4327 p1 = 3413234182.301700539091292655+xsq*p1;
4328 p1 = -1322983480332.126453125473247+xsq*p1;
4329 p1 = 290879526383477.5409737601689+xsq*p1;
4330 p1 = -35888175699101060.50743641413+xsq*p1;
4331 p1 = 2316433580634002297.931815435+xsq*p1;
4332 p1 = -66721065689249162980.20941484+xsq*p1;
4333 p1 = 581199354001606143928.050809+xsq*p1;
4335 q1 = 1606.931573481487801970916749+xsq*
q1;
4336 q1 = 1501793.594998585505921097578+xsq*
q1;
4337 q1 = 1013863514.358673989967045588+xsq*
q1;
4338 q1 = 524371026216.7649715406728642+xsq*
q1;
4339 q1 = 208166122130760.7351240184229+xsq*
q1;
4340 q1 = 60920613989175217.46105196863+xsq*
q1;
4341 q1 = 11857707121903209998.37113348+xsq*
q1;
4342 q1 = 1162398708003212287858.529400+xsq*
q1;
4457 pkm2 = (pkm1*r-pk*
x)/x;
4518 bessel_besselasympt0(x, &pzero, &qzero, _state);
4524 p4 = -41370.35497933148554125235152;
4525 p4 = 59152134.65686889654273830069+xsq*p4;
4526 p4 = -34363712229.79040378171030138+xsq*p4;
4527 p4 = 10255208596863.94284509167421+xsq*p4;
4528 p4 = -1648605817185729.473122082537+xsq*p4;
4529 p4 = 137562431639934407.8571335453+xsq*p4;
4530 p4 = -5247065581112764941.297350814+xsq*p4;
4531 p4 = 65874732757195549259.99402049+xsq*p4;
4532 p4 = -27502866786291095837.01933175+xsq*p4;
4534 q4 = 1282.452772478993804176329391+xsq*q4;
4535 q4 = 1001702.641288906265666651753+xsq*q4;
4536 q4 = 579512264.0700729537480087915+xsq*q4;
4537 q4 = 261306575504.1081249568482092+xsq*q4;
4538 q4 = 91620380340751.85262489147968+xsq*q4;
4539 q4 = 23928830434997818.57439356652+xsq*q4;
4540 q4 = 4192417043410839973.904769661+xsq*q4;
4541 q4 = 372645883898616588198.9980+xsq*q4;
4581 bessel_besselasympt1(x, &pzero, &qzero, _state);
4587 p4 = -2108847.540133123652824139923;
4588 p4 = 3639488548.124002058278999428+xsq*p4;
4589 p4 = -2580681702194.450950541426399+xsq*p4;
4590 p4 = 956993023992168.3481121552788+xsq*p4;
4591 p4 = -196588746272214065.8820322248+xsq*p4;
4592 p4 = 21931073399177975921.11427556+xsq*p4;
4593 p4 = -1212297555414509577913.561535+xsq*p4;
4594 p4 = 26554738314348543268942.48968+xsq*p4;
4595 p4 = -99637534243069222259967.44354+xsq*p4;
4597 q4 = 1612.361029677000859332072312+xsq*q4;
4598 q4 = 1563282.754899580604737366452+xsq*q4;
4599 q4 = 1128686837.169442121732366891+xsq*q4;
4600 q4 = 646534088126.5275571961681500+xsq*q4;
4601 q4 = 297663212564727.6729292742282+xsq*q4;
4602 q4 = 108225825940881955.2553850180+xsq*q4;
4603 q4 = 29549879358971486742.90758119+xsq*q4;
4604 q4 = 5435310377188854170800.653097+xsq*q4;
4605 q4 = 508206736694124324531442.4152+xsq*q4;
4664 for(i=1; i<=n-1; i++)
4714 bessel_besselmfirstcheb(-4.41534164647933937950E-18, &b0, &b1, &b2, _state);
4715 bessel_besselmnextcheb(y, 3.33079451882223809783E-17, &b0, &b1, &b2, _state);
4716 bessel_besselmnextcheb(y, -2.43127984654795469359E-16, &b0, &b1, &b2, _state);
4717 bessel_besselmnextcheb(y, 1.71539128555513303061E-15, &b0, &b1, &b2, _state);
4718 bessel_besselmnextcheb(y, -1.16853328779934516808E-14, &b0, &b1, &b2, _state);
4719 bessel_besselmnextcheb(y, 7.67618549860493561688E-14, &b0, &b1, &b2, _state);
4720 bessel_besselmnextcheb(y, -4.85644678311192946090E-13, &b0, &b1, &b2, _state);
4721 bessel_besselmnextcheb(y, 2.95505266312963983461E-12, &b0, &b1, &b2, _state);
4722 bessel_besselmnextcheb(y, -1.72682629144155570723E-11, &b0, &b1, &b2, _state);
4723 bessel_besselmnextcheb(y, 9.67580903537323691224E-11, &b0, &b1, &b2, _state);
4724 bessel_besselmnextcheb(y, -5.18979560163526290666E-10, &b0, &b1, &b2, _state);
4725 bessel_besselmnextcheb(y, 2.65982372468238665035E-9, &b0, &b1, &b2, _state);
4726 bessel_besselmnextcheb(y, -1.30002500998624804212E-8, &b0, &b1, &b2, _state);
4727 bessel_besselmnextcheb(y, 6.04699502254191894932E-8, &b0, &b1, &b2, _state);
4728 bessel_besselmnextcheb(y, -2.67079385394061173391E-7, &b0, &b1, &b2, _state);
4729 bessel_besselmnextcheb(y, 1.11738753912010371815E-6, &b0, &b1, &b2, _state);
4730 bessel_besselmnextcheb(y, -4.41673835845875056359E-6, &b0, &b1, &b2, _state);
4731 bessel_besselmnextcheb(y, 1.64484480707288970893E-5, &b0, &b1, &b2, _state);
4732 bessel_besselmnextcheb(y, -5.75419501008210370398E-5, &b0, &b1, &b2, _state);
4733 bessel_besselmnextcheb(y, 1.88502885095841655729E-4, &b0, &b1, &b2, _state);
4734 bessel_besselmnextcheb(y, -5.76375574538582365885E-4, &b0, &b1, &b2, _state);
4735 bessel_besselmnextcheb(y, 1.63947561694133579842E-3, &b0, &b1, &b2, _state);
4736 bessel_besselmnextcheb(y, -4.32430999505057594430E-3, &b0, &b1, &b2, _state);
4737 bessel_besselmnextcheb(y, 1.05464603945949983183E-2, &b0, &b1, &b2, _state);
4738 bessel_besselmnextcheb(y, -2.37374148058994688156E-2, &b0, &b1, &b2, _state);
4739 bessel_besselmnextcheb(y, 4.93052842396707084878E-2, &b0, &b1, &b2, _state);
4740 bessel_besselmnextcheb(y, -9.49010970480476444210E-2, &b0, &b1, &b2, _state);
4741 bessel_besselmnextcheb(y, 1.71620901522208775349E-1, &b0, &b1, &b2, _state);
4742 bessel_besselmnextcheb(y, -3.04682672343198398683E-1, &b0, &b1, &b2, _state);
4743 bessel_besselmnextcheb(y, 6.76795274409476084995E-1, &b0, &b1, &b2, _state);
4745 result =
ae_exp(x, _state)*v;
4749 bessel_besselmfirstcheb(-7.23318048787475395456E-18, &b0, &b1, &b2, _state);
4750 bessel_besselmnextcheb(z, -4.83050448594418207126E-18, &b0, &b1, &b2, _state);
4751 bessel_besselmnextcheb(z, 4.46562142029675999901E-17, &b0, &b1, &b2, _state);
4752 bessel_besselmnextcheb(z, 3.46122286769746109310E-17, &b0, &b1, &b2, _state);
4753 bessel_besselmnextcheb(z, -2.82762398051658348494E-16, &b0, &b1, &b2, _state);
4754 bessel_besselmnextcheb(z, -3.42548561967721913462E-16, &b0, &b1, &b2, _state);
4755 bessel_besselmnextcheb(z, 1.77256013305652638360E-15, &b0, &b1, &b2, _state);
4756 bessel_besselmnextcheb(z, 3.81168066935262242075E-15, &b0, &b1, &b2, _state);
4757 bessel_besselmnextcheb(z, -9.55484669882830764870E-15, &b0, &b1, &b2, _state);
4758 bessel_besselmnextcheb(z, -4.15056934728722208663E-14, &b0, &b1, &b2, _state);
4759 bessel_besselmnextcheb(z, 1.54008621752140982691E-14, &b0, &b1, &b2, _state);
4760 bessel_besselmnextcheb(z, 3.85277838274214270114E-13, &b0, &b1, &b2, _state);
4761 bessel_besselmnextcheb(z, 7.18012445138366623367E-13, &b0, &b1, &b2, _state);
4762 bessel_besselmnextcheb(z, -1.79417853150680611778E-12, &b0, &b1, &b2, _state);
4763 bessel_besselmnextcheb(z, -1.32158118404477131188E-11, &b0, &b1, &b2, _state);
4764 bessel_besselmnextcheb(z, -3.14991652796324136454E-11, &b0, &b1, &b2, _state);
4765 bessel_besselmnextcheb(z, 1.18891471078464383424E-11, &b0, &b1, &b2, _state);
4766 bessel_besselmnextcheb(z, 4.94060238822496958910E-10, &b0, &b1, &b2, _state);
4767 bessel_besselmnextcheb(z, 3.39623202570838634515E-9, &b0, &b1, &b2, _state);
4768 bessel_besselmnextcheb(z, 2.26666899049817806459E-8, &b0, &b1, &b2, _state);
4769 bessel_besselmnextcheb(z, 2.04891858946906374183E-7, &b0, &b1, &b2, _state);
4770 bessel_besselmnextcheb(z, 2.89137052083475648297E-6, &b0, &b1, &b2, _state);
4771 bessel_besselmnextcheb(z, 6.88975834691682398426E-5, &b0, &b1, &b2, _state);
4772 bessel_besselmnextcheb(z, 3.36911647825569408990E-3, &b0, &b1, &b2, _state);
4773 bessel_besselmnextcheb(z, 8.04490411014108831608E-1, &b0, &b1, &b2, _state);
4816 bessel_besselm1firstcheb(2.77791411276104639959E-18, &b0, &b1, &b2, _state);
4817 bessel_besselm1nextcheb(y, -2.11142121435816608115E-17, &b0, &b1, &b2, _state);
4818 bessel_besselm1nextcheb(y, 1.55363195773620046921E-16, &b0, &b1, &b2, _state);
4819 bessel_besselm1nextcheb(y, -1.10559694773538630805E-15, &b0, &b1, &b2, _state);
4820 bessel_besselm1nextcheb(y, 7.60068429473540693410E-15, &b0, &b1, &b2, _state);
4821 bessel_besselm1nextcheb(y, -5.04218550472791168711E-14, &b0, &b1, &b2, _state);
4822 bessel_besselm1nextcheb(y, 3.22379336594557470981E-13, &b0, &b1, &b2, _state);
4823 bessel_besselm1nextcheb(y, -1.98397439776494371520E-12, &b0, &b1, &b2, _state);
4824 bessel_besselm1nextcheb(y, 1.17361862988909016308E-11, &b0, &b1, &b2, _state);
4825 bessel_besselm1nextcheb(y, -6.66348972350202774223E-11, &b0, &b1, &b2, _state);
4826 bessel_besselm1nextcheb(y, 3.62559028155211703701E-10, &b0, &b1, &b2, _state);
4827 bessel_besselm1nextcheb(y, -1.88724975172282928790E-9, &b0, &b1, &b2, _state);
4828 bessel_besselm1nextcheb(y, 9.38153738649577178388E-9, &b0, &b1, &b2, _state);
4829 bessel_besselm1nextcheb(y, -4.44505912879632808065E-8, &b0, &b1, &b2, _state);
4830 bessel_besselm1nextcheb(y, 2.00329475355213526229E-7, &b0, &b1, &b2, _state);
4831 bessel_besselm1nextcheb(y, -8.56872026469545474066E-7, &b0, &b1, &b2, _state);
4832 bessel_besselm1nextcheb(y, 3.47025130813767847674E-6, &b0, &b1, &b2, _state);
4833 bessel_besselm1nextcheb(y, -1.32731636560394358279E-5, &b0, &b1, &b2, _state);
4834 bessel_besselm1nextcheb(y, 4.78156510755005422638E-5, &b0, &b1, &b2, _state);
4835 bessel_besselm1nextcheb(y, -1.61760815825896745588E-4, &b0, &b1, &b2, _state);
4836 bessel_besselm1nextcheb(y, 5.12285956168575772895E-4, &b0, &b1, &b2, _state);
4837 bessel_besselm1nextcheb(y, -1.51357245063125314899E-3, &b0, &b1, &b2, _state);
4838 bessel_besselm1nextcheb(y, 4.15642294431288815669E-3, &b0, &b1, &b2, _state);
4839 bessel_besselm1nextcheb(y, -1.05640848946261981558E-2, &b0, &b1, &b2, _state);
4840 bessel_besselm1nextcheb(y, 2.47264490306265168283E-2, &b0, &b1, &b2, _state);
4841 bessel_besselm1nextcheb(y, -5.29459812080949914269E-2, &b0, &b1, &b2, _state);
4842 bessel_besselm1nextcheb(y, 1.02643658689847095384E-1, &b0, &b1, &b2, _state);
4843 bessel_besselm1nextcheb(y, -1.76416518357834055153E-1, &b0, &b1, &b2, _state);
4844 bessel_besselm1nextcheb(y, 2.52587186443633654823E-1, &b0, &b1, &b2, _state);
4846 z = v*z*
ae_exp(z, _state);
4851 bessel_besselm1firstcheb(7.51729631084210481353E-18, &b0, &b1, &b2, _state);
4852 bessel_besselm1nextcheb(y, 4.41434832307170791151E-18, &b0, &b1, &b2, _state);
4853 bessel_besselm1nextcheb(y, -4.65030536848935832153E-17, &b0, &b1, &b2, _state);
4854 bessel_besselm1nextcheb(y, -3.20952592199342395980E-17, &b0, &b1, &b2, _state);
4855 bessel_besselm1nextcheb(y, 2.96262899764595013876E-16, &b0, &b1, &b2, _state);
4856 bessel_besselm1nextcheb(y, 3.30820231092092828324E-16, &b0, &b1, &b2, _state);
4857 bessel_besselm1nextcheb(y, -1.88035477551078244854E-15, &b0, &b1, &b2, _state);
4858 bessel_besselm1nextcheb(y, -3.81440307243700780478E-15, &b0, &b1, &b2, _state);
4859 bessel_besselm1nextcheb(y, 1.04202769841288027642E-14, &b0, &b1, &b2, _state);
4860 bessel_besselm1nextcheb(y, 4.27244001671195135429E-14, &b0, &b1, &b2, _state);
4861 bessel_besselm1nextcheb(y, -2.10154184277266431302E-14, &b0, &b1, &b2, _state);
4862 bessel_besselm1nextcheb(y, -4.08355111109219731823E-13, &b0, &b1, &b2, _state);
4863 bessel_besselm1nextcheb(y, -7.19855177624590851209E-13, &b0, &b1, &b2, _state);
4864 bessel_besselm1nextcheb(y, 2.03562854414708950722E-12, &b0, &b1, &b2, _state);
4865 bessel_besselm1nextcheb(y, 1.41258074366137813316E-11, &b0, &b1, &b2, _state);
4866 bessel_besselm1nextcheb(y, 3.25260358301548823856E-11, &b0, &b1, &b2, _state);
4867 bessel_besselm1nextcheb(y, -1.89749581235054123450E-11, &b0, &b1, &b2, _state);
4868 bessel_besselm1nextcheb(y, -5.58974346219658380687E-10, &b0, &b1, &b2, _state);
4869 bessel_besselm1nextcheb(y, -3.83538038596423702205E-9, &b0, &b1, &b2, _state);
4870 bessel_besselm1nextcheb(y, -2.63146884688951950684E-8, &b0, &b1, &b2, _state);
4871 bessel_besselm1nextcheb(y, -2.51223623787020892529E-7, &b0, &b1, &b2, _state);
4872 bessel_besselm1nextcheb(y, -3.88256480887769039346E-6, &b0, &b1, &b2, _state);
4873 bessel_besselm1nextcheb(y, -1.10588938762623716291E-4, &b0, &b1, &b2, _state);
4874 bessel_besselm1nextcheb(y, -9.76109749136146840777E-3, &b0, &b1, &b2, _state);
4875 bessel_besselm1nextcheb(y, 7.78576235018280120474E-1, &b0, &b1, &b2, _state);
4924 bessel_besselmfirstcheb(1.37446543561352307156E-16, &b0, &b1, &b2, _state);
4925 bessel_besselmnextcheb(y, 4.25981614279661018399E-14, &b0, &b1, &b2, _state);
4926 bessel_besselmnextcheb(y, 1.03496952576338420167E-11, &b0, &b1, &b2, _state);
4927 bessel_besselmnextcheb(y, 1.90451637722020886025E-9, &b0, &b1, &b2, _state);
4928 bessel_besselmnextcheb(y, 2.53479107902614945675E-7, &b0, &b1, &b2, _state);
4929 bessel_besselmnextcheb(y, 2.28621210311945178607E-5, &b0, &b1, &b2, _state);
4930 bessel_besselmnextcheb(y, 1.26461541144692592338E-3, &b0, &b1, &b2, _state);
4931 bessel_besselmnextcheb(y, 3.59799365153615016266E-2, &b0, &b1, &b2, _state);
4932 bessel_besselmnextcheb(y, 3.44289899924628486886E-1, &b0, &b1, &b2, _state);
4933 bessel_besselmnextcheb(y, -5.35327393233902768720E-1, &b0, &b1, &b2, _state);
4940 bessel_besselmfirstcheb(5.30043377268626276149E-18, &b0, &b1, &b2, _state);
4941 bessel_besselmnextcheb(z, -1.64758043015242134646E-17, &b0, &b1, &b2, _state);
4942 bessel_besselmnextcheb(z, 5.21039150503902756861E-17, &b0, &b1, &b2, _state);
4943 bessel_besselmnextcheb(z, -1.67823109680541210385E-16, &b0, &b1, &b2, _state);
4944 bessel_besselmnextcheb(z, 5.51205597852431940784E-16, &b0, &b1, &b2, _state);
4945 bessel_besselmnextcheb(z, -1.84859337734377901440E-15, &b0, &b1, &b2, _state);
4946 bessel_besselmnextcheb(z, 6.34007647740507060557E-15, &b0, &b1, &b2, _state);
4947 bessel_besselmnextcheb(z, -2.22751332699166985548E-14, &b0, &b1, &b2, _state);
4948 bessel_besselmnextcheb(z, 8.03289077536357521100E-14, &b0, &b1, &b2, _state);
4949 bessel_besselmnextcheb(z, -2.98009692317273043925E-13, &b0, &b1, &b2, _state);
4950 bessel_besselmnextcheb(z, 1.14034058820847496303E-12, &b0, &b1, &b2, _state);
4951 bessel_besselmnextcheb(z, -4.51459788337394416547E-12, &b0, &b1, &b2, _state);
4952 bessel_besselmnextcheb(z, 1.85594911495471785253E-11, &b0, &b1, &b2, _state);
4953 bessel_besselmnextcheb(z, -7.95748924447710747776E-11, &b0, &b1, &b2, _state);
4954 bessel_besselmnextcheb(z, 3.57739728140030116597E-10, &b0, &b1, &b2, _state);
4955 bessel_besselmnextcheb(z, -1.69753450938905987466E-9, &b0, &b1, &b2, _state);
4956 bessel_besselmnextcheb(z, 8.57403401741422608519E-9, &b0, &b1, &b2, _state);
4957 bessel_besselmnextcheb(z, -4.66048989768794782956E-8, &b0, &b1, &b2, _state);
4958 bessel_besselmnextcheb(z, 2.76681363944501510342E-7, &b0, &b1, &b2, _state);
4959 bessel_besselmnextcheb(z, -1.83175552271911948767E-6, &b0, &b1, &b2, _state);
4960 bessel_besselmnextcheb(z, 1.39498137188764993662E-5, &b0, &b1, &b2, _state);
4961 bessel_besselmnextcheb(z, -1.28495495816278026384E-4, &b0, &b1, &b2, _state);
4962 bessel_besselmnextcheb(z, 1.56988388573005337491E-3, &b0, &b1, &b2, _state);
4963 bessel_besselmnextcheb(z, -3.14481013119645005427E-2, &b0, &b1, &b2, _state);
4964 bessel_besselmnextcheb(z, 2.44030308206595545468E0, &b0, &b1, &b2, _state);
5008 bessel_besselm1firstcheb(-7.02386347938628759343E-18, &b0, &b1, &b2, _state);
5009 bessel_besselm1nextcheb(y, -2.42744985051936593393E-15, &b0, &b1, &b2, _state);
5010 bessel_besselm1nextcheb(y, -6.66690169419932900609E-13, &b0, &b1, &b2, _state);
5011 bessel_besselm1nextcheb(y, -1.41148839263352776110E-10, &b0, &b1, &b2, _state);
5012 bessel_besselm1nextcheb(y, -2.21338763073472585583E-8, &b0, &b1, &b2, _state);
5013 bessel_besselm1nextcheb(y, -2.43340614156596823496E-6, &b0, &b1, &b2, _state);
5014 bessel_besselm1nextcheb(y, -1.73028895751305206302E-4, &b0, &b1, &b2, _state);
5015 bessel_besselm1nextcheb(y, -6.97572385963986435018E-3, &b0, &b1, &b2, _state);
5016 bessel_besselm1nextcheb(y, -1.22611180822657148235E-1, &b0, &b1, &b2, _state);
5017 bessel_besselm1nextcheb(y, -3.53155960776544875667E-1, &b0, &b1, &b2, _state);
5018 bessel_besselm1nextcheb(y, 1.52530022733894777053E0, &b0, &b1, &b2, _state);
5025 bessel_besselm1firstcheb(-5.75674448366501715755E-18, &b0, &b1, &b2, _state);
5026 bessel_besselm1nextcheb(y, 1.79405087314755922667E-17, &b0, &b1, &b2, _state);
5027 bessel_besselm1nextcheb(y, -5.68946255844285935196E-17, &b0, &b1, &b2, _state);
5028 bessel_besselm1nextcheb(y, 1.83809354436663880070E-16, &b0, &b1, &b2, _state);
5029 bessel_besselm1nextcheb(y, -6.05704724837331885336E-16, &b0, &b1, &b2, _state);
5030 bessel_besselm1nextcheb(y, 2.03870316562433424052E-15, &b0, &b1, &b2, _state);
5031 bessel_besselm1nextcheb(y, -7.01983709041831346144E-15, &b0, &b1, &b2, _state);
5032 bessel_besselm1nextcheb(y, 2.47715442448130437068E-14, &b0, &b1, &b2, _state);
5033 bessel_besselm1nextcheb(y, -8.97670518232499435011E-14, &b0, &b1, &b2, _state);
5034 bessel_besselm1nextcheb(y, 3.34841966607842919884E-13, &b0, &b1, &b2, _state);
5035 bessel_besselm1nextcheb(y, -1.28917396095102890680E-12, &b0, &b1, &b2, _state);
5036 bessel_besselm1nextcheb(y, 5.13963967348173025100E-12, &b0, &b1, &b2, _state);
5037 bessel_besselm1nextcheb(y, -2.12996783842756842877E-11, &b0, &b1, &b2, _state);
5038 bessel_besselm1nextcheb(y, 9.21831518760500529508E-11, &b0, &b1, &b2, _state);
5039 bessel_besselm1nextcheb(y, -4.19035475934189648750E-10, &b0, &b1, &b2, _state);
5040 bessel_besselm1nextcheb(y, 2.01504975519703286596E-9, &b0, &b1, &b2, _state);
5041 bessel_besselm1nextcheb(y, -1.03457624656780970260E-8, &b0, &b1, &b2, _state);
5042 bessel_besselm1nextcheb(y, 5.74108412545004946722E-8, &b0, &b1, &b2, _state);
5043 bessel_besselm1nextcheb(y, -3.50196060308781257119E-7, &b0, &b1, &b2, _state);
5044 bessel_besselm1nextcheb(y, 2.40648494783721712015E-6, &b0, &b1, &b2, _state);
5045 bessel_besselm1nextcheb(y, -1.93619797416608296024E-5, &b0, &b1, &b2, _state);
5046 bessel_besselm1nextcheb(y, 1.95215518471351631108E-4, &b0, &b1, &b2, _state);
5047 bessel_besselm1nextcheb(y, -2.85781685962277938680E-3, &b0, &b1, &b2, _state);
5048 bessel_besselm1nextcheb(y, 1.03923736576817238437E-1, &b0, &b1, &b2, _state);
5049 bessel_besselm1nextcheb(y, 2.72062619048444266945E0, &b0, &b1, &b2, _state);
5103 eul = 5.772156649015328606065e-1;
5112 ae_assert(n<=31,
"Overflow in BesselKN", _state);
5126 for(i=1; i<=n-1; i++)
5144 for(i=1; i<=n-1; i++)
5162 tlg = 2.0*
ae_log(0.5*x, _state);
5178 t = t*(z0/(k*(k+
n)));
5181 s = s+(pk+pn-tlg)*t;
5235 static void bessel_besselmfirstcheb(
double c,
5255 static void bessel_besselmnextcheb(
double x,
5266 *b0 = x*(*b1)-(*b2)+
c;
5276 static void bessel_besselm1firstcheb(
double c,
5296 static void bessel_besselm1nextcheb(
double x,
5307 *b0 = x*(*b1)-(*b2)+
c;
5311 static void bessel_besselasympt0(
double x,
5327 p2 = 2485.271928957404011288128951+xsq*p2;
5328 p2 = 153982.6532623911470917825993+xsq*p2;
5329 p2 = 2016135.283049983642487182349+xsq*p2;
5330 p2 = 8413041.456550439208464315611+xsq*p2;
5331 p2 = 12332384.76817638145232406055+xsq*p2;
5332 p2 = 5393485.083869438325262122897+xsq*p2;
5334 q2 = 2615.700736920839685159081813+xsq*
q2;
5335 q2 = 156001.7276940030940592769933+xsq*
q2;
5336 q2 = 2025066.801570134013891035236+xsq*
q2;
5337 q2 = 8426449.050629797331554404810+xsq*
q2;
5338 q2 = 12338310.22786324960844856182+xsq*
q2;
5339 q2 = 5393485.083869438325560444960+xsq*
q2;
5341 p3 = -4.887199395841261531199129300+xsq*p3;
5342 p3 = -226.2630641933704113967255053+xsq*p3;
5343 p3 = -2365.956170779108192723612816+xsq*p3;
5344 p3 = -8239.066313485606568803548860+xsq*p3;
5345 p3 = -10381.41698748464093880530341+xsq*p3;
5346 p3 = -3984.617357595222463506790588+xsq*p3;
5348 q3 = 408.7714673983499223402830260+xsq*
q3;
5349 q3 = 15704.89191515395519392882766+xsq*
q3;
5350 q3 = 156021.3206679291652539287109+xsq*
q3;
5351 q3 = 533291.3634216897168722255057+xsq*
q3;
5352 q3 = 666745.4239319826986004038103+xsq*
q3;
5353 q3 = 255015.5108860942382983170882+xsq*
q3;
5359 static void bessel_besselasympt1(
double x,
5374 p2 = -1611.616644324610116477412898;
5375 p2 = -109824.0554345934672737413139+xsq*p2;
5376 p2 = -1523529.351181137383255105722+xsq*p2;
5377 p2 = -6603373.248364939109255245434+xsq*p2;
5378 p2 = -9942246.505077641195658377899+xsq*p2;
5379 p2 = -4435757.816794127857114720794+xsq*p2;
5381 q2 = -1455.009440190496182453565068+xsq*
q2;
5382 q2 = -107263.8599110382011903063867+xsq*
q2;
5383 q2 = -1511809.506634160881644546358+xsq*
q2;
5384 q2 = -6585339.479723087072826915069+xsq*
q2;
5385 q2 = -9934124.389934585658967556309+xsq*
q2;
5386 q2 = -4435757.816794127856828016962+xsq*
q2;
5387 p3 = 35.26513384663603218592175580;
5388 p3 = 1706.375429020768002061283546+xsq*p3;
5389 p3 = 18494.26287322386679652009819+xsq*p3;
5390 p3 = 66178.83658127083517939992166+xsq*p3;
5391 p3 = 85145.16067533570196555001171+xsq*p3;
5392 p3 = 33220.91340985722351859704442+xsq*p3;
5394 q3 = 863.8367769604990967475517183+xsq*
q3;
5395 q3 = 37890.22974577220264142952256+xsq*
q3;
5396 q3 = 400294.4358226697511708610813+xsq*
q3;
5397 q3 = 1419460.669603720892855755253+xsq*
q3;
5398 q3 = 1819458.042243997298924553839+xsq*
q3;
5399 q3 = 708712.8194102874357377502472+xsq*
q3;
5450 result = sg*
ae_exp(y, _state);
5529 big = 4.503599627370496e15;
5530 biginv = 2.22044604925031308085e-16;
5531 maxgam = 171.624376956302725;
5549 result = ibetaf_incompletebetaps(a, b, x, maxgam, _state);
5568 t = ibetaf_incompletebetaps(a, b, x, maxgam, _state);
5579 y = x*(a+b-2.0)-(a-1.0);
5582 w = ibetaf_incompletebetafe(a, b, x, big, biginv, _state);
5586 w = ibetaf_incompletebetafe2(a, b, x, big, biginv, _state)/
xc;
5589 t = b*
ae_log(xc, _state);
5592 t =
ae_pow(xc, b, _state);
5593 t = t*
ae_pow(x, a, _state);
5615 y = y+
ae_log(w/a, _state);
5746 breakihalvecycle = 6;
5757 if( mainlooppos==0 )
5768 mainlooppos = ihalve;
5791 lgm = (yp*yp-3.0)/6.0;
5792 x = 2.0/(1.0/(2.0*aaa-1.0)+1.0/(2.0*bbb-1.0));
5793 d = yp*
ae_sqrt(x+lgm, _state)/x-(1.0/(2.0*bbb-1.0)-1.0/(2.0*aaa-1.0))*(lgm+5.0/6.0-2.0/(3.0*x));
5800 x = aaa/(aaa+bbb*
ae_exp(d, _state));
5808 mainlooppos = ihalve;
5815 if( mainlooppos==ihalve )
5820 mainlooppos = ihalvecycle;
5827 if( mainlooppos==ihalvecycle )
5848 yp = (x1-
x0)/(x1+x0);
5874 di = 1.0-(1.0-
di)*(1.0-di);
5884 di = (y0-yyy)/(yh-yl);
5911 mainlooppos = ihalve;
5943 di = (yyy-
y0)/(yh-yl);
5950 mainlooppos = ihalvecycle;
5955 mainlooppos = breakihalvecycle;
5963 if( mainlooppos==breakihalvecycle )
5982 if( mainlooppos==newt )
5991 mainlooppos = newtcycle;
5998 if( mainlooppos==newtcycle )
6034 mainlooppos = breaknewtcycle;
6037 d = (aaa-1.0)*
ae_log(x, _state)+(bbb-1.0)*
ae_log(1.0-x, _state)+lgm;
6044 mainlooppos = breaknewtcycle;
6052 yyy = (x-
x0)/(x1-x0);
6053 xt = x0+0.5*yyy*(x-
x0);
6056 mainlooppos = breaknewtcycle;
6062 yyy = (x1-
x)/(x1-x0);
6063 xt = x1-0.5*yyy*(x1-
x);
6066 mainlooppos = breaknewtcycle;
6076 mainlooppos = newtcycle;
6081 mainlooppos = breaknewtcycle;
6089 if( mainlooppos==breaknewtcycle )
6092 mainlooppos = ihalve;
6122 static double ibetaf_incompletebetafe(
double a,
6170 xk = -x*k1*k2/(k3*k4);
6177 xk = x*k5*k6/(k7*k8);
6190 t =
ae_fabs((ans-r)/r, _state);
6238 static double ibetaf_incompletebetafe2(
double a,
6288 xk = -z*k1*k2/(k3*k4);
6295 xk = z*k5*k6/(k7*k8);
6308 t =
ae_fabs((ans-r)/r, _state);
6356 static double ibetaf_incompletebetaps(
double a,
6396 s = s*t*
ae_pow(x, a, _state);
6459 ae_assert(k>=-1&&k<=n,
"Domain error in BinomialDistribution", _state);
6473 dk =
ae_pow(1.0-p, dn, _state);
6529 ae_assert(k>=-1&&k<=n,
"Domain error in BinomialDistributionC", _state);
6549 dk = 1.0-
ae_pow(1.0-p, dn, _state);
6601 ae_assert(k>=0&&k<n,
"Domain error in InvBinomialDistribution", _state);
6611 p = 1.0-
ae_pow(y, 1.0/dn, _state);
6778 for(i=0; i<=n/2-1; i++)
7028 an = 1.13681498971755972054E-11;
7029 an = an*x2+8.49262267667473811108E-10;
7030 an = an*x2+1.94434204175553054283E-8;
7031 an = an*x2+9.53151741254484363489E-7;
7032 an = an*x2+3.07828309874913200438E-6;
7033 an = an*x2+3.52513368520288738649E-4;
7034 an = an*x2+(-8.50149846724410912031E-4);
7035 an = an*x2+4.22618223005546594270E-2;
7036 an = an*x2+(-9.17480371773452345351E-2);
7037 an = an*x2+9.99999999999999994612E-1;
7038 ad = 2.40372073066762605484E-11;
7039 ad = ad*x2+1.48864681368493396752E-9;
7040 ad = ad*x2+5.21265281010541664570E-8;
7041 ad = ad*x2+1.27258478273186970203E-6;
7042 ad = ad*x2+2.32490249820789513991E-5;
7043 ad = ad*x2+3.25524741826057911661E-4;
7044 ad = ad*x2+3.48805814657162590916E-3;
7045 ad = ad*x2+2.79448531198828973716E-2;
7046 ad = ad*x2+1.58874241960120565368E-1;
7047 ad = ad*x2+5.74918629489320327824E-1;
7048 ad = ad*x2+1.00000000000000000539E0;
7056 bn = 5.08955156417900903354E-1;
7057 bn = bn*x2-2.44754418142697847934E-1;
7058 bn = bn*x2+9.41512335303534411857E-2;
7059 bn = bn*x2-2.18711255142039025206E-2;
7060 bn = bn*x2+3.66207612329569181322E-3;
7061 bn = bn*x2-4.23209114460388756528E-4;
7062 bn = bn*x2+3.59641304793896631888E-5;
7063 bn = bn*x2-2.14640351719968974225E-6;
7064 bn = bn*x2+9.10010780076391431042E-8;
7065 bn = bn*x2-2.40274520828250956942E-9;
7066 bn = bn*x2+3.59233385440928410398E-11;
7067 bd = 1.00000000000000000000E0;
7068 bd = bd*x2-6.31839869873368190192E-1;
7069 bd = bd*x2+2.36706788228248691528E-1;
7070 bd = bd*x2-5.31806367003223277662E-2;
7071 bd = bd*x2+8.48041718586295374409E-3;
7072 bd = bd*x2-9.47996768486665330168E-4;
7073 bd = bd*x2+7.81025592944552338085E-5;
7074 bd = bd*x2-4.55875153252442634831E-6;
7075 bd = bd*x2+1.89100358111421846170E-7;
7076 bd = bd*x2-4.91324691331920606875E-9;
7077 bd = bd*x2+7.18466403235734541950E-11;
7078 y = 1.0/x+x2*bn/(bd*
x);
7087 cn = -5.90592860534773254987E-1;
7088 cn = cn*x2+6.29235242724368800674E-1;
7089 cn = cn*x2-1.72858975380388136411E-1;
7090 cn = cn*x2+1.64837047825189632310E-2;
7091 cn = cn*x2-4.86827613020462700845E-4;
7092 cd = 1.00000000000000000000E0;
7093 cd = cd*x2-2.69820057197544900361E0;
7094 cd = cd*x2+1.73270799045947845857E0;
7095 cd = cd*x2-3.93708582281939493482E-1;
7096 cd = cd*x2+3.44278924041233391079E-2;
7097 cd = cd*x2-9.73655226040941223894E-4;
7098 y = 1.0/x+x2*cn/(cd*
x);
7190 result = 1.3862943611198906188E0-0.5*
ae_log(m1, _state);
7194 p = 1.37982864606273237150E-4;
7195 p = p*m1+2.28025724005875567385E-3;
7196 p = p*m1+7.97404013220415179367E-3;
7197 p = p*m1+9.85821379021226008714E-3;
7198 p = p*m1+6.87489687449949877925E-3;
7199 p = p*m1+6.18901033637687613229E-3;
7200 p = p*m1+8.79078273952743772254E-3;
7201 p = p*m1+1.49380448916805252718E-2;
7202 p = p*m1+3.08851465246711995998E-2;
7203 p = p*m1+9.65735902811690126535E-2;
7204 p = p*m1+1.38629436111989062502E0;
7205 q = 2.94078955048598507511E-5;
7206 q = q*m1+9.14184723865917226571E-4;
7207 q = q*m1+5.94058303753167793257E-3;
7208 q = q*m1+1.54850516649762399335E-2;
7209 q = q*m1+2.39089602715924892727E-2;
7210 q = q*m1+3.01204715227604046988E-2;
7211 q = q*m1+3.73774314173823228969E-2;
7212 q = q*m1+4.88280347570998239232E-2;
7213 q = q*m1+7.03124996963957469739E-2;
7214 q = q*m1+1.24999999999870820058E-1;
7215 q = q*m1+4.99999999999999999821E-1;
7216 result = p-q*
ae_log(m1, _state);
7273 pio2 = 1.57079632679489661923;
7282 result =
ae_log(
ae_tan(0.5*(pio2+phi), _state), _state);
7293 phi = phi-npio2*pio2;
7325 result = temp+npio2*
k;
7337 md =
ae_trunc((phi+pio2)/ae_pi, _state);
7338 t = t*(1.0+temp)/(1.0-temp*t*t);
7350 result = temp+npio2*
k;
7396 p = 1.53552577301013293365E-4;
7397 p = p*m+2.50888492163602060990E-3;
7398 p = p*m+8.68786816565889628429E-3;
7399 p = p*m+1.07350949056076193403E-2;
7400 p = p*m+7.77395492516787092951E-3;
7401 p = p*m+7.58395289413514708519E-3;
7402 p = p*m+1.15688436810574127319E-2;
7403 p = p*m+2.18317996015557253103E-2;
7404 p = p*m+5.68051945617860553470E-2;
7405 p = p*m+4.43147180560990850618E-1;
7406 p = p*m+1.00000000000000000299E0;
7407 q = 3.27954898576485872656E-5;
7408 q = q*m+1.00962792679356715133E-3;
7409 q = q*m+6.50609489976927491433E-3;
7410 q = q*m+1.68862163993311317300E-2;
7411 q = q*m+2.61769742454493659583E-2;
7412 q = q*m+3.34833904888224918614E-2;
7413 q = q*m+4.27180926518931511717E-2;
7414 q = q*m+5.85936634471101055642E-2;
7415 q = q*m+9.37499997197644278445E-2;
7416 q = q*m+2.49999999999888314361E-1;
7417 result = p-q*m*
ae_log(m, _state);
7470 pio2 = 1.57079632679489661923;
7482 lphi = lphi-npio2*pio2;
7496 temp =
ae_sin(lphi, _state);
7501 result = temp+npio2*ebig;
7504 t =
ae_tan(lphi, _state);
7530 result = temp+npio2*ebig;
7543 md =
ae_trunc((lphi+pio2)/ae_pi, _state);
7544 t = t*(1.0+temp)/(1.0-temp*t*t);
7550 e = e+c*
ae_sin(lphi, _state);
7559 result = temp+npio2*ebig;
7601 eul = 0.5772156649015328606065;
7609 f1 = -5.350447357812542947283;
7610 f1 = f1*x+218.5049168816613393830;
7611 f1 = f1*x-4176.572384826693777058;
7612 f1 = f1*x+55411.76756393557601232;
7613 f1 = f1*x-331338.1331178144034309;
7614 f1 = f1*x+1592627.163384945414220;
7615 f2 = 1.000000000000000000000;
7616 f2 = f2*x-52.50547959112862969197;
7617 f2 = f2*x+1259.616186786790571525;
7618 f2 = f2*x-17565.49581973534652631;
7619 f2 = f2*x+149306.2117002725991967;
7620 f2 = f2*x-729494.9239640527645655;
7621 f2 = f2*x+1592627.163384945429726;
7623 result = eul+
ae_log(x, _state)+x*
f;
7629 f1 = 1.981808503259689673238E-2;
7630 f1 = f1*w-1.271645625984917501326;
7631 f1 = f1*w-2.088160335681228318920;
7632 f1 = f1*w+2.755544509187936721172;
7633 f1 = f1*w-4.409507048701600257171E-1;
7634 f1 = f1*w+4.665623805935891391017E-2;
7635 f1 = f1*w-1.545042679673485262580E-3;
7636 f1 = f1*w+7.059980605299617478514E-5;
7637 f2 = 1.000000000000000000000;
7638 f2 = f2*w+1.476498670914921440652;
7639 f2 = f2*w+5.629177174822436244827E-1;
7640 f2 = f2*w+1.699017897879307263248E-1;
7641 f2 = f2*w+2.291647179034212017463E-2;
7642 f2 = f2*w+4.450150439728752875043E-3;
7643 f2 = f2*w+1.727439612206521482874E-4;
7644 f2 = f2*w+3.953167195549672482304E-5;
7646 result =
ae_exp(x, _state)*w*(1+w*
f);
7652 f1 = -1.373215375871208729803;
7653 f1 = f1*w-7.084559133740838761406E-1;
7654 f1 = f1*w+1.580806855547941010501;
7655 f1 = f1*w-2.601500427425622944234E-1;
7656 f1 = f1*w+2.994674694113713763365E-2;
7657 f1 = f1*w-1.038086040188744005513E-3;
7658 f1 = f1*w+4.371064420753005429514E-5;
7659 f1 = f1*w+2.141783679522602903795E-6;
7660 f2 = 1.000000000000000000000;
7661 f2 = f2*w+8.585231423622028380768E-1;
7662 f2 = f2*w+4.483285822873995129957E-1;
7663 f2 = f2*w+7.687932158124475434091E-2;
7664 f2 = f2*w+2.449868241021887685904E-2;
7665 f2 = f2*w+8.832165941927796567926E-4;
7666 f2 = f2*w+4.590952299511353531215E-4;
7667 f2 = f2*w+(-4.729848351866523044863E-6);
7668 f2 = f2*w+2.665195537390710170105E-6;
7670 result =
ae_exp(x, _state)*w*(1+w*
f);
7676 f1 = -2.106934601691916512584;
7677 f1 = f1*w+1.732733869664688041885;
7678 f1 = f1*w-2.423619178935841904839E-1;
7679 f1 = f1*w+2.322724180937565842585E-2;
7680 f1 = f1*w+2.372880440493179832059E-4;
7681 f1 = f1*w-8.343219561192552752335E-5;
7682 f1 = f1*w+1.363408795605250394881E-5;
7683 f1 = f1*w-3.655412321999253963714E-7;
7684 f1 = f1*w+1.464941733975961318456E-8;
7685 f1 = f1*w+6.176407863710360207074E-10;
7686 f2 = 1.000000000000000000000;
7687 f2 = f2*w-2.298062239901678075778E-1;
7688 f2 = f2*w+1.105077041474037862347E-1;
7689 f2 = f2*w-1.566542966630792353556E-2;
7690 f2 = f2*w+2.761106850817352773874E-3;
7691 f2 = f2*w-2.089148012284048449115E-4;
7692 f2 = f2*w+1.708528938807675304186E-5;
7693 f2 = f2*w-4.459311796356686423199E-7;
7694 f2 = f2*w+1.394634930353847498145E-8;
7695 f2 = f2*w+6.150865933977338354138E-10;
7697 result =
ae_exp(x, _state)*w*(1+w*
f);
7703 f1 = -2.458119367674020323359E-1;
7704 f1 = f1*w-1.483382253322077687183E-1;
7705 f1 = f1*w+7.248291795735551591813E-2;
7706 f1 = f1*w-1.348315687380940523823E-2;
7707 f1 = f1*w+1.342775069788636972294E-3;
7708 f1 = f1*w-7.942465637159712264564E-5;
7709 f1 = f1*w+2.644179518984235952241E-6;
7710 f1 = f1*w-4.239473659313765177195E-8;
7711 f2 = 1.000000000000000000000;
7712 f2 = f2*w-1.044225908443871106315E-1;
7713 f2 = f2*w-2.676453128101402655055E-1;
7714 f2 = f2*w+9.695000254621984627876E-2;
7715 f2 = f2*w-1.601745692712991078208E-2;
7716 f2 = f2*w+1.496414899205908021882E-3;
7717 f2 = f2*w-8.462452563778485013756E-5;
7718 f2 = f2*w+2.728938403476726394024E-6;
7719 f2 = f2*w-4.239462431819542051337E-8;
7721 result =
ae_exp(x, _state)*w*(1+w*
f);
7727 f1 = 1.212561118105456670844E-1;
7728 f1 = f1*w-5.823133179043894485122E-1;
7729 f1 = f1*w+2.348887314557016779211E-1;
7730 f1 = f1*w-3.040034318113248237280E-2;
7731 f1 = f1*w+1.510082146865190661777E-3;
7732 f1 = f1*w-2.523137095499571377122E-5;
7733 f2 = 1.000000000000000000000;
7734 f2 = f2*w-1.002252150365854016662;
7735 f2 = f2*w+2.928709694872224144953E-1;
7736 f2 = f2*w-3.337004338674007801307E-2;
7737 f2 = f2*w+1.560544881127388842819E-3;
7738 f2 = f2*w-2.523137093603234562648E-5;
7740 result =
ae_exp(x, _state)*w*(1+w*
f);
7744 f1 = -7.657847078286127362028E-1;
7745 f1 = f1*w+6.886192415566705051750E-1;
7746 f1 = f1*w-2.132598113545206124553E-1;
7747 f1 = f1*w+3.346107552384193813594E-2;
7748 f1 = f1*w-3.076541477344756050249E-3;
7749 f1 = f1*w+1.747119316454907477380E-4;
7750 f1 = f1*w-6.103711682274170530369E-6;
7751 f1 = f1*w+1.218032765428652199087E-7;
7752 f1 = f1*w-1.086076102793290233007E-9;
7753 f2 = 1.000000000000000000000;
7754 f2 = f2*w-1.888802868662308731041;
7755 f2 = f2*w+1.066691687211408896850;
7756 f2 = f2*w-2.751915982306380647738E-1;
7757 f2 = f2*w+3.930852688233823569726E-2;
7758 f2 = f2*w-3.414684558602365085394E-3;
7759 f2 = f2*w+1.866844370703555398195E-4;
7760 f2 = f2*w-6.345146083130515357861E-6;
7761 f2 = f2*w+1.239754287483206878024E-7;
7762 f2 = f2*w-1.086076102793126632978E-9;
7764 result =
ae_exp(x, _state)*w*(1+w*
f);
7821 eul = 0.57721566490153286060;
7822 big = 1.44115188075855872*
ae_pow(10, 17, _state);
7830 result = (double)1/(
double)(n-1);
7835 result =
ae_exp(-x, _state)/
x;
7843 result = yk*t*(6*x*x-8*t*x+t*t);
7844 result = yk*(result+t*(t-2.0*
x));
7845 result = yk*(result+t);
7846 result = (result+1)*
ae_exp(-x, _state)/xk;
7851 psi = -eul-
ae_log(x, _state);
7852 for(i=1; i<=n-1; i++)
7854 psi = psi+(double)1/(
double)
i;
7875 result = result+yk/pk;
7879 t =
ae_fabs(yk/result, _state);
7888 for(i=1; i<=n-1; i++)
7892 result = psi*t-result;
7909 xk = n+(double)(k-1)/(double)2;
7914 xk = (double)k/(
double)2;
7916 pk = pkm1*yk+pkm2*xk;
7917 qk = qkm1*yk+qkm2*xk;
7921 t =
ae_fabs((result-r)/r, _state);
7941 result = result*
ae_exp(-x, _state);
8102 result = (b-b*
w)/(a*w);
8107 result = b*w/(a*(1.0-
w));
8179 mpi = 3.14159265358979323846;
8180 mpio2 = 1.57079632679489661923;
8187 sn = -2.99181919401019853726E3;
8188 sn = sn*t+7.08840045257738576863E5;
8189 sn = sn*t-6.29741486205862506537E7;
8190 sn = sn*t+2.54890880573376359104E9;
8191 sn = sn*t-4.42979518059697779103E10;
8192 sn = sn*t+3.18016297876567817986E11;
8193 sd = 1.00000000000000000000E0;
8194 sd = sd*t+2.81376268889994315696E2;
8195 sd = sd*t+4.55847810806532581675E4;
8196 sd = sd*t+5.17343888770096400730E6;
8197 sd = sd*t+4.19320245898111231129E8;
8198 sd = sd*t+2.24411795645340920940E10;
8199 sd = sd*t+6.07366389490084639049E11;
8200 cn = -4.98843114573573548651E-8;
8201 cn = cn*t+9.50428062829859605134E-6;
8202 cn = cn*t-6.45191435683965050962E-4;
8203 cn = cn*t+1.88843319396703850064E-2;
8204 cn = cn*t-2.05525900955013891793E-1;
8205 cn = cn*t+9.99999999999999998822E-1;
8206 cd = 3.99982968972495980367E-12;
8207 cd = cd*t+9.15439215774657478799E-10;
8208 cd = cd*t+1.25001862479598821474E-7;
8209 cd = cd*t+1.22262789024179030997E-5;
8210 cd = cd*t+8.68029542941784300606E-4;
8211 cd = cd*t+4.12142090722199792936E-2;
8212 cd = cd*t+1.00000000000000000118E0;
8213 *s =
ae_sign(xxa, _state)*x*x2*sn/sd;
8219 *c =
ae_sign(xxa, _state)*0.5;
8220 *s =
ae_sign(xxa, _state)*0.5;
8227 fn = 4.21543555043677546506E-1;
8228 fn = fn*u+1.43407919780758885261E-1;
8229 fn = fn*u+1.15220955073585758835E-2;
8230 fn = fn*u+3.45017939782574027900E-4;
8231 fn = fn*u+4.63613749287867322088E-6;
8232 fn = fn*u+3.05568983790257605827E-8;
8233 fn = fn*u+1.02304514164907233465E-10;
8234 fn = fn*u+1.72010743268161828879E-13;
8235 fn = fn*u+1.34283276233062758925E-16;
8236 fn = fn*u+3.76329711269987889006E-20;
8237 fd = 1.00000000000000000000E0;
8238 fd = fd*u+7.51586398353378947175E-1;
8239 fd = fd*u+1.16888925859191382142E-1;
8240 fd = fd*u+6.44051526508858611005E-3;
8241 fd = fd*u+1.55934409164153020873E-4;
8242 fd = fd*u+1.84627567348930545870E-6;
8243 fd = fd*u+1.12699224763999035261E-8;
8244 fd = fd*u+3.60140029589371370404E-11;
8245 fd = fd*u+5.88754533621578410010E-14;
8246 fd = fd*u+4.52001434074129701496E-17;
8247 fd = fd*u+1.25443237090011264384E-20;
8248 gn = 5.04442073643383265887E-1;
8249 gn = gn*u+1.97102833525523411709E-1;
8250 gn = gn*u+1.87648584092575249293E-2;
8251 gn = gn*u+6.84079380915393090172E-4;
8252 gn = gn*u+1.15138826111884280931E-5;
8253 gn = gn*u+9.82852443688422223854E-8;
8254 gn = gn*u+4.45344415861750144738E-10;
8255 gn = gn*u+1.08268041139020870318E-12;
8256 gn = gn*u+1.37555460633261799868E-15;
8257 gn = gn*u+8.36354435630677421531E-19;
8258 gn = gn*u+1.86958710162783235106E-22;
8259 gd = 1.00000000000000000000E0;
8260 gd = gd*u+1.47495759925128324529E0;
8261 gd = gd*u+3.37748989120019970451E-1;
8262 gd = gd*u+2.53603741420338795122E-2;
8263 gd = gd*u+8.14679107184306179049E-4;
8264 gd = gd*u+1.27545075667729118702E-5;
8265 gd = gd*u+1.04314589657571990585E-7;
8266 gd = gd*u+4.60680728146520428211E-10;
8267 gd = gd*u+1.10273215066240270757E-12;
8268 gd = gd*u+1.38796531259578871258E-15;
8269 gd = gd*u+8.39158816283118707363E-19;
8270 gd = gd*u+1.86958710162783236342E-22;
8277 *c = 0.5+(f*ss-g*cc)/t;
8278 *s = 0.5-(f*cc+g*ss)/t;
8331 result = 2*x*b-2*(i-1)*a;
8399 for(i=0; i<=n/2-1; i++)
8481 ai = 0.25*m*(u-t*
b);
8485 *dn = 1.0-0.5*m*t*t;
8496 *sn = t+ai*(twon-
u)/(b*b);
8497 *ph = 2.0*
ae_atan(
ae_exp(u, _state), _state)-1.57079632679489661923+ai*(twon-
u)/b;
8499 *cn = phi-ai*(twon-
u);
8500 *dn = phi+ai*(twon+
u);
8529 phi = (
ae_asin(t, _state)+phi)/2.0;
8533 *sn =
ae_sin(phi, _state);
8536 *dn = t/
ae_cos(phi-b, _state);
8572 result = ((2*i-1-
x)*b-(i-1)*
a)/i;
8609 result = (2*i+1-
x)*b1/(i+1)-(i+1)*b2/(i+2)+c->
ptr.
p_double[
i];
8636 for(i=0; i<=n-1; i++)
8678 result = ((2*i-1)*x*b-(i-1)*
a)/i;
8715 result = (2*i+1)*x*b1/(i+1)-(i+1)*b2/(i+2)+c->
ptr.
p_double[
i];
8750 for(i=0; i<=n/2-1; i++)
8952 for(i=1; i<=n-1; i++)
8957 y = y-0.57721566490153286061;
8971 polv = 8.33333333333333333333E-2;
8972 polv = polv*z-2.10927960927960927961E-2;
8973 polv = polv*z+7.57575757575757575758E-3;
8974 polv = polv*z-4.16666666666666666667E-3;
8975 polv = polv*z+3.96825396825396825397E-3;
8976 polv = polv*z-8.33333333333333333333E-3;
8977 polv = polv*z+8.33333333333333333333E-2;
8984 y =
ae_log(s, _state)-0.5/s-y-
w;
9051 ae_assert(k>0,
"Domain error in StudentTDistribution", _state);
9085 tz = tz*((j-1)/(z*j));
9100 tz = tz*((j-1)/(z*j));
9104 p = f*x/
ae_sqrt(z*rk, _state);
9152 t =
ae_sqrt(rk*z/(1.0-z), _state);
9259 *si = 1.570796326794896619-
ae_cos(x, _state)/
x;
9266 sn = -8.39167827910303881427E-11;
9267 sn = sn*z+4.62591714427012837309E-8;
9268 sn = sn*z-9.75759303843632795789E-6;
9269 sn = sn*z+9.76945438170435310816E-4;
9270 sn = sn*z-4.13470316229406538752E-2;
9271 sn = sn*z+1.00000000000000000302E0;
9272 sd = 2.03269266195951942049E-12;
9273 sd = sd*z+1.27997891179943299903E-9;
9274 sd = sd*z+4.41827842801218905784E-7;
9275 sd = sd*z+9.96412122043875552487E-5;
9276 sd = sd*z+1.42085239326149893930E-2;
9277 sd = sd*z+9.99999999999999996984E-1;
9279 cn = 2.02524002389102268789E-11;
9280 cn = cn*z-1.35249504915790756375E-8;
9281 cn = cn*z+3.59325051419993077021E-6;
9282 cn = cn*z-4.74007206873407909465E-4;
9283 cn = cn*z+2.89159652607555242092E-2;
9284 cn = cn*z-1.00000000000000000080E0;
9285 cd = 4.07746040061880559506E-12;
9286 cd = cd*z+3.06780997581887812692E-9;
9287 cd = cd*z+1.23210355685883423679E-6;
9288 cd = cd*z+3.17442024775032769882E-4;
9289 cd = cd*z+5.10028056236446052392E-2;
9290 cd = cd*z+4.00000000000000000080E0;
9297 *ci = 0.57721566490153286061+
ae_log(x, _state)+
c;
9305 fn = 4.23612862892216586994E0;
9306 fn = fn*z+5.45937717161812843388E0;
9307 fn = fn*z+1.62083287701538329132E0;
9308 fn = fn*z+1.67006611831323023771E-1;
9309 fn = fn*z+6.81020132472518137426E-3;
9310 fn = fn*z+1.08936580650328664411E-4;
9311 fn = fn*z+5.48900223421373614008E-7;
9312 fd = 1.00000000000000000000E0;
9313 fd = fd*z+8.16496634205391016773E0;
9314 fd = fd*z+7.30828822505564552187E0;
9315 fd = fd*z+1.86792257950184183883E0;
9316 fd = fd*z+1.78792052963149907262E-1;
9317 fd = fd*z+7.01710668322789753610E-3;
9318 fd = fd*z+1.10034357153915731354E-4;
9319 fd = fd*z+5.48900252756255700982E-7;
9321 gn = 8.71001698973114191777E-2;
9322 gn = gn*z+6.11379109952219284151E-1;
9323 gn = gn*z+3.97180296392337498885E-1;
9324 gn = gn*z+7.48527737628469092119E-2;
9325 gn = gn*z+5.38868681462177273157E-3;
9326 gn = gn*z+1.61999794598934024525E-4;
9327 gn = gn*z+1.97963874140963632189E-6;
9328 gn = gn*z+7.82579040744090311069E-9;
9329 gd = 1.00000000000000000000E0;
9330 gd = gd*z+1.64402202413355338886E0;
9331 gd = gd*z+6.66296701268987968381E-1;
9332 gd = gd*z+9.88771761277688796203E-2;
9333 gd = gd*z+6.22396345441768420760E-3;
9334 gd = gd*z+1.73221081474177119497E-4;
9335 gd = gd*z+2.02659182086343991969E-6;
9336 gd = gd*z+7.82579218933534490868E-9;
9341 fn = 4.55880873470465315206E-1;
9342 fn = fn*z+7.13715274100146711374E-1;
9343 fn = fn*z+1.60300158222319456320E-1;
9344 fn = fn*z+1.16064229408124407915E-2;
9345 fn = fn*z+3.49556442447859055605E-4;
9346 fn = fn*z+4.86215430826454749482E-6;
9347 fn = fn*z+3.20092790091004902806E-8;
9348 fn = fn*z+9.41779576128512936592E-11;
9349 fn = fn*z+9.70507110881952024631E-14;
9350 fd = 1.00000000000000000000E0;
9351 fd = fd*z+9.17463611873684053703E-1;
9352 fd = fd*z+1.78685545332074536321E-1;
9353 fd = fd*z+1.22253594771971293032E-2;
9354 fd = fd*z+3.58696481881851580297E-4;
9355 fd = fd*z+4.92435064317881464393E-6;
9356 fd = fd*z+3.21956939101046018377E-8;
9357 fd = fd*z+9.43720590350276732376E-11;
9358 fd = fd*z+9.70507110881952025725E-14;
9360 gn = 6.97359953443276214934E-1;
9361 gn = gn*z+3.30410979305632063225E-1;
9362 gn = gn*z+3.84878767649974295920E-2;
9363 gn = gn*z+1.71718239052347903558E-3;
9364 gn = gn*z+3.48941165502279436777E-5;
9365 gn = gn*z+3.47131167084116673800E-7;
9366 gn = gn*z+1.70404452782044526189E-9;
9367 gn = gn*z+3.85945925430276600453E-12;
9368 gn = gn*z+3.14040098946363334640E-15;
9369 gd = 1.00000000000000000000E0;
9370 gd = gd*z+1.68548898811011640017E0;
9371 gd = gd*z+4.87852258695304967486E-1;
9372 gd = gd*z+4.67913194259625806320E-2;
9373 gd = gd*z+1.90284426674399523638E-3;
9374 gd = gd*z+3.68475504442561108162E-5;
9375 gd = gd*z+3.57043223443740838771E-7;
9376 gd = gd*z+1.72693748966316146736E-9;
9377 gd = gd*z+3.87830166023954706752E-12;
9378 gd = gd*z+3.14040098946363335242E-15;
9381 *si = 1.570796326794896619-f*c-g*s;
9486 a = (576.0/x-52.0)/10.0;
9488 b0 = 1.83889230173399459482E-17;
9490 trigintegrals_chebiterationshichi(a, -9.55485532279655569575E-17, &b0, &b1, &b2, _state);
9491 trigintegrals_chebiterationshichi(a, 2.04326105980879882648E-16, &b0, &b1, &b2, _state);
9492 trigintegrals_chebiterationshichi(a, 1.09896949074905343022E-15, &b0, &b1, &b2, _state);
9493 trigintegrals_chebiterationshichi(a, -1.31313534344092599234E-14, &b0, &b1, &b2, _state);
9494 trigintegrals_chebiterationshichi(a, 5.93976226264314278932E-14, &b0, &b1, &b2, _state);
9495 trigintegrals_chebiterationshichi(a, -3.47197010497749154755E-14, &b0, &b1, &b2, _state);
9496 trigintegrals_chebiterationshichi(a, -1.40059764613117131000E-12, &b0, &b1, &b2, _state);
9497 trigintegrals_chebiterationshichi(a, 9.49044626224223543299E-12, &b0, &b1, &b2, _state);
9498 trigintegrals_chebiterationshichi(a, -1.61596181145435454033E-11, &b0, &b1, &b2, _state);
9499 trigintegrals_chebiterationshichi(a, -1.77899784436430310321E-10, &b0, &b1, &b2, _state);
9500 trigintegrals_chebiterationshichi(a, 1.35455469767246947469E-9, &b0, &b1, &b2, _state);
9501 trigintegrals_chebiterationshichi(a, -1.03257121792819495123E-9, &b0, &b1, &b2, _state);
9502 trigintegrals_chebiterationshichi(a, -3.56699611114982536845E-8, &b0, &b1, &b2, _state);
9503 trigintegrals_chebiterationshichi(a, 1.44818877384267342057E-7, &b0, &b1, &b2, _state);
9504 trigintegrals_chebiterationshichi(a, 7.82018215184051295296E-7, &b0, &b1, &b2, _state);
9505 trigintegrals_chebiterationshichi(a, -5.39919118403805073710E-6, &b0, &b1, &b2, _state);
9506 trigintegrals_chebiterationshichi(a, -3.12458202168959833422E-5, &b0, &b1, &b2, _state);
9507 trigintegrals_chebiterationshichi(a, 8.90136741950727517826E-5, &b0, &b1, &b2, _state);
9508 trigintegrals_chebiterationshichi(a, 2.02558474743846862168E-3, &b0, &b1, &b2, _state);
9509 trigintegrals_chebiterationshichi(a, 2.96064440855633256972E-2, &b0, &b1, &b2, _state);
9510 trigintegrals_chebiterationshichi(a, 1.11847751047257036625E0, &b0, &b1, &b2, _state);
9512 b0 = -8.12435385225864036372E-18;
9514 trigintegrals_chebiterationshichi(a, 2.17586413290339214377E-17, &b0, &b1, &b2, _state);
9515 trigintegrals_chebiterationshichi(a, 5.22624394924072204667E-17, &b0, &b1, &b2, _state);
9516 trigintegrals_chebiterationshichi(a, -9.48812110591690559363E-16, &b0, &b1, &b2, _state);
9517 trigintegrals_chebiterationshichi(a, 5.35546311647465209166E-15, &b0, &b1, &b2, _state);
9518 trigintegrals_chebiterationshichi(a, -1.21009970113732918701E-14, &b0, &b1, &b2, _state);
9519 trigintegrals_chebiterationshichi(a, -6.00865178553447437951E-14, &b0, &b1, &b2, _state);
9520 trigintegrals_chebiterationshichi(a, 7.16339649156028587775E-13, &b0, &b1, &b2, _state);
9521 trigintegrals_chebiterationshichi(a, -2.93496072607599856104E-12, &b0, &b1, &b2, _state);
9522 trigintegrals_chebiterationshichi(a, -1.40359438136491256904E-12, &b0, &b1, &b2, _state);
9523 trigintegrals_chebiterationshichi(a, 8.76302288609054966081E-11, &b0, &b1, &b2, _state);
9524 trigintegrals_chebiterationshichi(a, -4.40092476213282340617E-10, &b0, &b1, &b2, _state);
9525 trigintegrals_chebiterationshichi(a, -1.87992075640569295479E-10, &b0, &b1, &b2, _state);
9526 trigintegrals_chebiterationshichi(a, 1.31458150989474594064E-8, &b0, &b1, &b2, _state);
9527 trigintegrals_chebiterationshichi(a, -4.75513930924765465590E-8, &b0, &b1, &b2, _state);
9528 trigintegrals_chebiterationshichi(a, -2.21775018801848880741E-7, &b0, &b1, &b2, _state);
9529 trigintegrals_chebiterationshichi(a, 1.94635531373272490962E-6, &b0, &b1, &b2, _state);
9530 trigintegrals_chebiterationshichi(a, 4.33505889257316408893E-6, &b0, &b1, &b2, _state);
9531 trigintegrals_chebiterationshichi(a, -6.13387001076494349496E-5, &b0, &b1, &b2, _state);
9532 trigintegrals_chebiterationshichi(a, -3.13085477492997465138E-4, &b0, &b1, &b2, _state);
9533 trigintegrals_chebiterationshichi(a, 4.97164789823116062801E-4, &b0, &b1, &b2, _state);
9534 trigintegrals_chebiterationshichi(a, 2.64347496031374526641E-2, &b0, &b1, &b2, _state);
9535 trigintegrals_chebiterationshichi(a, 1.11446150876699213025E0, &b0, &b1, &b2, _state);
9542 a = (6336.0/x-212.0)/70.0;
9544 b0 = -1.05311574154850938805E-17;
9546 trigintegrals_chebiterationshichi(a, 2.62446095596355225821E-17, &b0, &b1, &b2, _state);
9547 trigintegrals_chebiterationshichi(a, 8.82090135625368160657E-17, &b0, &b1, &b2, _state);
9548 trigintegrals_chebiterationshichi(a, -3.38459811878103047136E-16, &b0, &b1, &b2, _state);
9549 trigintegrals_chebiterationshichi(a, -8.30608026366935789136E-16, &b0, &b1, &b2, _state);
9550 trigintegrals_chebiterationshichi(a, 3.93397875437050071776E-15, &b0, &b1, &b2, _state);
9551 trigintegrals_chebiterationshichi(a, 1.01765565969729044505E-14, &b0, &b1, &b2, _state);
9552 trigintegrals_chebiterationshichi(a, -4.21128170307640802703E-14, &b0, &b1, &b2, _state);
9553 trigintegrals_chebiterationshichi(a, -1.60818204519802480035E-13, &b0, &b1, &b2, _state);
9554 trigintegrals_chebiterationshichi(a, 3.34714954175994481761E-13, &b0, &b1, &b2, _state);
9555 trigintegrals_chebiterationshichi(a, 2.72600352129153073807E-12, &b0, &b1, &b2, _state);
9556 trigintegrals_chebiterationshichi(a, 1.66894954752839083608E-12, &b0, &b1, &b2, _state);
9557 trigintegrals_chebiterationshichi(a, -3.49278141024730899554E-11, &b0, &b1, &b2, _state);
9558 trigintegrals_chebiterationshichi(a, -1.58580661666482709598E-10, &b0, &b1, &b2, _state);
9559 trigintegrals_chebiterationshichi(a, -1.79289437183355633342E-10, &b0, &b1, &b2, _state);
9560 trigintegrals_chebiterationshichi(a, 1.76281629144264523277E-9, &b0, &b1, &b2, _state);
9561 trigintegrals_chebiterationshichi(a, 1.69050228879421288846E-8, &b0, &b1, &b2, _state);
9562 trigintegrals_chebiterationshichi(a, 1.25391771228487041649E-7, &b0, &b1, &b2, _state);
9563 trigintegrals_chebiterationshichi(a, 1.16229947068677338732E-6, &b0, &b1, &b2, _state);
9564 trigintegrals_chebiterationshichi(a, 1.61038260117376323993E-5, &b0, &b1, &b2, _state);
9565 trigintegrals_chebiterationshichi(a, 3.49810375601053973070E-4, &b0, &b1, &b2, _state);
9566 trigintegrals_chebiterationshichi(a, 1.28478065259647610779E-2, &b0, &b1, &b2, _state);
9567 trigintegrals_chebiterationshichi(a, 1.03665722588798326712E0, &b0, &b1, &b2, _state);
9569 b0 = 8.06913408255155572081E-18;
9571 trigintegrals_chebiterationshichi(a, -2.08074168180148170312E-17, &b0, &b1, &b2, _state);
9572 trigintegrals_chebiterationshichi(a, -5.98111329658272336816E-17, &b0, &b1, &b2, _state);
9573 trigintegrals_chebiterationshichi(a, 2.68533951085945765591E-16, &b0, &b1, &b2, _state);
9574 trigintegrals_chebiterationshichi(a, 4.52313941698904694774E-16, &b0, &b1, &b2, _state);
9575 trigintegrals_chebiterationshichi(a, -3.10734917335299464535E-15, &b0, &b1, &b2, _state);
9576 trigintegrals_chebiterationshichi(a, -4.42823207332531972288E-15, &b0, &b1, &b2, _state);
9577 trigintegrals_chebiterationshichi(a, 3.49639695410806959872E-14, &b0, &b1, &b2, _state);
9578 trigintegrals_chebiterationshichi(a, 6.63406731718911586609E-14, &b0, &b1, &b2, _state);
9579 trigintegrals_chebiterationshichi(a, -3.71902448093119218395E-13, &b0, &b1, &b2, _state);
9580 trigintegrals_chebiterationshichi(a, -1.27135418132338309016E-12, &b0, &b1, &b2, _state);
9581 trigintegrals_chebiterationshichi(a, 2.74851141935315395333E-12, &b0, &b1, &b2, _state);
9582 trigintegrals_chebiterationshichi(a, 2.33781843985453438400E-11, &b0, &b1, &b2, _state);
9583 trigintegrals_chebiterationshichi(a, 2.71436006377612442764E-11, &b0, &b1, &b2, _state);
9584 trigintegrals_chebiterationshichi(a, -2.56600180000355990529E-10, &b0, &b1, &b2, _state);
9585 trigintegrals_chebiterationshichi(a, -1.61021375163803438552E-9, &b0, &b1, &b2, _state);
9586 trigintegrals_chebiterationshichi(a, -4.72543064876271773512E-9, &b0, &b1, &b2, _state);
9587 trigintegrals_chebiterationshichi(a, -3.00095178028681682282E-9, &b0, &b1, &b2, _state);
9588 trigintegrals_chebiterationshichi(a, 7.79387474390914922337E-8, &b0, &b1, &b2, _state);
9589 trigintegrals_chebiterationshichi(a, 1.06942765566401507066E-6, &b0, &b1, &b2, _state);
9590 trigintegrals_chebiterationshichi(a, 1.59503164802313196374E-5, &b0, &b1, &b2, _state);
9591 trigintegrals_chebiterationshichi(a, 3.49592575153777996871E-4, &b0, &b1, &b2, _state);
9592 trigintegrals_chebiterationshichi(a, 1.28475387530065247392E-2, &b0, &b1, &b2, _state);
9593 trigintegrals_chebiterationshichi(a, 1.03665693917934275131E0, &b0, &b1, &b2, _state);
9616 *chi = 0.57721566490153286061+
ae_log(x, _state)+
c;
9620 static void trigintegrals_chebiterationshichi(
double x,
9631 *b0 = x*(*b1)-(*b2)+
c;
struct alglib_impl::ae_state ae_state
void hyperbolicsinecosineintegrals(double x, double *shi, double *chi, ae_state *_state)
double beta(double a, double b, ae_state *_state)
ae_bool ae_fp_greater_eq(double v1, double v2)
double incompleteellipticintegralk(double phi, double m, ae_state *_state)
double invbinomialdistribution(ae_int_t k, ae_int_t n, double y, ae_state *_state)
double poissoncdistribution(ae_int_t k, double m, ae_state *_state)
void airy(double x, double *ai, double *aip, double *bi, double *bip, ae_state *_state)
double invfdistribution(ae_int_t a, ae_int_t b, double y, ae_state *_state)
double besselj1(double x, ae_state *_state)
double binomialcdistribution(ae_int_t k, ae_int_t n, double p, ae_state *_state)
double ae_sin(double x, ae_state *state)
double bessely1(double x, ae_state *_state)
double ae_fabs(double x, ae_state *state)
double errorfunction(double x, ae_state *_state)
double ae_tan(double x, ae_state *state)
double chebyshevsum(ae_vector *c, ae_int_t r, ae_int_t n, double x, ae_state *_state)
double ae_pow(double x, double y, ae_state *state)
void ae_frame_make(ae_state *state, ae_frame *tmp)
void fresnelintegral(double x, double *c, double *s, ae_state *_state)
void fromchebyshev(ae_vector *a, ae_int_t n, ae_vector *b, ae_state *_state)
void jacobianellipticfunctions(double u, double m, double *sn, double *cn, double *dn, double *ph, ae_state *_state)
double psi(double x, ae_state *_state)
void hermitecoefficients(ae_int_t n, ae_vector *c, ae_state *_state)
double invincompletegammac(double a, double y0, ae_state *_state)
double besseli0(double x, ae_state *_state)
ae_int_t ae_sign(double x, ae_state *state)
double normaldistribution(double x, ae_state *_state)
void ae_state_clear(ae_state *state)
double hermitecalculate(ae_int_t n, double x, ae_state *_state)
ae_bool ae_fp_eq(double v1, double v2)
double exponentialintegralen(double x, ae_int_t n, ae_state *_state)
double fcdistribution(ae_int_t a, ae_int_t b, double x, ae_state *_state)
double incompletebeta(double a, double b, double x, ae_state *_state)
double ae_cos(double x, ae_state *state)
ql0001_ & k(htemp+1),(cvec+1),(atemp+1),(bj+1),(bl+1),(bu+1),(x+1),(clamda+1), &iout, infoqp, &zero,(w+1), &lenw,(iw+1), &leniw, &glob_grd.epsmac
double studenttdistribution(ae_int_t k, double t, ae_state *_state)
void laguerrecoefficients(ae_int_t n, ae_vector *c, ae_state *_state)
double ae_asin(double x, ae_state *state)
double chisquaredistribution(double v, double x, ae_state *_state)
double chisquarecdistribution(double v, double x, ae_state *_state)
double fdistribution(ae_int_t a, ae_int_t b, double x, ae_state *_state)
double incompletegammac(double a, double x, ae_state *_state)
double laguerrecalculate(ae_int_t n, double x, ae_state *_state)
void ae_vector_clear(ae_vector *dst)
double invchisquaredistribution(double v, double y, ae_state *_state)
double incompleteellipticintegrale(double phi, double m, ae_state *_state)
double errorfunctionc(double x, ae_state *_state)
ae_bool ae_fp_less(double v1, double v2)
double exponentialintegralei(double x, ae_state *_state)
double invpoissondistribution(ae_int_t k, double y, ae_state *_state)
ae_bool ae_fp_neq(double v1, double v2)
double dawsonintegral(double x, ae_state *_state)
void chebyshevcoefficients(ae_int_t n, ae_vector *c, ae_state *_state)
double ellipticintegralkhighprecision(double m1, ae_state *_state)
double besselk0(double x, ae_state *_state)
double bessely0(double x, ae_state *_state)
ae_bool ae_vector_set_length(ae_vector *dst, ae_int_t newsize, ae_state *state)
double chebyshevcalculate(ae_int_t r, ae_int_t n, double x, ae_state *_state)
double besselj0(double x, ae_state *_state)
double ae_log(double x, ae_state *state)
double gammafunction(double x, ae_state *_state)
double invincompletebeta(double a, double b, double y, ae_state *_state)
double nuexpm1(double x, ae_state *_state)
const alglib_impl::ae_vector * c_ptr() const
double invnormaldistribution(double y0, ae_state *_state)
double ellipticintegralk(double m, ae_state *_state)
double invstudenttdistribution(ae_int_t k, double p, ae_state *_state)
double ae_sinh(double x, ae_state *state)
ae_int_t ae_ifloor(double x, ae_state *state)
void ae_state_init(ae_state *state)
double ae_sqrt(double x, ae_state *state)
void ae_assert(ae_bool cond, const char *msg, ae_state *state)
union alglib_impl::ae_vector::@11 ptr
const char *volatile error_msg
double hermitesum(ae_vector *c, ae_int_t n, double x, ae_state *_state)
double ae_atan(double x, ae_state *state)
double besselk1(double x, ae_state *_state)
double nulog1p(double x, ae_state *_state)
#define ae_machineepsilon
double ae_cosh(double x, ae_state *state)
double ae_exp(double x, ae_state *state)
double ae_tanh(double x, ae_state *state)
void legendrecoefficients(ae_int_t n, ae_vector *c, ae_state *_state)
double besselyn(ae_int_t n, double x, ae_state *_state)
double besseljn(ae_int_t n, double x, ae_state *_state)
ae_bool ae_vector_init(ae_vector *dst, ae_int_t size, ae_datatype datatype, ae_state *state, ae_bool make_automatic)
double ae_sqr(double x, ae_state *state)
void sinecosineintegrals(double x, double *si, double *ci, ae_state *_state)
double incompletegamma(double a, double x, ae_state *_state)
double binomialdistribution(ae_int_t k, ae_int_t n, double p, ae_state *_state)
double ellipticintegrale(double m, ae_state *_state)
double inverf(double e, ae_state *_state)
ae_bool ae_fp_less_eq(double v1, double v2)
ae_int_t ae_round(double x, ae_state *state)
double lngamma(double x, double *sgngam, ae_state *_state)
alglib_impl::ae_int_t ae_int_t
void ae_frame_leave(ae_state *state)
ae_bool ae_fp_greater(double v1, double v2)
double besselkn(ae_int_t nn, double x, ae_state *_state)
double legendresum(ae_vector *c, ae_int_t n, double x, ae_state *_state)
double besseli1(double x, ae_state *_state)
ae_int_t ae_trunc(double x, ae_state *state)
double laguerresum(ae_vector *c, ae_int_t n, double x, ae_state *_state)
double legendrecalculate(ae_int_t n, double x, ae_state *_state)
double poissondistribution(ae_int_t k, double m, ae_state *_state)