\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 -11.259168060094455:\\
\;\;\;\;\left(\frac{0.5}{x} + \frac{0.15298196345929327}{\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)}\right) + \frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x}\\
\mathbf{elif}\;x \le 10.953272553467178:\\
\;\;\;\;\left(\left(\left(0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\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(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right)\right)\right) - \left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(2 \cdot 0.0001789971\right)\right) \cdot x\right) \cdot \log \left(e^{\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(\left(x \cdot x\right) \cdot 0.1049934947 + 1\right)\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.0005064034\right) + 0.0001789971 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right)}{\left(0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\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(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right)\right)\right) \cdot \left(0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\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(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right)\right)\right) - \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(2 \cdot 0.0001789971\right)\right) \cdot \left(\left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(2 \cdot 0.0001789971\right)\right)}}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{0.5}{x} + \frac{0.15298196345929327}{\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)}\right) + \frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x}\\
\end{array}double f(double x) {
double r5551332 = 1.0;
double r5551333 = 0.1049934947;
double r5551334 = x;
double r5551335 = r5551334 * r5551334;
double r5551336 = r5551333 * r5551335;
double r5551337 = r5551332 + r5551336;
double r5551338 = 0.0424060604;
double r5551339 = r5551335 * r5551335;
double r5551340 = r5551338 * r5551339;
double r5551341 = r5551337 + r5551340;
double r5551342 = 0.0072644182;
double r5551343 = r5551339 * r5551335;
double r5551344 = r5551342 * r5551343;
double r5551345 = r5551341 + r5551344;
double r5551346 = 0.0005064034;
double r5551347 = r5551343 * r5551335;
double r5551348 = r5551346 * r5551347;
double r5551349 = r5551345 + r5551348;
double r5551350 = 0.0001789971;
double r5551351 = r5551347 * r5551335;
double r5551352 = r5551350 * r5551351;
double r5551353 = r5551349 + r5551352;
double r5551354 = 0.7715471019;
double r5551355 = r5551354 * r5551335;
double r5551356 = r5551332 + r5551355;
double r5551357 = 0.2909738639;
double r5551358 = r5551357 * r5551339;
double r5551359 = r5551356 + r5551358;
double r5551360 = 0.0694555761;
double r5551361 = r5551360 * r5551343;
double r5551362 = r5551359 + r5551361;
double r5551363 = 0.0140005442;
double r5551364 = r5551363 * r5551347;
double r5551365 = r5551362 + r5551364;
double r5551366 = 0.0008327945;
double r5551367 = r5551366 * r5551351;
double r5551368 = r5551365 + r5551367;
double r5551369 = 2.0;
double r5551370 = r5551369 * r5551350;
double r5551371 = r5551351 * r5551335;
double r5551372 = r5551370 * r5551371;
double r5551373 = r5551368 + r5551372;
double r5551374 = r5551353 / r5551373;
double r5551375 = r5551374 * r5551334;
return r5551375;
}
double f(double x) {
double r5551376 = x;
double r5551377 = -11.259168060094455;
bool r5551378 = r5551376 <= r5551377;
double r5551379 = 0.5;
double r5551380 = r5551379 / r5551376;
double r5551381 = 0.15298196345929327;
double r5551382 = r5551376 * r5551376;
double r5551383 = r5551382 * r5551376;
double r5551384 = r5551382 * r5551383;
double r5551385 = r5551381 / r5551384;
double r5551386 = r5551380 + r5551385;
double r5551387 = 0.2514179000665375;
double r5551388 = r5551387 / r5551383;
double r5551389 = r5551386 + r5551388;
double r5551390 = 10.953272553467178;
bool r5551391 = r5551376 <= r5551390;
double r5551392 = 0.0008327945;
double r5551393 = r5551382 * r5551382;
double r5551394 = r5551393 * r5551382;
double r5551395 = r5551382 * r5551394;
double r5551396 = r5551395 * r5551382;
double r5551397 = r5551392 * r5551396;
double r5551398 = 0.0140005442;
double r5551399 = r5551395 * r5551398;
double r5551400 = 0.0694555761;
double r5551401 = r5551394 * r5551400;
double r5551402 = 0.2909738639;
double r5551403 = r5551402 * r5551393;
double r5551404 = 0.7715471019;
double r5551405 = r5551382 * r5551404;
double r5551406 = 1.0;
double r5551407 = r5551405 + r5551406;
double r5551408 = r5551403 + r5551407;
double r5551409 = r5551401 + r5551408;
double r5551410 = r5551399 + r5551409;
double r5551411 = r5551397 + r5551410;
double r5551412 = r5551396 * r5551382;
double r5551413 = 2.0;
double r5551414 = 0.0001789971;
double r5551415 = r5551413 * r5551414;
double r5551416 = r5551412 * r5551415;
double r5551417 = r5551411 - r5551416;
double r5551418 = r5551417 * r5551376;
double r5551419 = 0.0072644182;
double r5551420 = r5551419 * r5551394;
double r5551421 = 0.0424060604;
double r5551422 = r5551421 * r5551393;
double r5551423 = 0.1049934947;
double r5551424 = r5551382 * r5551423;
double r5551425 = r5551424 + r5551406;
double r5551426 = r5551422 + r5551425;
double r5551427 = r5551420 + r5551426;
double r5551428 = 0.0005064034;
double r5551429 = r5551395 * r5551428;
double r5551430 = r5551427 + r5551429;
double r5551431 = r5551414 * r5551396;
double r5551432 = r5551430 + r5551431;
double r5551433 = r5551411 * r5551411;
double r5551434 = r5551416 * r5551416;
double r5551435 = r5551433 - r5551434;
double r5551436 = r5551432 / r5551435;
double r5551437 = exp(r5551436);
double r5551438 = log(r5551437);
double r5551439 = r5551418 * r5551438;
double r5551440 = r5551391 ? r5551439 : r5551389;
double r5551441 = r5551378 ? r5551389 : r5551440;
return r5551441;
}



Bits error versus x
Results
if x < -11.259168060094455 or 10.953272553467178 < x Initial program 57.8
Taylor expanded around inf 0.2
Simplified0.2
if -11.259168060094455 < x < 10.953272553467178Initial program 0.0
rmApplied flip-+0.0
Applied associate-/r/0.0
Applied associate-*l*0.0
rmApplied add-log-exp0.1
Final simplification0.1
herbie shell --seed 2019129
(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))