\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 -784.29567757355517 \lor \neg \left(x \le 5774466.4635821618\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}:\\
\;\;\;\;\frac{\left(\left(\left({x}^{2} \cdot \left({x}^{2} \cdot \left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right)\right)\right) \cdot 1.789971 \cdot 10^{-4} + \left(\left(\left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right) \cdot 0.00726441819999999999 + 1\right) + {x}^{2} \cdot \left(0.1049934947 + 0.042406060400000001 \cdot {x}^{2}\right)\right)\right) + \left({x}^{2} \cdot \left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right)\right) \cdot 5.0640340000000002 \cdot 10^{-4}\right) \cdot x}{\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)}\\
\end{array}double f(double x) {
double r268720 = 1.0;
double r268721 = 0.1049934947;
double r268722 = x;
double r268723 = r268722 * r268722;
double r268724 = r268721 * r268723;
double r268725 = r268720 + r268724;
double r268726 = 0.0424060604;
double r268727 = r268723 * r268723;
double r268728 = r268726 * r268727;
double r268729 = r268725 + r268728;
double r268730 = 0.0072644182;
double r268731 = r268727 * r268723;
double r268732 = r268730 * r268731;
double r268733 = r268729 + r268732;
double r268734 = 0.0005064034;
double r268735 = r268731 * r268723;
double r268736 = r268734 * r268735;
double r268737 = r268733 + r268736;
double r268738 = 0.0001789971;
double r268739 = r268735 * r268723;
double r268740 = r268738 * r268739;
double r268741 = r268737 + r268740;
double r268742 = 0.7715471019;
double r268743 = r268742 * r268723;
double r268744 = r268720 + r268743;
double r268745 = 0.2909738639;
double r268746 = r268745 * r268727;
double r268747 = r268744 + r268746;
double r268748 = 0.0694555761;
double r268749 = r268748 * r268731;
double r268750 = r268747 + r268749;
double r268751 = 0.0140005442;
double r268752 = r268751 * r268735;
double r268753 = r268750 + r268752;
double r268754 = 0.0008327945;
double r268755 = r268754 * r268739;
double r268756 = r268753 + r268755;
double r268757 = 2.0;
double r268758 = r268757 * r268738;
double r268759 = r268739 * r268723;
double r268760 = r268758 * r268759;
double r268761 = r268756 + r268760;
double r268762 = r268741 / r268761;
double r268763 = r268762 * r268722;
return r268763;
}
double f(double x) {
double r268764 = x;
double r268765 = -784.2956775735552;
bool r268766 = r268764 <= r268765;
double r268767 = 5774466.463582162;
bool r268768 = r268764 <= r268767;
double r268769 = !r268768;
bool r268770 = r268766 || r268769;
double r268771 = 0.2514179000665375;
double r268772 = 1.0;
double r268773 = 3.0;
double r268774 = pow(r268764, r268773);
double r268775 = r268772 / r268774;
double r268776 = r268771 * r268775;
double r268777 = 0.15298196345929327;
double r268778 = 5.0;
double r268779 = pow(r268764, r268778);
double r268780 = r268772 / r268779;
double r268781 = r268777 * r268780;
double r268782 = 0.5;
double r268783 = r268772 / r268764;
double r268784 = r268782 * r268783;
double r268785 = r268781 + r268784;
double r268786 = r268776 + r268785;
double r268787 = 2.0;
double r268788 = pow(r268764, r268787);
double r268789 = r268764 * r268774;
double r268790 = r268788 * r268789;
double r268791 = r268788 * r268790;
double r268792 = r268788 * r268791;
double r268793 = 0.0001789971;
double r268794 = r268792 * r268793;
double r268795 = 0.0072644182;
double r268796 = r268790 * r268795;
double r268797 = 1.0;
double r268798 = r268796 + r268797;
double r268799 = 0.1049934947;
double r268800 = 0.0424060604;
double r268801 = r268800 * r268788;
double r268802 = r268799 + r268801;
double r268803 = r268788 * r268802;
double r268804 = r268798 + r268803;
double r268805 = r268794 + r268804;
double r268806 = 0.0005064034;
double r268807 = r268791 * r268806;
double r268808 = r268805 + r268807;
double r268809 = r268808 * r268764;
double r268810 = 0.7715471019;
double r268811 = r268764 * r268764;
double r268812 = r268810 * r268811;
double r268813 = r268797 + r268812;
double r268814 = 0.2909738639;
double r268815 = r268811 * r268811;
double r268816 = r268814 * r268815;
double r268817 = r268813 + r268816;
double r268818 = 0.0694555761;
double r268819 = r268815 * r268811;
double r268820 = r268818 * r268819;
double r268821 = r268817 + r268820;
double r268822 = 0.0140005442;
double r268823 = r268819 * r268811;
double r268824 = r268822 * r268823;
double r268825 = r268821 + r268824;
double r268826 = 0.0008327945;
double r268827 = r268823 * r268811;
double r268828 = r268826 * r268827;
double r268829 = r268825 + r268828;
double r268830 = 2.0;
double r268831 = r268830 * r268793;
double r268832 = r268827 * r268811;
double r268833 = r268831 * r268832;
double r268834 = r268829 + r268833;
double r268835 = r268809 / r268834;
double r268836 = r268770 ? r268786 : r268835;
return r268836;
}



Bits error versus x
Results
if x < -784.2956775735552 or 5774466.463582162 < x Initial program 59.8
Taylor expanded around inf 0.0
if -784.2956775735552 < x < 5774466.463582162Initial program 0.0
rmApplied associate-*l/0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2020083
(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))