\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\begin{array}{l}
\mathbf{if}\;\left(t - 1\right) \cdot \log a \le -5.03697642486000799 \cdot 10^{32} \lor \neg \left(\left(t - 1\right) \cdot \log a \le 65.394405318574371\right):\\
\;\;\;\;\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{{1}^{1}}}{\frac{\frac{e^{b - y \cdot \log z}}{{a}^{t}}}{\frac{\frac{1}{y}}{{a}^{1}}}}\\
\end{array}double f(double x, double y, double z, double t, double a, double b) {
double r100468 = x;
double r100469 = y;
double r100470 = z;
double r100471 = log(r100470);
double r100472 = r100469 * r100471;
double r100473 = t;
double r100474 = 1.0;
double r100475 = r100473 - r100474;
double r100476 = a;
double r100477 = log(r100476);
double r100478 = r100475 * r100477;
double r100479 = r100472 + r100478;
double r100480 = b;
double r100481 = r100479 - r100480;
double r100482 = exp(r100481);
double r100483 = r100468 * r100482;
double r100484 = r100483 / r100469;
return r100484;
}
double f(double x, double y, double z, double t, double a, double b) {
double r100485 = t;
double r100486 = 1.0;
double r100487 = r100485 - r100486;
double r100488 = a;
double r100489 = log(r100488);
double r100490 = r100487 * r100489;
double r100491 = -5.036976424860008e+32;
bool r100492 = r100490 <= r100491;
double r100493 = 65.39440531857437;
bool r100494 = r100490 <= r100493;
double r100495 = !r100494;
bool r100496 = r100492 || r100495;
double r100497 = x;
double r100498 = y;
double r100499 = z;
double r100500 = log(r100499);
double r100501 = r100498 * r100500;
double r100502 = r100501 + r100490;
double r100503 = b;
double r100504 = r100502 - r100503;
double r100505 = exp(r100504);
double r100506 = r100497 * r100505;
double r100507 = r100506 / r100498;
double r100508 = 1.0;
double r100509 = pow(r100508, r100486);
double r100510 = r100497 / r100509;
double r100511 = r100503 - r100501;
double r100512 = exp(r100511);
double r100513 = pow(r100488, r100485);
double r100514 = r100512 / r100513;
double r100515 = r100508 / r100498;
double r100516 = pow(r100488, r100486);
double r100517 = r100515 / r100516;
double r100518 = r100514 / r100517;
double r100519 = r100510 / r100518;
double r100520 = r100496 ? r100507 : r100519;
return r100520;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b
Results
if (* (- t 1.0) (log a)) < -5.036976424860008e+32 or 65.39440531857437 < (* (- t 1.0) (log a)) Initial program 0.5
if -5.036976424860008e+32 < (* (- t 1.0) (log a)) < 65.39440531857437Initial program 4.7
Taylor expanded around inf 4.7
Simplified8.5
rmApplied *-un-lft-identity8.5
Applied unpow-prod-down8.5
Applied div-inv8.5
Applied times-frac2.6
Applied associate-/l*0.9
Simplified0.9
Final simplification0.6
herbie shell --seed 2020046
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
:precision binary64
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1) (log a))) b))) y))