\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 -1972.5019371328153:\\
\;\;\;\;\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 1119.150968415439:\\
\;\;\;\;\left(x \cdot \frac{\left(0.0005064034 \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(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) + 0.0072644182 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) + 0.0001789971 \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)}{\sqrt{\left(\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 + \left(\left(0.0694555761 \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.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)\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.0140005442\right)\right) + \left(0.0001789971 \cdot 2\right) \cdot \left(\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 \left(x \cdot x\right)\right)}}\right) \cdot \frac{1}{\sqrt{\left(\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 + \left(\left(0.0694555761 \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.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)\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.0140005442\right)\right) + \left(0.0001789971 \cdot 2\right) \cdot \left(\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 \left(x \cdot x\right)\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 r3596254 = 1.0;
double r3596255 = 0.1049934947;
double r3596256 = x;
double r3596257 = r3596256 * r3596256;
double r3596258 = r3596255 * r3596257;
double r3596259 = r3596254 + r3596258;
double r3596260 = 0.0424060604;
double r3596261 = r3596257 * r3596257;
double r3596262 = r3596260 * r3596261;
double r3596263 = r3596259 + r3596262;
double r3596264 = 0.0072644182;
double r3596265 = r3596261 * r3596257;
double r3596266 = r3596264 * r3596265;
double r3596267 = r3596263 + r3596266;
double r3596268 = 0.0005064034;
double r3596269 = r3596265 * r3596257;
double r3596270 = r3596268 * r3596269;
double r3596271 = r3596267 + r3596270;
double r3596272 = 0.0001789971;
double r3596273 = r3596269 * r3596257;
double r3596274 = r3596272 * r3596273;
double r3596275 = r3596271 + r3596274;
double r3596276 = 0.7715471019;
double r3596277 = r3596276 * r3596257;
double r3596278 = r3596254 + r3596277;
double r3596279 = 0.2909738639;
double r3596280 = r3596279 * r3596261;
double r3596281 = r3596278 + r3596280;
double r3596282 = 0.0694555761;
double r3596283 = r3596282 * r3596265;
double r3596284 = r3596281 + r3596283;
double r3596285 = 0.0140005442;
double r3596286 = r3596285 * r3596269;
double r3596287 = r3596284 + r3596286;
double r3596288 = 0.0008327945;
double r3596289 = r3596288 * r3596273;
double r3596290 = r3596287 + r3596289;
double r3596291 = 2.0;
double r3596292 = r3596291 * r3596272;
double r3596293 = r3596273 * r3596257;
double r3596294 = r3596292 * r3596293;
double r3596295 = r3596290 + r3596294;
double r3596296 = r3596275 / r3596295;
double r3596297 = r3596296 * r3596256;
return r3596297;
}
double f(double x) {
double r3596298 = x;
double r3596299 = -1972.5019371328153;
bool r3596300 = r3596298 <= r3596299;
double r3596301 = 0.15298196345929327;
double r3596302 = 5.0;
double r3596303 = pow(r3596298, r3596302);
double r3596304 = r3596301 / r3596303;
double r3596305 = 0.2514179000665375;
double r3596306 = r3596298 * r3596298;
double r3596307 = r3596298 * r3596306;
double r3596308 = r3596305 / r3596307;
double r3596309 = 0.5;
double r3596310 = r3596309 / r3596298;
double r3596311 = r3596308 + r3596310;
double r3596312 = r3596304 + r3596311;
double r3596313 = 1119.150968415439;
bool r3596314 = r3596298 <= r3596313;
double r3596315 = 0.0005064034;
double r3596316 = r3596306 * r3596306;
double r3596317 = r3596306 * r3596316;
double r3596318 = r3596306 * r3596317;
double r3596319 = r3596315 * r3596318;
double r3596320 = 0.0424060604;
double r3596321 = r3596320 * r3596316;
double r3596322 = 0.1049934947;
double r3596323 = r3596322 * r3596306;
double r3596324 = 1.0;
double r3596325 = r3596323 + r3596324;
double r3596326 = r3596321 + r3596325;
double r3596327 = 0.0072644182;
double r3596328 = r3596327 * r3596317;
double r3596329 = r3596326 + r3596328;
double r3596330 = r3596319 + r3596329;
double r3596331 = 0.0001789971;
double r3596332 = r3596306 * r3596318;
double r3596333 = r3596331 * r3596332;
double r3596334 = r3596330 + r3596333;
double r3596335 = 0.0008327945;
double r3596336 = r3596332 * r3596335;
double r3596337 = 0.0694555761;
double r3596338 = r3596337 * r3596317;
double r3596339 = 0.7715471019;
double r3596340 = r3596339 * r3596306;
double r3596341 = r3596324 + r3596340;
double r3596342 = 0.2909738639;
double r3596343 = r3596342 * r3596316;
double r3596344 = r3596341 + r3596343;
double r3596345 = r3596338 + r3596344;
double r3596346 = 0.0140005442;
double r3596347 = r3596318 * r3596346;
double r3596348 = r3596345 + r3596347;
double r3596349 = r3596336 + r3596348;
double r3596350 = 2.0;
double r3596351 = r3596331 * r3596350;
double r3596352 = r3596332 * r3596306;
double r3596353 = r3596351 * r3596352;
double r3596354 = r3596349 + r3596353;
double r3596355 = sqrt(r3596354);
double r3596356 = r3596334 / r3596355;
double r3596357 = r3596298 * r3596356;
double r3596358 = r3596324 / r3596355;
double r3596359 = r3596357 * r3596358;
double r3596360 = r3596314 ? r3596359 : r3596312;
double r3596361 = r3596300 ? r3596312 : r3596360;
return r3596361;
}



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