\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 -9087113337.796595:\\
\;\;\;\;\frac{0.15298196345929327}{{x}^{5}} + \left(\frac{0.2514179000665375}{x \cdot \left(x \cdot x\right)} + \frac{0.5}{x}\right)\\
\mathbf{elif}\;x \le 720.524117169935:\\
\;\;\;\;\frac{\left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) \cdot 0.0001789971 + \left(\left(0.0072644182 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + \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)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot 0.0005064034\right)\right) \cdot x}{\left(0.0008327945 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.0694555761 + \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.2909738639 + \left(1 + \left(x \cdot x\right) \cdot 0.7715471019\right)\right)\right) + 0.0140005442 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right) \cdot \left(0.0008327945 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.0694555761 + \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.2909738639 + \left(1 + \left(x \cdot x\right) \cdot 0.7715471019\right)\right)\right) + 0.0140005442 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right) - \left(\left(0.0001789971 \cdot 2\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(\left(0.0001789971 \cdot 2\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right)\right)} \cdot \left(\left(0.0008327945 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.0694555761 + \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.2909738639 + \left(1 + \left(x \cdot x\right) \cdot 0.7715471019\right)\right)\right) + 0.0140005442 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right) - \left(0.0001789971 \cdot 2\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{0.15298196345929327}{{x}^{5}} + \left(\frac{0.2514179000665375}{x \cdot \left(x \cdot x\right)} + \frac{0.5}{x}\right)\\
\end{array}double f(double x) {
double r6279712 = 1.0;
double r6279713 = 0.1049934947;
double r6279714 = x;
double r6279715 = r6279714 * r6279714;
double r6279716 = r6279713 * r6279715;
double r6279717 = r6279712 + r6279716;
double r6279718 = 0.0424060604;
double r6279719 = r6279715 * r6279715;
double r6279720 = r6279718 * r6279719;
double r6279721 = r6279717 + r6279720;
double r6279722 = 0.0072644182;
double r6279723 = r6279719 * r6279715;
double r6279724 = r6279722 * r6279723;
double r6279725 = r6279721 + r6279724;
double r6279726 = 0.0005064034;
double r6279727 = r6279723 * r6279715;
double r6279728 = r6279726 * r6279727;
double r6279729 = r6279725 + r6279728;
double r6279730 = 0.0001789971;
double r6279731 = r6279727 * r6279715;
double r6279732 = r6279730 * r6279731;
double r6279733 = r6279729 + r6279732;
double r6279734 = 0.7715471019;
double r6279735 = r6279734 * r6279715;
double r6279736 = r6279712 + r6279735;
double r6279737 = 0.2909738639;
double r6279738 = r6279737 * r6279719;
double r6279739 = r6279736 + r6279738;
double r6279740 = 0.0694555761;
double r6279741 = r6279740 * r6279723;
double r6279742 = r6279739 + r6279741;
double r6279743 = 0.0140005442;
double r6279744 = r6279743 * r6279727;
double r6279745 = r6279742 + r6279744;
double r6279746 = 0.0008327945;
double r6279747 = r6279746 * r6279731;
double r6279748 = r6279745 + r6279747;
double r6279749 = 2.0;
double r6279750 = r6279749 * r6279730;
double r6279751 = r6279731 * r6279715;
double r6279752 = r6279750 * r6279751;
double r6279753 = r6279748 + r6279752;
double r6279754 = r6279733 / r6279753;
double r6279755 = r6279754 * r6279714;
return r6279755;
}
double f(double x) {
double r6279756 = x;
double r6279757 = -9087113337.796595;
bool r6279758 = r6279756 <= r6279757;
double r6279759 = 0.15298196345929327;
double r6279760 = 5.0;
double r6279761 = pow(r6279756, r6279760);
double r6279762 = r6279759 / r6279761;
double r6279763 = 0.2514179000665375;
double r6279764 = r6279756 * r6279756;
double r6279765 = r6279756 * r6279764;
double r6279766 = r6279763 / r6279765;
double r6279767 = 0.5;
double r6279768 = r6279767 / r6279756;
double r6279769 = r6279766 + r6279768;
double r6279770 = r6279762 + r6279769;
double r6279771 = 720.524117169935;
bool r6279772 = r6279756 <= r6279771;
double r6279773 = r6279764 * r6279764;
double r6279774 = r6279764 * r6279773;
double r6279775 = r6279764 * r6279774;
double r6279776 = r6279764 * r6279775;
double r6279777 = 0.0001789971;
double r6279778 = r6279776 * r6279777;
double r6279779 = 0.0072644182;
double r6279780 = r6279779 * r6279774;
double r6279781 = 1.0;
double r6279782 = 0.1049934947;
double r6279783 = r6279782 * r6279764;
double r6279784 = r6279781 + r6279783;
double r6279785 = 0.0424060604;
double r6279786 = r6279785 * r6279773;
double r6279787 = r6279784 + r6279786;
double r6279788 = r6279780 + r6279787;
double r6279789 = 0.0005064034;
double r6279790 = r6279775 * r6279789;
double r6279791 = r6279788 + r6279790;
double r6279792 = r6279778 + r6279791;
double r6279793 = r6279792 * r6279756;
double r6279794 = 0.0008327945;
double r6279795 = r6279794 * r6279776;
double r6279796 = 0.0694555761;
double r6279797 = r6279774 * r6279796;
double r6279798 = 0.2909738639;
double r6279799 = r6279773 * r6279798;
double r6279800 = 0.7715471019;
double r6279801 = r6279764 * r6279800;
double r6279802 = r6279781 + r6279801;
double r6279803 = r6279799 + r6279802;
double r6279804 = r6279797 + r6279803;
double r6279805 = 0.0140005442;
double r6279806 = r6279805 * r6279775;
double r6279807 = r6279804 + r6279806;
double r6279808 = r6279795 + r6279807;
double r6279809 = r6279808 * r6279808;
double r6279810 = 2.0;
double r6279811 = r6279777 * r6279810;
double r6279812 = r6279776 * r6279764;
double r6279813 = r6279811 * r6279812;
double r6279814 = r6279813 * r6279813;
double r6279815 = r6279809 - r6279814;
double r6279816 = r6279793 / r6279815;
double r6279817 = r6279808 - r6279813;
double r6279818 = r6279816 * r6279817;
double r6279819 = r6279772 ? r6279818 : r6279770;
double r6279820 = r6279758 ? r6279770 : r6279819;
return r6279820;
}



Bits error versus x
Results
if x < -9087113337.796595 or 720.524117169935 < x Initial program 59.0
Taylor expanded around inf 0.0
Simplified0.0
if -9087113337.796595 < x < 720.524117169935Initial program 0.0
rmApplied associate-*l/0.0
rmApplied flip-+0.0
Applied associate-/r/0.0
Final simplification0.0
herbie shell --seed 2019141
(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))