Use the --timeout flag to change the timeout.
\left(\left(\sqrt{\pi \cdot 2} \cdot {\left(\left(\left(z - 1\right) + 7\right) + 0.5\right)}^{\left(\left(z - 1\right) + 0.5\right)}\right) \cdot e^{-\left(\left(\left(z - 1\right) + 7\right) + 0.5\right)}\right) \cdot \left(\left(\left(\left(\left(\left(\left(\left(0.99999999999980993 + \frac{676.520368121885099}{\left(z - 1\right) + 1}\right) + \frac{-1259.13921672240281}{\left(z - 1\right) + 2}\right) + \frac{771.32342877765313}{\left(z - 1\right) + 3}\right) + \frac{-176.615029162140587}{\left(z - 1\right) + 4}\right) + \frac{12.5073432786869052}{\left(z - 1\right) + 5}\right) + \frac{-0.138571095265720118}{\left(z - 1\right) + 6}\right) + \frac{9.98436957801957158 \cdot 10^{-6}}{\left(z - 1\right) + 7}\right) + \frac{1.50563273514931162 \cdot 10^{-7}}{\left(z - 1\right) + 8}\right)double f(double z) {
double r166615 = atan2(1.0, 0.0);
double r166616 = 2.0;
double r166617 = r166615 * r166616;
double r166618 = sqrt(r166617);
double r166619 = z;
double r166620 = 1.0;
double r166621 = r166619 - r166620;
double r166622 = 7.0;
double r166623 = r166621 + r166622;
double r166624 = 0.5;
double r166625 = r166623 + r166624;
double r166626 = r166621 + r166624;
double r166627 = pow(r166625, r166626);
double r166628 = r166618 * r166627;
double r166629 = -r166625;
double r166630 = exp(r166629);
double r166631 = r166628 * r166630;
double r166632 = 0.9999999999998099;
double r166633 = 676.5203681218851;
double r166634 = r166621 + r166620;
double r166635 = r166633 / r166634;
double r166636 = r166632 + r166635;
double r166637 = -1259.1392167224028;
double r166638 = r166621 + r166616;
double r166639 = r166637 / r166638;
double r166640 = r166636 + r166639;
double r166641 = 771.3234287776531;
double r166642 = 3.0;
double r166643 = r166621 + r166642;
double r166644 = r166641 / r166643;
double r166645 = r166640 + r166644;
double r166646 = -176.6150291621406;
double r166647 = 4.0;
double r166648 = r166621 + r166647;
double r166649 = r166646 / r166648;
double r166650 = r166645 + r166649;
double r166651 = 12.507343278686905;
double r166652 = 5.0;
double r166653 = r166621 + r166652;
double r166654 = r166651 / r166653;
double r166655 = r166650 + r166654;
double r166656 = -0.13857109526572012;
double r166657 = 6.0;
double r166658 = r166621 + r166657;
double r166659 = r166656 / r166658;
double r166660 = r166655 + r166659;
double r166661 = 9.984369578019572e-06;
double r166662 = r166661 / r166623;
double r166663 = r166660 + r166662;
double r166664 = 1.5056327351493116e-07;
double r166665 = 8.0;
double r166666 = r166621 + r166665;
double r166667 = r166664 / r166666;
double r166668 = r166663 + r166667;
double r166669 = r166631 * r166668;
return r166669;
}
herbie shell --seed 2020083 +o rules:numerics
(FPCore (z)
:name "Jmat.Real.gamma, branch z greater than 0.5"
:precision binary64
(* (* (* (sqrt (* PI 2)) (pow (+ (+ (- z 1) 7) 0.5) (+ (- z 1) 0.5))) (exp (- (+ (+ (- z 1) 7) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (+ (- z 1) 1))) (/ -1259.1392167224028 (+ (- z 1) 2))) (/ 771.3234287776531 (+ (- z 1) 3))) (/ -176.6150291621406 (+ (- z 1) 4))) (/ 12.507343278686905 (+ (- z 1) 5))) (/ -0.13857109526572012 (+ (- z 1) 6))) (/ 9.984369578019572e-06 (+ (- z 1) 7))) (/ 1.5056327351493116e-07 (+ (- z 1) 8)))))