\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 -1368.36871296981394 \lor \neg \left(x \le 638.42769757308974\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{\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)} \cdot \left(\sqrt{\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)} \cdot \frac{x}{\left(\left(\left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right) \cdot 0.069455576099999999 + 1\right) + {x}^{2} \cdot \left(0.77154710189999998 + 0.29097386390000002 \cdot {x}^{2}\right)\right) + \left(\left({x}^{2} \cdot \left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right)\right) \cdot 0.014000544199999999 + {x}^{2} \cdot \left(8.32794500000000044 \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) + \left(2 \cdot 1.789971 \cdot 10^{-4}\right) \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)\right)}\right)\\
\end{array}double f(double x) {
double r176722 = 1.0;
double r176723 = 0.1049934947;
double r176724 = x;
double r176725 = r176724 * r176724;
double r176726 = r176723 * r176725;
double r176727 = r176722 + r176726;
double r176728 = 0.0424060604;
double r176729 = r176725 * r176725;
double r176730 = r176728 * r176729;
double r176731 = r176727 + r176730;
double r176732 = 0.0072644182;
double r176733 = r176729 * r176725;
double r176734 = r176732 * r176733;
double r176735 = r176731 + r176734;
double r176736 = 0.0005064034;
double r176737 = r176733 * r176725;
double r176738 = r176736 * r176737;
double r176739 = r176735 + r176738;
double r176740 = 0.0001789971;
double r176741 = r176737 * r176725;
double r176742 = r176740 * r176741;
double r176743 = r176739 + r176742;
double r176744 = 0.7715471019;
double r176745 = r176744 * r176725;
double r176746 = r176722 + r176745;
double r176747 = 0.2909738639;
double r176748 = r176747 * r176729;
double r176749 = r176746 + r176748;
double r176750 = 0.0694555761;
double r176751 = r176750 * r176733;
double r176752 = r176749 + r176751;
double r176753 = 0.0140005442;
double r176754 = r176753 * r176737;
double r176755 = r176752 + r176754;
double r176756 = 0.0008327945;
double r176757 = r176756 * r176741;
double r176758 = r176755 + r176757;
double r176759 = 2.0;
double r176760 = r176759 * r176740;
double r176761 = r176741 * r176725;
double r176762 = r176760 * r176761;
double r176763 = r176758 + r176762;
double r176764 = r176743 / r176763;
double r176765 = r176764 * r176724;
return r176765;
}
double f(double x) {
double r176766 = x;
double r176767 = -1368.368712969814;
bool r176768 = r176766 <= r176767;
double r176769 = 638.4276975730897;
bool r176770 = r176766 <= r176769;
double r176771 = !r176770;
bool r176772 = r176768 || r176771;
double r176773 = 0.2514179000665375;
double r176774 = 1.0;
double r176775 = 3.0;
double r176776 = pow(r176766, r176775);
double r176777 = r176774 / r176776;
double r176778 = r176773 * r176777;
double r176779 = 0.15298196345929327;
double r176780 = 5.0;
double r176781 = pow(r176766, r176780);
double r176782 = r176774 / r176781;
double r176783 = r176779 * r176782;
double r176784 = 0.5;
double r176785 = r176774 / r176766;
double r176786 = r176784 * r176785;
double r176787 = r176783 + r176786;
double r176788 = r176778 + r176787;
double r176789 = 1.0;
double r176790 = 0.1049934947;
double r176791 = r176766 * r176766;
double r176792 = r176790 * r176791;
double r176793 = r176789 + r176792;
double r176794 = 0.0424060604;
double r176795 = r176791 * r176791;
double r176796 = r176794 * r176795;
double r176797 = r176793 + r176796;
double r176798 = 0.0072644182;
double r176799 = r176795 * r176791;
double r176800 = r176798 * r176799;
double r176801 = r176797 + r176800;
double r176802 = 0.0005064034;
double r176803 = r176799 * r176791;
double r176804 = r176802 * r176803;
double r176805 = r176801 + r176804;
double r176806 = 0.0001789971;
double r176807 = r176803 * r176791;
double r176808 = r176806 * r176807;
double r176809 = r176805 + r176808;
double r176810 = sqrt(r176809);
double r176811 = 2.0;
double r176812 = pow(r176766, r176811);
double r176813 = r176766 * r176776;
double r176814 = r176812 * r176813;
double r176815 = 0.0694555761;
double r176816 = r176814 * r176815;
double r176817 = r176816 + r176789;
double r176818 = 0.7715471019;
double r176819 = 0.2909738639;
double r176820 = r176819 * r176812;
double r176821 = r176818 + r176820;
double r176822 = r176812 * r176821;
double r176823 = r176817 + r176822;
double r176824 = r176812 * r176814;
double r176825 = 0.0140005442;
double r176826 = r176824 * r176825;
double r176827 = 0.0008327945;
double r176828 = r176827 * r176803;
double r176829 = 2.0;
double r176830 = r176829 * r176806;
double r176831 = r176830 * r176807;
double r176832 = r176828 + r176831;
double r176833 = r176812 * r176832;
double r176834 = r176826 + r176833;
double r176835 = r176823 + r176834;
double r176836 = r176766 / r176835;
double r176837 = r176810 * r176836;
double r176838 = r176810 * r176837;
double r176839 = r176772 ? r176788 : r176838;
return r176839;
}



Bits error versus x
Results
if x < -1368.368712969814 or 638.4276975730897 < x Initial program 59.7
Taylor expanded around inf 0.0
if -1368.368712969814 < x < 638.4276975730897Initial program 0.0
rmApplied div-inv0.0
Applied associate-*l*0.0
Simplified0.0
rmApplied add-sqr-sqrt0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2020021
(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))