\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 -977.674404093754788 \lor \neg \left(x \le 905.501617915186557\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{1}{\frac{\left(\left({x}^{2} \cdot \left({x}^{2} \cdot \left({x}^{2} \cdot \left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right)\right)\right)\right) \cdot \left(2 \cdot 1.789971 \cdot 10^{-4}\right) + \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)\right) + {x}^{2} \cdot \left(0.014000544199999999 \cdot \left(\left(\left(x \cdot x\right) \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)\right) + 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)\right)}{\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}}} \cdot x\\
\end{array}double f(double x) {
double r205977 = 1.0;
double r205978 = 0.1049934947;
double r205979 = x;
double r205980 = r205979 * r205979;
double r205981 = r205978 * r205980;
double r205982 = r205977 + r205981;
double r205983 = 0.0424060604;
double r205984 = r205980 * r205980;
double r205985 = r205983 * r205984;
double r205986 = r205982 + r205985;
double r205987 = 0.0072644182;
double r205988 = r205984 * r205980;
double r205989 = r205987 * r205988;
double r205990 = r205986 + r205989;
double r205991 = 0.0005064034;
double r205992 = r205988 * r205980;
double r205993 = r205991 * r205992;
double r205994 = r205990 + r205993;
double r205995 = 0.0001789971;
double r205996 = r205992 * r205980;
double r205997 = r205995 * r205996;
double r205998 = r205994 + r205997;
double r205999 = 0.7715471019;
double r206000 = r205999 * r205980;
double r206001 = r205977 + r206000;
double r206002 = 0.2909738639;
double r206003 = r206002 * r205984;
double r206004 = r206001 + r206003;
double r206005 = 0.0694555761;
double r206006 = r206005 * r205988;
double r206007 = r206004 + r206006;
double r206008 = 0.0140005442;
double r206009 = r206008 * r205992;
double r206010 = r206007 + r206009;
double r206011 = 0.0008327945;
double r206012 = r206011 * r205996;
double r206013 = r206010 + r206012;
double r206014 = 2.0;
double r206015 = r206014 * r205995;
double r206016 = r205996 * r205980;
double r206017 = r206015 * r206016;
double r206018 = r206013 + r206017;
double r206019 = r205998 / r206018;
double r206020 = r206019 * r205979;
return r206020;
}
double f(double x) {
double r206021 = x;
double r206022 = -977.6744040937548;
bool r206023 = r206021 <= r206022;
double r206024 = 905.5016179151866;
bool r206025 = r206021 <= r206024;
double r206026 = !r206025;
bool r206027 = r206023 || r206026;
double r206028 = 0.2514179000665375;
double r206029 = 1.0;
double r206030 = 3.0;
double r206031 = pow(r206021, r206030);
double r206032 = r206029 / r206031;
double r206033 = r206028 * r206032;
double r206034 = 0.15298196345929327;
double r206035 = 5.0;
double r206036 = pow(r206021, r206035);
double r206037 = r206029 / r206036;
double r206038 = r206034 * r206037;
double r206039 = 0.5;
double r206040 = r206029 / r206021;
double r206041 = r206039 * r206040;
double r206042 = r206038 + r206041;
double r206043 = r206033 + r206042;
double r206044 = 2.0;
double r206045 = pow(r206021, r206044);
double r206046 = r206021 * r206031;
double r206047 = r206045 * r206046;
double r206048 = r206045 * r206047;
double r206049 = r206045 * r206048;
double r206050 = r206045 * r206049;
double r206051 = 2.0;
double r206052 = 0.0001789971;
double r206053 = r206051 * r206052;
double r206054 = r206050 * r206053;
double r206055 = 0.0694555761;
double r206056 = r206047 * r206055;
double r206057 = 1.0;
double r206058 = r206056 + r206057;
double r206059 = 0.7715471019;
double r206060 = 0.2909738639;
double r206061 = r206060 * r206045;
double r206062 = r206059 + r206061;
double r206063 = r206045 * r206062;
double r206064 = r206058 + r206063;
double r206065 = r206054 + r206064;
double r206066 = 0.0140005442;
double r206067 = r206021 * r206021;
double r206068 = r206067 * r206021;
double r206069 = r206068 * r206068;
double r206070 = r206066 * r206069;
double r206071 = 0.0008327945;
double r206072 = r206067 * r206067;
double r206073 = r206072 * r206067;
double r206074 = r206073 * r206067;
double r206075 = r206071 * r206074;
double r206076 = r206070 + r206075;
double r206077 = r206045 * r206076;
double r206078 = r206065 + r206077;
double r206079 = r206049 * r206052;
double r206080 = 0.0072644182;
double r206081 = r206047 * r206080;
double r206082 = r206081 + r206057;
double r206083 = 0.1049934947;
double r206084 = 0.0424060604;
double r206085 = r206084 * r206045;
double r206086 = r206083 + r206085;
double r206087 = r206045 * r206086;
double r206088 = r206082 + r206087;
double r206089 = r206079 + r206088;
double r206090 = 0.0005064034;
double r206091 = r206048 * r206090;
double r206092 = r206089 + r206091;
double r206093 = r206078 / r206092;
double r206094 = r206029 / r206093;
double r206095 = r206094 * r206021;
double r206096 = r206027 ? r206043 : r206095;
return r206096;
}



Bits error versus x
Results
if x < -977.6744040937548 or 905.5016179151866 < x Initial program 59.2
Taylor expanded around inf 0.0
if -977.6744040937548 < x < 905.5016179151866Initial program 0.0
rmApplied clear-num0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2020100
(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))