\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 -1207.5387837908731:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\mathbf{elif}\;x \le 727.2141028040012:\\
\;\;\;\;x \cdot \frac{\frac{\left(\left(0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.1049934947 \cdot \left(x \cdot x\right) + 1\right)\right)\right) + \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 0.0005064034\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)}{\sqrt{\left(\left(x \cdot x\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) \cdot \left(0.0001789971 \cdot 2\right) + \left(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) + \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 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}}}{\sqrt{\left(\left(x \cdot x\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) \cdot \left(0.0001789971 \cdot 2\right) + \left(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) + \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 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\end{array}double f(double x) {
double r82746754 = 1.0;
double r82746755 = 0.1049934947;
double r82746756 = x;
double r82746757 = r82746756 * r82746756;
double r82746758 = r82746755 * r82746757;
double r82746759 = r82746754 + r82746758;
double r82746760 = 0.0424060604;
double r82746761 = r82746757 * r82746757;
double r82746762 = r82746760 * r82746761;
double r82746763 = r82746759 + r82746762;
double r82746764 = 0.0072644182;
double r82746765 = r82746761 * r82746757;
double r82746766 = r82746764 * r82746765;
double r82746767 = r82746763 + r82746766;
double r82746768 = 0.0005064034;
double r82746769 = r82746765 * r82746757;
double r82746770 = r82746768 * r82746769;
double r82746771 = r82746767 + r82746770;
double r82746772 = 0.0001789971;
double r82746773 = r82746769 * r82746757;
double r82746774 = r82746772 * r82746773;
double r82746775 = r82746771 + r82746774;
double r82746776 = 0.7715471019;
double r82746777 = r82746776 * r82746757;
double r82746778 = r82746754 + r82746777;
double r82746779 = 0.2909738639;
double r82746780 = r82746779 * r82746761;
double r82746781 = r82746778 + r82746780;
double r82746782 = 0.0694555761;
double r82746783 = r82746782 * r82746765;
double r82746784 = r82746781 + r82746783;
double r82746785 = 0.0140005442;
double r82746786 = r82746785 * r82746769;
double r82746787 = r82746784 + r82746786;
double r82746788 = 0.0008327945;
double r82746789 = r82746788 * r82746773;
double r82746790 = r82746787 + r82746789;
double r82746791 = 2.0;
double r82746792 = r82746791 * r82746772;
double r82746793 = r82746773 * r82746757;
double r82746794 = r82746792 * r82746793;
double r82746795 = r82746790 + r82746794;
double r82746796 = r82746775 / r82746795;
double r82746797 = r82746796 * r82746756;
return r82746797;
}
double f(double x) {
double r82746798 = x;
double r82746799 = -1207.5387837908731;
bool r82746800 = r82746798 <= r82746799;
double r82746801 = 0.5;
double r82746802 = 0.2514179000665375;
double r82746803 = r82746798 * r82746798;
double r82746804 = r82746802 / r82746803;
double r82746805 = r82746801 + r82746804;
double r82746806 = r82746805 / r82746798;
double r82746807 = 0.15298196345929327;
double r82746808 = 5.0;
double r82746809 = pow(r82746798, r82746808);
double r82746810 = r82746807 / r82746809;
double r82746811 = r82746806 + r82746810;
double r82746812 = 727.2141028040012;
bool r82746813 = r82746798 <= r82746812;
double r82746814 = 0.0072644182;
double r82746815 = r82746803 * r82746803;
double r82746816 = r82746815 * r82746803;
double r82746817 = r82746814 * r82746816;
double r82746818 = 0.0424060604;
double r82746819 = r82746818 * r82746815;
double r82746820 = 0.1049934947;
double r82746821 = r82746820 * r82746803;
double r82746822 = 1.0;
double r82746823 = r82746821 + r82746822;
double r82746824 = r82746819 + r82746823;
double r82746825 = r82746817 + r82746824;
double r82746826 = r82746816 * r82746803;
double r82746827 = 0.0005064034;
double r82746828 = r82746826 * r82746827;
double r82746829 = r82746825 + r82746828;
double r82746830 = 0.0001789971;
double r82746831 = r82746826 * r82746803;
double r82746832 = r82746830 * r82746831;
double r82746833 = r82746829 + r82746832;
double r82746834 = r82746803 * r82746831;
double r82746835 = 2.0;
double r82746836 = r82746830 * r82746835;
double r82746837 = r82746834 * r82746836;
double r82746838 = 0.0008327945;
double r82746839 = r82746838 * r82746831;
double r82746840 = 0.0140005442;
double r82746841 = r82746826 * r82746840;
double r82746842 = 0.0694555761;
double r82746843 = r82746816 * r82746842;
double r82746844 = 0.2909738639;
double r82746845 = r82746844 * r82746815;
double r82746846 = 0.7715471019;
double r82746847 = r82746846 * r82746803;
double r82746848 = r82746847 + r82746822;
double r82746849 = r82746845 + r82746848;
double r82746850 = r82746843 + r82746849;
double r82746851 = r82746841 + r82746850;
double r82746852 = r82746839 + r82746851;
double r82746853 = r82746837 + r82746852;
double r82746854 = sqrt(r82746853);
double r82746855 = r82746833 / r82746854;
double r82746856 = r82746855 / r82746854;
double r82746857 = r82746798 * r82746856;
double r82746858 = r82746813 ? r82746857 : r82746811;
double r82746859 = r82746800 ? r82746811 : r82746858;
return r82746859;
}



Bits error versus x
Results
if x < -1207.5387837908731 or 727.2141028040012 < x Initial program 58.0
rmApplied add-sqr-sqrt58.0
Applied associate-/r*58.0
Taylor expanded around -inf 0.0
Simplified0.0
if -1207.5387837908731 < x < 727.2141028040012Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied associate-/r*0.0
Final simplification0.0
herbie shell --seed 2019107 +o rules:numerics
(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))