\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 -5.395689961147953 \cdot 10^{+18}:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\mathbf{elif}\;x \le 14413682.003027312:\\
\;\;\;\;\left(x \cdot \frac{1}{\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)}\right) \cdot \left(\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)\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 r37301013 = 1.0;
double r37301014 = 0.1049934947;
double r37301015 = x;
double r37301016 = r37301015 * r37301015;
double r37301017 = r37301014 * r37301016;
double r37301018 = r37301013 + r37301017;
double r37301019 = 0.0424060604;
double r37301020 = r37301016 * r37301016;
double r37301021 = r37301019 * r37301020;
double r37301022 = r37301018 + r37301021;
double r37301023 = 0.0072644182;
double r37301024 = r37301020 * r37301016;
double r37301025 = r37301023 * r37301024;
double r37301026 = r37301022 + r37301025;
double r37301027 = 0.0005064034;
double r37301028 = r37301024 * r37301016;
double r37301029 = r37301027 * r37301028;
double r37301030 = r37301026 + r37301029;
double r37301031 = 0.0001789971;
double r37301032 = r37301028 * r37301016;
double r37301033 = r37301031 * r37301032;
double r37301034 = r37301030 + r37301033;
double r37301035 = 0.7715471019;
double r37301036 = r37301035 * r37301016;
double r37301037 = r37301013 + r37301036;
double r37301038 = 0.2909738639;
double r37301039 = r37301038 * r37301020;
double r37301040 = r37301037 + r37301039;
double r37301041 = 0.0694555761;
double r37301042 = r37301041 * r37301024;
double r37301043 = r37301040 + r37301042;
double r37301044 = 0.0140005442;
double r37301045 = r37301044 * r37301028;
double r37301046 = r37301043 + r37301045;
double r37301047 = 0.0008327945;
double r37301048 = r37301047 * r37301032;
double r37301049 = r37301046 + r37301048;
double r37301050 = 2.0;
double r37301051 = r37301050 * r37301031;
double r37301052 = r37301032 * r37301016;
double r37301053 = r37301051 * r37301052;
double r37301054 = r37301049 + r37301053;
double r37301055 = r37301034 / r37301054;
double r37301056 = r37301055 * r37301015;
return r37301056;
}
double f(double x) {
double r37301057 = x;
double r37301058 = -5.395689961147953e+18;
bool r37301059 = r37301057 <= r37301058;
double r37301060 = 0.5;
double r37301061 = 0.2514179000665375;
double r37301062 = r37301057 * r37301057;
double r37301063 = r37301061 / r37301062;
double r37301064 = r37301060 + r37301063;
double r37301065 = r37301064 / r37301057;
double r37301066 = 0.15298196345929327;
double r37301067 = 5.0;
double r37301068 = pow(r37301057, r37301067);
double r37301069 = r37301066 / r37301068;
double r37301070 = r37301065 + r37301069;
double r37301071 = 14413682.003027312;
bool r37301072 = r37301057 <= r37301071;
double r37301073 = 1.0;
double r37301074 = r37301062 * r37301062;
double r37301075 = r37301074 * r37301062;
double r37301076 = r37301075 * r37301062;
double r37301077 = r37301076 * r37301062;
double r37301078 = r37301062 * r37301077;
double r37301079 = 0.0001789971;
double r37301080 = 2.0;
double r37301081 = r37301079 * r37301080;
double r37301082 = r37301078 * r37301081;
double r37301083 = 0.0008327945;
double r37301084 = r37301083 * r37301077;
double r37301085 = 0.0140005442;
double r37301086 = r37301076 * r37301085;
double r37301087 = 0.0694555761;
double r37301088 = r37301075 * r37301087;
double r37301089 = 0.2909738639;
double r37301090 = r37301089 * r37301074;
double r37301091 = 0.7715471019;
double r37301092 = r37301091 * r37301062;
double r37301093 = r37301092 + r37301073;
double r37301094 = r37301090 + r37301093;
double r37301095 = r37301088 + r37301094;
double r37301096 = r37301086 + r37301095;
double r37301097 = r37301084 + r37301096;
double r37301098 = r37301082 + r37301097;
double r37301099 = r37301073 / r37301098;
double r37301100 = r37301057 * r37301099;
double r37301101 = 0.0072644182;
double r37301102 = r37301101 * r37301075;
double r37301103 = 0.0424060604;
double r37301104 = r37301103 * r37301074;
double r37301105 = 0.1049934947;
double r37301106 = r37301105 * r37301062;
double r37301107 = r37301106 + r37301073;
double r37301108 = r37301104 + r37301107;
double r37301109 = r37301102 + r37301108;
double r37301110 = 0.0005064034;
double r37301111 = r37301076 * r37301110;
double r37301112 = r37301109 + r37301111;
double r37301113 = r37301079 * r37301077;
double r37301114 = r37301112 + r37301113;
double r37301115 = r37301100 * r37301114;
double r37301116 = r37301072 ? r37301115 : r37301070;
double r37301117 = r37301059 ? r37301070 : r37301116;
return r37301117;
}



Bits error versus x
Results
if x < -5.395689961147953e+18 or 14413682.003027312 < x Initial program 60.5
rmApplied div-inv60.5
Applied associate-*l*60.5
Taylor expanded around -inf 0.0
Simplified0.0
if -5.395689961147953e+18 < x < 14413682.003027312Initial program 0.0
rmApplied div-inv0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2019121
(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))