\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 -228768963.19517976:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\mathbf{elif}\;x \le 23326.04538430404:\\
\;\;\;\;x \cdot \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(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) + \sqrt{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)} \cdot \sqrt{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(2 \cdot 0.0001789971\right) \cdot \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) + \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 0.0008327945 + \left(\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) + 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)\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 r42111333 = 1.0;
double r42111334 = 0.1049934947;
double r42111335 = x;
double r42111336 = r42111335 * r42111335;
double r42111337 = r42111334 * r42111336;
double r42111338 = r42111333 + r42111337;
double r42111339 = 0.0424060604;
double r42111340 = r42111336 * r42111336;
double r42111341 = r42111339 * r42111340;
double r42111342 = r42111338 + r42111341;
double r42111343 = 0.0072644182;
double r42111344 = r42111340 * r42111336;
double r42111345 = r42111343 * r42111344;
double r42111346 = r42111342 + r42111345;
double r42111347 = 0.0005064034;
double r42111348 = r42111344 * r42111336;
double r42111349 = r42111347 * r42111348;
double r42111350 = r42111346 + r42111349;
double r42111351 = 0.0001789971;
double r42111352 = r42111348 * r42111336;
double r42111353 = r42111351 * r42111352;
double r42111354 = r42111350 + r42111353;
double r42111355 = 0.7715471019;
double r42111356 = r42111355 * r42111336;
double r42111357 = r42111333 + r42111356;
double r42111358 = 0.2909738639;
double r42111359 = r42111358 * r42111340;
double r42111360 = r42111357 + r42111359;
double r42111361 = 0.0694555761;
double r42111362 = r42111361 * r42111344;
double r42111363 = r42111360 + r42111362;
double r42111364 = 0.0140005442;
double r42111365 = r42111364 * r42111348;
double r42111366 = r42111363 + r42111365;
double r42111367 = 0.0008327945;
double r42111368 = r42111367 * r42111352;
double r42111369 = r42111366 + r42111368;
double r42111370 = 2.0;
double r42111371 = r42111370 * r42111351;
double r42111372 = r42111352 * r42111336;
double r42111373 = r42111371 * r42111372;
double r42111374 = r42111369 + r42111373;
double r42111375 = r42111354 / r42111374;
double r42111376 = r42111375 * r42111335;
return r42111376;
}
double f(double x) {
double r42111377 = x;
double r42111378 = -228768963.19517976;
bool r42111379 = r42111377 <= r42111378;
double r42111380 = 0.5;
double r42111381 = 0.2514179000665375;
double r42111382 = r42111377 * r42111377;
double r42111383 = r42111381 / r42111382;
double r42111384 = r42111380 + r42111383;
double r42111385 = r42111384 / r42111377;
double r42111386 = 0.15298196345929327;
double r42111387 = 5.0;
double r42111388 = pow(r42111377, r42111387);
double r42111389 = r42111386 / r42111388;
double r42111390 = r42111385 + r42111389;
double r42111391 = 23326.04538430404;
bool r42111392 = r42111377 <= r42111391;
double r42111393 = 0.0072644182;
double r42111394 = r42111382 * r42111382;
double r42111395 = r42111394 * r42111382;
double r42111396 = r42111393 * r42111395;
double r42111397 = 0.0424060604;
double r42111398 = r42111397 * r42111394;
double r42111399 = 0.1049934947;
double r42111400 = r42111399 * r42111382;
double r42111401 = 1.0;
double r42111402 = r42111400 + r42111401;
double r42111403 = r42111398 + r42111402;
double r42111404 = r42111396 + r42111403;
double r42111405 = r42111395 * r42111382;
double r42111406 = 0.0005064034;
double r42111407 = r42111405 * r42111406;
double r42111408 = r42111404 + r42111407;
double r42111409 = 0.0001789971;
double r42111410 = r42111405 * r42111382;
double r42111411 = r42111409 * r42111410;
double r42111412 = sqrt(r42111411);
double r42111413 = r42111412 * r42111412;
double r42111414 = r42111408 + r42111413;
double r42111415 = 2.0;
double r42111416 = r42111415 * r42111409;
double r42111417 = r42111382 * r42111410;
double r42111418 = r42111416 * r42111417;
double r42111419 = 0.0008327945;
double r42111420 = r42111410 * r42111419;
double r42111421 = 0.0694555761;
double r42111422 = r42111395 * r42111421;
double r42111423 = 0.2909738639;
double r42111424 = r42111423 * r42111394;
double r42111425 = 0.7715471019;
double r42111426 = r42111382 * r42111425;
double r42111427 = r42111426 + r42111401;
double r42111428 = r42111424 + r42111427;
double r42111429 = r42111422 + r42111428;
double r42111430 = 0.0140005442;
double r42111431 = r42111430 * r42111405;
double r42111432 = r42111429 + r42111431;
double r42111433 = r42111420 + r42111432;
double r42111434 = r42111418 + r42111433;
double r42111435 = r42111414 / r42111434;
double r42111436 = r42111377 * r42111435;
double r42111437 = r42111392 ? r42111436 : r42111390;
double r42111438 = r42111379 ? r42111390 : r42111437;
return r42111438;
}



Bits error versus x
Results
if x < -228768963.19517976 or 23326.04538430404 < x Initial program 59.3
Taylor expanded around inf 0.0
Simplified0.0
if -228768963.19517976 < x < 23326.04538430404Initial program 0.0
rmApplied add-sqr-sqrt0.0
Final simplification0.0
herbie shell --seed 2019125
(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))