\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 -22331647.500376165:\\
\;\;\;\;\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 1595.5609010574753:\\
\;\;\;\;\left(x \cdot \frac{\sqrt{\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)}}{\left(\left(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) + \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(0.7715471019 \cdot \left(x \cdot x\right) + 1\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(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) \cdot 0.0008327945\right) + \left(0.0001789971 \cdot 2\right) \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(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right)\right)}\right) \cdot \sqrt{\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)}\\
\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 r3654125 = 1.0;
double r3654126 = 0.1049934947;
double r3654127 = x;
double r3654128 = r3654127 * r3654127;
double r3654129 = r3654126 * r3654128;
double r3654130 = r3654125 + r3654129;
double r3654131 = 0.0424060604;
double r3654132 = r3654128 * r3654128;
double r3654133 = r3654131 * r3654132;
double r3654134 = r3654130 + r3654133;
double r3654135 = 0.0072644182;
double r3654136 = r3654132 * r3654128;
double r3654137 = r3654135 * r3654136;
double r3654138 = r3654134 + r3654137;
double r3654139 = 0.0005064034;
double r3654140 = r3654136 * r3654128;
double r3654141 = r3654139 * r3654140;
double r3654142 = r3654138 + r3654141;
double r3654143 = 0.0001789971;
double r3654144 = r3654140 * r3654128;
double r3654145 = r3654143 * r3654144;
double r3654146 = r3654142 + r3654145;
double r3654147 = 0.7715471019;
double r3654148 = r3654147 * r3654128;
double r3654149 = r3654125 + r3654148;
double r3654150 = 0.2909738639;
double r3654151 = r3654150 * r3654132;
double r3654152 = r3654149 + r3654151;
double r3654153 = 0.0694555761;
double r3654154 = r3654153 * r3654136;
double r3654155 = r3654152 + r3654154;
double r3654156 = 0.0140005442;
double r3654157 = r3654156 * r3654140;
double r3654158 = r3654155 + r3654157;
double r3654159 = 0.0008327945;
double r3654160 = r3654159 * r3654144;
double r3654161 = r3654158 + r3654160;
double r3654162 = 2.0;
double r3654163 = r3654162 * r3654143;
double r3654164 = r3654144 * r3654128;
double r3654165 = r3654163 * r3654164;
double r3654166 = r3654161 + r3654165;
double r3654167 = r3654146 / r3654166;
double r3654168 = r3654167 * r3654127;
return r3654168;
}
double f(double x) {
double r3654169 = x;
double r3654170 = -22331647.500376165;
bool r3654171 = r3654169 <= r3654170;
double r3654172 = 0.15298196345929327;
double r3654173 = 5.0;
double r3654174 = pow(r3654169, r3654173);
double r3654175 = r3654172 / r3654174;
double r3654176 = 0.2514179000665375;
double r3654177 = r3654169 * r3654169;
double r3654178 = r3654169 * r3654177;
double r3654179 = r3654176 / r3654178;
double r3654180 = 0.5;
double r3654181 = r3654180 / r3654169;
double r3654182 = r3654179 + r3654181;
double r3654183 = r3654175 + r3654182;
double r3654184 = 1595.5609010574753;
bool r3654185 = r3654169 <= r3654184;
double r3654186 = r3654177 * r3654177;
double r3654187 = r3654177 * r3654186;
double r3654188 = r3654177 * r3654187;
double r3654189 = r3654177 * r3654188;
double r3654190 = 0.0001789971;
double r3654191 = r3654189 * r3654190;
double r3654192 = 0.0072644182;
double r3654193 = r3654192 * r3654187;
double r3654194 = 1.0;
double r3654195 = 0.1049934947;
double r3654196 = r3654195 * r3654177;
double r3654197 = r3654194 + r3654196;
double r3654198 = 0.0424060604;
double r3654199 = r3654198 * r3654186;
double r3654200 = r3654197 + r3654199;
double r3654201 = r3654193 + r3654200;
double r3654202 = 0.0005064034;
double r3654203 = r3654188 * r3654202;
double r3654204 = r3654201 + r3654203;
double r3654205 = r3654191 + r3654204;
double r3654206 = sqrt(r3654205);
double r3654207 = 0.0140005442;
double r3654208 = r3654207 * r3654188;
double r3654209 = 0.0694555761;
double r3654210 = r3654187 * r3654209;
double r3654211 = 0.2909738639;
double r3654212 = r3654186 * r3654211;
double r3654213 = 0.7715471019;
double r3654214 = r3654213 * r3654177;
double r3654215 = r3654214 + r3654194;
double r3654216 = r3654212 + r3654215;
double r3654217 = r3654210 + r3654216;
double r3654218 = r3654208 + r3654217;
double r3654219 = 0.0008327945;
double r3654220 = r3654189 * r3654219;
double r3654221 = r3654218 + r3654220;
double r3654222 = 2.0;
double r3654223 = r3654190 * r3654222;
double r3654224 = r3654177 * r3654189;
double r3654225 = r3654223 * r3654224;
double r3654226 = r3654221 + r3654225;
double r3654227 = r3654206 / r3654226;
double r3654228 = r3654169 * r3654227;
double r3654229 = r3654228 * r3654206;
double r3654230 = r3654185 ? r3654229 : r3654183;
double r3654231 = r3654171 ? r3654183 : r3654230;
return r3654231;
}



Bits error versus x
Results
if x < -22331647.500376165 or 1595.5609010574753 < x Initial program 59.1
Taylor expanded around inf 0.0
Simplified0.0
if -22331647.500376165 < x < 1595.5609010574753Initial program 0.0
rmApplied *-un-lft-identity0.0
Applied add-sqr-sqrt0.0
Applied times-frac0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2019154
(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))