\frac{\left(\left(\left(\left(1 + 0.1049934947 \cdot \left(x \cdot x\right)\right) + 0.042406060400000001 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.00726441819999999999 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 5.0640340000000002 \cdot 10^{-4} \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) + 1.789971 \cdot 10^{-4} \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.77154710189999998 \cdot \left(x \cdot x\right)\right) + 0.29097386390000002 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.069455576099999999 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.014000544199999999 \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) + 8.32794500000000044 \cdot 10^{-4} \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 1.789971 \cdot 10^{-4}\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 -266320.39092243672 \lor \neg \left(x \le 4432.6458967592771\right):\\
\;\;\;\;\left(\frac{0.25141790006653753}{{x}^{3}} + \frac{0.5}{x}\right) + \frac{0.1529819634592933}{{x}^{5}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\left(2 \cdot 1.789971 \cdot 10^{-4}\right) \cdot {x}^{12} + \left({x}^{8} \cdot 0.014000544199999999 + \left(\left(0.069455576099999999 \cdot {x}^{6} + \left(\left(1 + 0.77154710189999998 \cdot \left(x \cdot x\right)\right) + {x}^{4} \cdot 0.29097386390000002\right)\right) + 8.32794500000000044 \cdot 10^{-4} \cdot \left({x}^{6} \cdot {x}^{4}\right)\right)\right)}{\left(1.789971 \cdot 10^{-4} \cdot \left({x}^{6} \cdot {x}^{4}\right) + \left({x}^{6} \cdot 0.00726441819999999999 + \left({x}^{4} \cdot 0.042406060400000001 + \left(1 + \left(\sqrt[3]{0.1049934947 \cdot \left(x \cdot x\right)} \cdot \sqrt[3]{0.1049934947 \cdot \left(x \cdot x\right)}\right) \cdot \sqrt[3]{0.1049934947 \cdot \left(x \cdot x\right)}\right)\right)\right)\right) + 5.0640340000000002 \cdot 10^{-4} \cdot {x}^{8}}} \cdot x\\
\end{array}double f(double x) {
double r563553 = 1.0;
double r563554 = 0.1049934947;
double r563555 = x;
double r563556 = r563555 * r563555;
double r563557 = r563554 * r563556;
double r563558 = r563553 + r563557;
double r563559 = 0.0424060604;
double r563560 = r563556 * r563556;
double r563561 = r563559 * r563560;
double r563562 = r563558 + r563561;
double r563563 = 0.0072644182;
double r563564 = r563560 * r563556;
double r563565 = r563563 * r563564;
double r563566 = r563562 + r563565;
double r563567 = 0.0005064034;
double r563568 = r563564 * r563556;
double r563569 = r563567 * r563568;
double r563570 = r563566 + r563569;
double r563571 = 0.0001789971;
double r563572 = r563568 * r563556;
double r563573 = r563571 * r563572;
double r563574 = r563570 + r563573;
double r563575 = 0.7715471019;
double r563576 = r563575 * r563556;
double r563577 = r563553 + r563576;
double r563578 = 0.2909738639;
double r563579 = r563578 * r563560;
double r563580 = r563577 + r563579;
double r563581 = 0.0694555761;
double r563582 = r563581 * r563564;
double r563583 = r563580 + r563582;
double r563584 = 0.0140005442;
double r563585 = r563584 * r563568;
double r563586 = r563583 + r563585;
double r563587 = 0.0008327945;
double r563588 = r563587 * r563572;
double r563589 = r563586 + r563588;
double r563590 = 2.0;
double r563591 = r563590 * r563571;
double r563592 = r563572 * r563556;
double r563593 = r563591 * r563592;
double r563594 = r563589 + r563593;
double r563595 = r563574 / r563594;
double r563596 = r563595 * r563555;
return r563596;
}
double f(double x) {
double r563597 = x;
double r563598 = -266320.3909224367;
bool r563599 = r563597 <= r563598;
double r563600 = 4432.645896759277;
bool r563601 = r563597 <= r563600;
double r563602 = !r563601;
bool r563603 = r563599 || r563602;
double r563604 = 0.2514179000665375;
double r563605 = 3.0;
double r563606 = pow(r563597, r563605);
double r563607 = r563604 / r563606;
double r563608 = 0.5;
double r563609 = r563608 / r563597;
double r563610 = r563607 + r563609;
double r563611 = 0.15298196345929327;
double r563612 = 5.0;
double r563613 = pow(r563597, r563612);
double r563614 = r563611 / r563613;
double r563615 = r563610 + r563614;
double r563616 = 1.0;
double r563617 = 2.0;
double r563618 = 0.0001789971;
double r563619 = r563617 * r563618;
double r563620 = 12.0;
double r563621 = pow(r563597, r563620);
double r563622 = r563619 * r563621;
double r563623 = 8.0;
double r563624 = pow(r563597, r563623);
double r563625 = 0.0140005442;
double r563626 = r563624 * r563625;
double r563627 = 0.0694555761;
double r563628 = 6.0;
double r563629 = pow(r563597, r563628);
double r563630 = r563627 * r563629;
double r563631 = 1.0;
double r563632 = 0.7715471019;
double r563633 = r563597 * r563597;
double r563634 = r563632 * r563633;
double r563635 = r563631 + r563634;
double r563636 = 4.0;
double r563637 = pow(r563597, r563636);
double r563638 = 0.2909738639;
double r563639 = r563637 * r563638;
double r563640 = r563635 + r563639;
double r563641 = r563630 + r563640;
double r563642 = 0.0008327945;
double r563643 = r563629 * r563637;
double r563644 = r563642 * r563643;
double r563645 = r563641 + r563644;
double r563646 = r563626 + r563645;
double r563647 = r563622 + r563646;
double r563648 = r563618 * r563643;
double r563649 = 0.0072644182;
double r563650 = r563629 * r563649;
double r563651 = 0.0424060604;
double r563652 = r563637 * r563651;
double r563653 = 0.1049934947;
double r563654 = r563653 * r563633;
double r563655 = cbrt(r563654);
double r563656 = r563655 * r563655;
double r563657 = r563656 * r563655;
double r563658 = r563631 + r563657;
double r563659 = r563652 + r563658;
double r563660 = r563650 + r563659;
double r563661 = r563648 + r563660;
double r563662 = 0.0005064034;
double r563663 = r563662 * r563624;
double r563664 = r563661 + r563663;
double r563665 = r563647 / r563664;
double r563666 = r563616 / r563665;
double r563667 = r563666 * r563597;
double r563668 = r563603 ? r563615 : r563667;
return r563668;
}



Bits error versus x
Results
if x < -266320.3909224367 or 4432.645896759277 < x Initial program 59.2
Taylor expanded around inf 0.0
Simplified0.0
if -266320.3909224367 < x < 4432.645896759277Initial program 0.0
rmApplied clear-num0.0
Simplified0.0
rmApplied add-cube-cbrt0.0
Final simplification0.0
herbie shell --seed 2020042
(FPCore (x)
:name "Jmat.Real.dawson"
:precision binary64
(* (/ (+ (+ (+ (+ (+ 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))