\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 -621.4550718514375:\\
\;\;\;\;\frac{\frac{0.2514179000665375}{x}}{x \cdot x} + \left(\frac{0.5}{x} + \frac{\frac{0.15298196345929327}{x \cdot x}}{\left(x \cdot x\right) \cdot x}\right)\\
\mathbf{elif}\;x \le 640.4757042222545:\\
\;\;\;\;x \cdot \sqrt[3]{\frac{\mathsf{fma}\left(\left(x \cdot x\right) \cdot 0.0001789971, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(0.0005064034, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0072644182, \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right), 0.0424060604, \mathsf{fma}\left(0.1049934947, x \cdot x, 1\right)\right)\right)\right)\right)}{\mathsf{fma}\left(2 \cdot 0.0001789971, \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0008327945 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0140005442 \cdot \left(x \cdot x\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0694555761, \mathsf{fma}\left(0.2909738639, \left(x \cdot x\right) \cdot \left(x \cdot x\right), \mathsf{fma}\left(x \cdot x, 0.7715471019, 1\right)\right)\right)\right)\right)\right)} \cdot \left(\frac{\mathsf{fma}\left(\left(x \cdot x\right) \cdot 0.0001789971, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(0.0005064034, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0072644182, \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right), 0.0424060604, \mathsf{fma}\left(0.1049934947, x \cdot x, 1\right)\right)\right)\right)\right)}{\mathsf{fma}\left(2 \cdot 0.0001789971, \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0008327945 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0140005442 \cdot \left(x \cdot x\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0694555761, \mathsf{fma}\left(0.2909738639, \left(x \cdot x\right) \cdot \left(x \cdot x\right), \mathsf{fma}\left(x \cdot x, 0.7715471019, 1\right)\right)\right)\right)\right)\right)} \cdot \frac{\mathsf{fma}\left(\left(x \cdot x\right) \cdot 0.0001789971, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(0.0005064034, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(x \cdot x, \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0072644182, \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right), 0.0424060604, \mathsf{fma}\left(0.1049934947, x \cdot x, 1\right)\right)\right)\right)\right)}{\mathsf{fma}\left(2 \cdot 0.0001789971, \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0008327945 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0140005442 \cdot \left(x \cdot x\right), \mathsf{fma}\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right), 0.0694555761, \mathsf{fma}\left(0.2909738639, \left(x \cdot x\right) \cdot \left(x \cdot x\right), \mathsf{fma}\left(x \cdot x, 0.7715471019, 1\right)\right)\right)\right)\right)\right)}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{0.2514179000665375}{x}}{x \cdot x} + \left(\frac{0.5}{x} + \frac{\frac{0.15298196345929327}{x \cdot x}}{\left(x \cdot x\right) \cdot x}\right)\\
\end{array}double f(double x) {
double r4819507 = 1.0;
double r4819508 = 0.1049934947;
double r4819509 = x;
double r4819510 = r4819509 * r4819509;
double r4819511 = r4819508 * r4819510;
double r4819512 = r4819507 + r4819511;
double r4819513 = 0.0424060604;
double r4819514 = r4819510 * r4819510;
double r4819515 = r4819513 * r4819514;
double r4819516 = r4819512 + r4819515;
double r4819517 = 0.0072644182;
double r4819518 = r4819514 * r4819510;
double r4819519 = r4819517 * r4819518;
double r4819520 = r4819516 + r4819519;
double r4819521 = 0.0005064034;
double r4819522 = r4819518 * r4819510;
double r4819523 = r4819521 * r4819522;
double r4819524 = r4819520 + r4819523;
double r4819525 = 0.0001789971;
double r4819526 = r4819522 * r4819510;
double r4819527 = r4819525 * r4819526;
double r4819528 = r4819524 + r4819527;
double r4819529 = 0.7715471019;
double r4819530 = r4819529 * r4819510;
double r4819531 = r4819507 + r4819530;
double r4819532 = 0.2909738639;
double r4819533 = r4819532 * r4819514;
double r4819534 = r4819531 + r4819533;
double r4819535 = 0.0694555761;
double r4819536 = r4819535 * r4819518;
double r4819537 = r4819534 + r4819536;
double r4819538 = 0.0140005442;
double r4819539 = r4819538 * r4819522;
double r4819540 = r4819537 + r4819539;
double r4819541 = 0.0008327945;
double r4819542 = r4819541 * r4819526;
double r4819543 = r4819540 + r4819542;
double r4819544 = 2.0;
double r4819545 = r4819544 * r4819525;
double r4819546 = r4819526 * r4819510;
double r4819547 = r4819545 * r4819546;
double r4819548 = r4819543 + r4819547;
double r4819549 = r4819528 / r4819548;
double r4819550 = r4819549 * r4819509;
return r4819550;
}
double f(double x) {
double r4819551 = x;
double r4819552 = -621.4550718514375;
bool r4819553 = r4819551 <= r4819552;
double r4819554 = 0.2514179000665375;
double r4819555 = r4819554 / r4819551;
double r4819556 = r4819551 * r4819551;
double r4819557 = r4819555 / r4819556;
double r4819558 = 0.5;
double r4819559 = r4819558 / r4819551;
double r4819560 = 0.15298196345929327;
double r4819561 = r4819560 / r4819556;
double r4819562 = r4819556 * r4819551;
double r4819563 = r4819561 / r4819562;
double r4819564 = r4819559 + r4819563;
double r4819565 = r4819557 + r4819564;
double r4819566 = 640.4757042222545;
bool r4819567 = r4819551 <= r4819566;
double r4819568 = 0.0001789971;
double r4819569 = r4819556 * r4819568;
double r4819570 = r4819556 * r4819556;
double r4819571 = r4819570 * r4819570;
double r4819572 = 0.0005064034;
double r4819573 = 0.0072644182;
double r4819574 = r4819570 * r4819573;
double r4819575 = 0.0424060604;
double r4819576 = 0.1049934947;
double r4819577 = 1.0;
double r4819578 = fma(r4819576, r4819556, r4819577);
double r4819579 = fma(r4819570, r4819575, r4819578);
double r4819580 = fma(r4819556, r4819574, r4819579);
double r4819581 = fma(r4819572, r4819571, r4819580);
double r4819582 = fma(r4819569, r4819571, r4819581);
double r4819583 = 2.0;
double r4819584 = r4819583 * r4819568;
double r4819585 = r4819556 * r4819570;
double r4819586 = r4819585 * r4819585;
double r4819587 = 0.0008327945;
double r4819588 = r4819587 * r4819570;
double r4819589 = 0.0140005442;
double r4819590 = r4819589 * r4819556;
double r4819591 = 0.0694555761;
double r4819592 = 0.2909738639;
double r4819593 = 0.7715471019;
double r4819594 = fma(r4819556, r4819593, r4819577);
double r4819595 = fma(r4819592, r4819570, r4819594);
double r4819596 = fma(r4819585, r4819591, r4819595);
double r4819597 = fma(r4819585, r4819590, r4819596);
double r4819598 = fma(r4819585, r4819588, r4819597);
double r4819599 = fma(r4819584, r4819586, r4819598);
double r4819600 = r4819582 / r4819599;
double r4819601 = r4819600 * r4819600;
double r4819602 = r4819600 * r4819601;
double r4819603 = cbrt(r4819602);
double r4819604 = r4819551 * r4819603;
double r4819605 = r4819567 ? r4819604 : r4819565;
double r4819606 = r4819553 ? r4819565 : r4819605;
return r4819606;
}



Bits error versus x
if x < -621.4550718514375 or 640.4757042222545 < x Initial program 58.2
Simplified58.2
Taylor expanded around inf 0.0
Simplified0.0
if -621.4550718514375 < x < 640.4757042222545Initial program 0.0
Simplified0.0
rmApplied add-cbrt-cube0.0
Final simplification0.0
herbie shell --seed 2019143 +o rules:numerics
(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))