97 1.60119522476751861407E-4,
98 1.19135147006586384913E-3,
99 1.04213797561761569935E-2,
100 4.76367800457137231464E-2,
101 2.07448227648435975150E-1,
102 4.94214826801497100753E-1,
103 9.99999999999999996796E-1
106 static double Q[] = {
107 -2.31581873324120129819E-5,
108 5.39605580493303397842E-4,
109 -4.45641913851797240494E-3,
110 1.18139785222060435552E-2,
111 3.58236398605498653373E-2,
112 -2.34591795718243348568E-1,
113 7.14304917030273074085E-2,
114 1.00000000000000000320E0
117 #define MAXGAM 171.624376956302725
118 static double LOGPI = 1.14472988584940017414;
122 7.87311395793093628397E-4,
123 -2.29549961613378126380E-4,
124 -2.68132617805781232825E-3,
125 3.47222221605458667310E-3,
126 8.33333333333482257126E-2,
129 #define MAXSTIR 143.01608
130 static double SQTPI = 2.50662827463100050242E0;
133 static double stirf(
double);
149 v =
pow(
x, 0.5 *
x - 0.25);
189 return (sgngam * INFINITY);
233 return (
z / ((1.0 + 0.5772156649015329 *
x) *
x));
241 static double A[] = {
242 8.11614167470508450300E-4,
243 -5.95061904284301438324E-4,
244 7.93650340457716943945E-4,
245 -2.77777777730099687205E-3,
246 8.33333333333331927722E-2
249 static double B[] = {
250 -1.37825152569120859100E3,
251 -3.88016315134637840924E4,
252 -3.31612992738871184744E5,
253 -1.16237097492762307383E6,
254 -1.72173700820839662146E6,
255 -8.53555664245765465627E5
258 static double C[] = {
260 -3.51815701436523470549E2,
261 -1.70642106651881159223E4,
262 -2.20528590553854454839E5,
263 -1.13933444367982507207E6,
264 -2.53252307177582951285E6,
265 -2.01889141433532773231E6
269 static double LS2PI = 0.91893853320467274178;
271 #define MAXLGM 2.556348e305
283 double p,
q, u,
w,
z;
349 return (*
sign * INFINITY);
358 q += ((7.9365079365079365079365e-4 *
p
359 - 2.7777777777777777777778e-3) *
p
360 + 0.0833333333333333333333) /
x;