27 static const double cof[6] = {76.18009172947146,-86.50532032941677,
28 24.01409824083091, -1.231739572450155,
29 0.1208650973866179e-2, -0.5395239384953e-5};
33 tmp -= (x + 0.5) * std::log(tmp);
34 double ser = 1.000000000190015;
36 for ( j = 0; j <= 5; j++ ) {
40 return -tmp + std::log(2.5066282746310005*ser);