\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 -3749.0578359442816 \lor \neg \left(x \le 667.593533562165362\right):\\
\;\;\;\;0.25141790006653753 \cdot \frac{1}{{x}^{3}} + \left(0.1529819634592933 \cdot \frac{1}{{x}^{5}} + 0.5 \cdot \frac{1}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\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 \left(\sqrt{\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\right)\\
\end{array}double f(double x) {
double r661 = 1.0;
double r662 = 0.1049934947;
double r663 = x;
double r664 = r663 * r663;
double r665 = r662 * r664;
double r666 = r661 + r665;
double r667 = 0.0424060604;
double r668 = r664 * r664;
double r669 = r667 * r668;
double r670 = r666 + r669;
double r671 = 0.0072644182;
double r672 = r668 * r664;
double r673 = r671 * r672;
double r674 = r670 + r673;
double r675 = 0.0005064034;
double r676 = r672 * r664;
double r677 = r675 * r676;
double r678 = r674 + r677;
double r679 = 0.0001789971;
double r680 = r676 * r664;
double r681 = r679 * r680;
double r682 = r678 + r681;
double r683 = 0.7715471019;
double r684 = r683 * r664;
double r685 = r661 + r684;
double r686 = 0.2909738639;
double r687 = r686 * r668;
double r688 = r685 + r687;
double r689 = 0.0694555761;
double r690 = r689 * r672;
double r691 = r688 + r690;
double r692 = 0.0140005442;
double r693 = r692 * r676;
double r694 = r691 + r693;
double r695 = 0.0008327945;
double r696 = r695 * r680;
double r697 = r694 + r696;
double r698 = 2.0;
double r699 = r698 * r679;
double r700 = r680 * r664;
double r701 = r699 * r700;
double r702 = r697 + r701;
double r703 = r682 / r702;
double r704 = r703 * r663;
return r704;
}
double f(double x) {
double r705 = x;
double r706 = -3749.0578359442816;
bool r707 = r705 <= r706;
double r708 = 667.5935335621654;
bool r709 = r705 <= r708;
double r710 = !r709;
bool r711 = r707 || r710;
double r712 = 0.2514179000665375;
double r713 = 1.0;
double r714 = 3.0;
double r715 = pow(r705, r714);
double r716 = r713 / r715;
double r717 = r712 * r716;
double r718 = 0.15298196345929327;
double r719 = 5.0;
double r720 = pow(r705, r719);
double r721 = r713 / r720;
double r722 = r718 * r721;
double r723 = 0.5;
double r724 = r713 / r705;
double r725 = r723 * r724;
double r726 = r722 + r725;
double r727 = r717 + r726;
double r728 = 1.0;
double r729 = 0.1049934947;
double r730 = r705 * r705;
double r731 = r729 * r730;
double r732 = r728 + r731;
double r733 = 0.0424060604;
double r734 = r730 * r730;
double r735 = r733 * r734;
double r736 = r732 + r735;
double r737 = 0.0072644182;
double r738 = r734 * r730;
double r739 = r737 * r738;
double r740 = r736 + r739;
double r741 = 0.0005064034;
double r742 = r738 * r730;
double r743 = r741 * r742;
double r744 = r740 + r743;
double r745 = 0.0001789971;
double r746 = r742 * r730;
double r747 = r745 * r746;
double r748 = r744 + r747;
double r749 = 0.7715471019;
double r750 = r749 * r730;
double r751 = r728 + r750;
double r752 = 0.2909738639;
double r753 = r752 * r734;
double r754 = r751 + r753;
double r755 = 0.0694555761;
double r756 = r755 * r738;
double r757 = r754 + r756;
double r758 = 0.0140005442;
double r759 = r758 * r742;
double r760 = r757 + r759;
double r761 = 0.0008327945;
double r762 = r761 * r746;
double r763 = r760 + r762;
double r764 = 2.0;
double r765 = r764 * r745;
double r766 = r746 * r730;
double r767 = r765 * r766;
double r768 = r763 + r767;
double r769 = r748 / r768;
double r770 = sqrt(r769);
double r771 = r770 * r705;
double r772 = r770 * r771;
double r773 = r711 ? r727 : r772;
return r773;
}



Bits error versus x
Results
if x < -3749.0578359442816 or 667.5935335621654 < x Initial program 59.2
Taylor expanded around inf 0.0
if -3749.0578359442816 < x < 667.5935335621654Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2020025
(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))