\frac{\left(\left(\left(\left(1 + 0.1049934947 \cdot \left(x \cdot x\right)\right) + 0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0005064034 \cdot \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0001789971 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(\left(\left(\left(\left(1 + 0.7715471019 \cdot \left(x \cdot x\right)\right) + 0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0694555761 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0140005442 \cdot \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.0008327945 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + \left(2 \cdot 0.0001789971\right) \cdot \left(\left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)} \cdot x\begin{array}{l}
\mathbf{if}\;x \le -4578.839371396249:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\mathbf{elif}\;x \le 1189.8967229595753:\\
\;\;\;\;x \cdot \sqrt[3]{\frac{\left(\left(0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.1049934947 \cdot \left(x \cdot x\right) + 1\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0005064034\right) + 0.0001789971 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(\left(x \cdot x\right) \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.0001789971 \cdot 2\right) + \left(0.0008327945 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)} \cdot \left(\left(\left(\sqrt[3]{\frac{\left(\left(0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.1049934947 \cdot \left(x \cdot x\right) + 1\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0005064034\right) + 0.0001789971 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(\left(x \cdot x\right) \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.0001789971 \cdot 2\right) + \left(0.0008327945 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}} \cdot \sqrt[3]{\frac{\left(\left(0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.1049934947 \cdot \left(x \cdot x\right) + 1\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0005064034\right) + 0.0001789971 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(\left(x \cdot x\right) \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.0001789971 \cdot 2\right) + \left(0.0008327945 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}}\right) \cdot \sqrt[3]{\frac{\left(\left(0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.1049934947 \cdot \left(x \cdot x\right) + 1\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0005064034\right) + 0.0001789971 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(\left(x \cdot x\right) \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.0001789971 \cdot 2\right) + \left(0.0008327945 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}}\right) \cdot \frac{\left(\left(0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.1049934947 \cdot \left(x \cdot x\right) + 1\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0005064034\right) + 0.0001789971 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(\left(x \cdot x\right) \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.0001789971 \cdot 2\right) + \left(0.0008327945 \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\end{array}double f(double x) {
double r40861503 = 1.0;
double r40861504 = 0.1049934947;
double r40861505 = x;
double r40861506 = r40861505 * r40861505;
double r40861507 = r40861504 * r40861506;
double r40861508 = r40861503 + r40861507;
double r40861509 = 0.0424060604;
double r40861510 = r40861506 * r40861506;
double r40861511 = r40861509 * r40861510;
double r40861512 = r40861508 + r40861511;
double r40861513 = 0.0072644182;
double r40861514 = r40861510 * r40861506;
double r40861515 = r40861513 * r40861514;
double r40861516 = r40861512 + r40861515;
double r40861517 = 0.0005064034;
double r40861518 = r40861514 * r40861506;
double r40861519 = r40861517 * r40861518;
double r40861520 = r40861516 + r40861519;
double r40861521 = 0.0001789971;
double r40861522 = r40861518 * r40861506;
double r40861523 = r40861521 * r40861522;
double r40861524 = r40861520 + r40861523;
double r40861525 = 0.7715471019;
double r40861526 = r40861525 * r40861506;
double r40861527 = r40861503 + r40861526;
double r40861528 = 0.2909738639;
double r40861529 = r40861528 * r40861510;
double r40861530 = r40861527 + r40861529;
double r40861531 = 0.0694555761;
double r40861532 = r40861531 * r40861514;
double r40861533 = r40861530 + r40861532;
double r40861534 = 0.0140005442;
double r40861535 = r40861534 * r40861518;
double r40861536 = r40861533 + r40861535;
double r40861537 = 0.0008327945;
double r40861538 = r40861537 * r40861522;
double r40861539 = r40861536 + r40861538;
double r40861540 = 2.0;
double r40861541 = r40861540 * r40861521;
double r40861542 = r40861522 * r40861506;
double r40861543 = r40861541 * r40861542;
double r40861544 = r40861539 + r40861543;
double r40861545 = r40861524 / r40861544;
double r40861546 = r40861545 * r40861505;
return r40861546;
}
double f(double x) {
double r40861547 = x;
double r40861548 = -4578.839371396249;
bool r40861549 = r40861547 <= r40861548;
double r40861550 = 0.5;
double r40861551 = 0.2514179000665375;
double r40861552 = r40861547 * r40861547;
double r40861553 = r40861551 / r40861552;
double r40861554 = r40861550 + r40861553;
double r40861555 = r40861554 / r40861547;
double r40861556 = 0.15298196345929327;
double r40861557 = 5.0;
double r40861558 = pow(r40861547, r40861557);
double r40861559 = r40861556 / r40861558;
double r40861560 = r40861555 + r40861559;
double r40861561 = 1189.8967229595753;
bool r40861562 = r40861547 <= r40861561;
double r40861563 = 0.0072644182;
double r40861564 = r40861552 * r40861552;
double r40861565 = r40861564 * r40861552;
double r40861566 = r40861563 * r40861565;
double r40861567 = 0.0424060604;
double r40861568 = r40861567 * r40861564;
double r40861569 = 0.1049934947;
double r40861570 = r40861569 * r40861552;
double r40861571 = 1.0;
double r40861572 = r40861570 + r40861571;
double r40861573 = r40861568 + r40861572;
double r40861574 = r40861566 + r40861573;
double r40861575 = r40861565 * r40861552;
double r40861576 = 0.0005064034;
double r40861577 = r40861575 * r40861576;
double r40861578 = r40861574 + r40861577;
double r40861579 = 0.0001789971;
double r40861580 = r40861575 * r40861552;
double r40861581 = r40861579 * r40861580;
double r40861582 = r40861578 + r40861581;
double r40861583 = r40861552 * r40861580;
double r40861584 = 2.0;
double r40861585 = r40861579 * r40861584;
double r40861586 = r40861583 * r40861585;
double r40861587 = 0.0008327945;
double r40861588 = r40861587 * r40861580;
double r40861589 = 0.0140005442;
double r40861590 = r40861575 * r40861589;
double r40861591 = 0.0694555761;
double r40861592 = r40861565 * r40861591;
double r40861593 = 0.2909738639;
double r40861594 = r40861593 * r40861564;
double r40861595 = 0.7715471019;
double r40861596 = r40861595 * r40861552;
double r40861597 = r40861596 + r40861571;
double r40861598 = r40861594 + r40861597;
double r40861599 = r40861592 + r40861598;
double r40861600 = r40861590 + r40861599;
double r40861601 = r40861588 + r40861600;
double r40861602 = r40861586 + r40861601;
double r40861603 = r40861582 / r40861602;
double r40861604 = cbrt(r40861603);
double r40861605 = r40861604 * r40861604;
double r40861606 = r40861605 * r40861604;
double r40861607 = r40861606 * r40861603;
double r40861608 = r40861603 * r40861607;
double r40861609 = cbrt(r40861608);
double r40861610 = r40861547 * r40861609;
double r40861611 = r40861562 ? r40861610 : r40861560;
double r40861612 = r40861549 ? r40861560 : r40861611;
return r40861612;
}



Bits error versus x
Results
if x < -4578.839371396249 or 1189.8967229595753 < x Initial program 58.3
Taylor expanded around -inf 0.0
Simplified0.0
if -4578.839371396249 < x < 1189.8967229595753Initial program 0.0
rmApplied add-cbrt-cube0.0
rmApplied add-cube-cbrt0.0
Final simplification0.0
herbie shell --seed 2019124
(FPCore (x)
:name "Jmat.Real.dawson"
(* (/ (+ (+ (+ (+ (+ 1 (* 0.1049934947 (* x x))) (* 0.0424060604 (* (* x x) (* x x)))) (* 0.0072644182 (* (* (* x x) (* x x)) (* x x)))) (* 0.0005064034 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (+ (+ (+ (+ (+ (+ 1 (* 0.7715471019 (* x x))) (* 0.2909738639 (* (* x x) (* x x)))) (* 0.0694555761 (* (* (* x x) (* x x)) (* x x)))) (* 0.0140005442 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0008327945 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (* 2 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x))