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.9999999999998099 + \frac{676.5203681218851}{\left(z - 1\right) + 1}\right) + \frac{-1259.1392167224028}{\left(z - 1\right) + 2}\right) + \frac{771.3234287776531}{\left(z - 1\right) + 3}\right) + \frac{-176.6150291621406}{\left(z - 1\right) + 4}\right) + \frac{12.507343278686905}{\left(z - 1\right) + 5}\right) + \frac{-0.13857109526572012}{\left(z - 1\right) + 6}\right) + \frac{9.984369578019572 \cdot 10^{-06}}{\left(z - 1\right) + 7}\right) + \frac{1.5056327351493116 \cdot 10^{-07}}{\left(z - 1\right) + 8}\right)double f(double z) {
double r74470561 = atan2(1.0, 0.0);
double r74470562 = 2.0;
double r74470563 = r74470561 * r74470562;
double r74470564 = sqrt(r74470563);
double r74470565 = z;
double r74470566 = 1.0;
double r74470567 = r74470565 - r74470566;
double r74470568 = 7.0;
double r74470569 = r74470567 + r74470568;
double r74470570 = 0.5;
double r74470571 = r74470569 + r74470570;
double r74470572 = r74470567 + r74470570;
double r74470573 = pow(r74470571, r74470572);
double r74470574 = r74470564 * r74470573;
double r74470575 = -r74470571;
double r74470576 = exp(r74470575);
double r74470577 = r74470574 * r74470576;
double r74470578 = 0.9999999999998099;
double r74470579 = 676.5203681218851;
double r74470580 = r74470567 + r74470566;
double r74470581 = r74470579 / r74470580;
double r74470582 = r74470578 + r74470581;
double r74470583 = -1259.1392167224028;
double r74470584 = r74470567 + r74470562;
double r74470585 = r74470583 / r74470584;
double r74470586 = r74470582 + r74470585;
double r74470587 = 771.3234287776531;
double r74470588 = 3.0;
double r74470589 = r74470567 + r74470588;
double r74470590 = r74470587 / r74470589;
double r74470591 = r74470586 + r74470590;
double r74470592 = -176.6150291621406;
double r74470593 = 4.0;
double r74470594 = r74470567 + r74470593;
double r74470595 = r74470592 / r74470594;
double r74470596 = r74470591 + r74470595;
double r74470597 = 12.507343278686905;
double r74470598 = 5.0;
double r74470599 = r74470567 + r74470598;
double r74470600 = r74470597 / r74470599;
double r74470601 = r74470596 + r74470600;
double r74470602 = -0.13857109526572012;
double r74470603 = 6.0;
double r74470604 = r74470567 + r74470603;
double r74470605 = r74470602 / r74470604;
double r74470606 = r74470601 + r74470605;
double r74470607 = 9.984369578019572e-06;
double r74470608 = r74470607 / r74470569;
double r74470609 = r74470606 + r74470608;
double r74470610 = 1.5056327351493116e-07;
double r74470611 = 8.0;
double r74470612 = r74470567 + r74470611;
double r74470613 = r74470610 / r74470612;
double r74470614 = r74470609 + r74470613;
double r74470615 = r74470577 * r74470614;
return r74470615;
}
herbie shell --seed 2019125 +o rules:numerics
(FPCore (z)
:name "Jmat.Real.gamma, branch z greater than 0.5"
(* (* (* (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)))))