\frac{\left(\left(\left(\left(1 + 0.1049934947 \cdot \left(x \cdot x\right)\right) + 0.042406060400000001 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.00726441819999999999 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 5.0640340000000002 \cdot 10^{-4} \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) + 1.789971 \cdot 10^{-4} \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.77154710189999998 \cdot \left(x \cdot x\right)\right) + 0.29097386390000002 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.069455576099999999 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.014000544199999999 \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) + 8.32794500000000044 \cdot 10^{-4} \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 1.789971 \cdot 10^{-4}\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 -6753.5048671236073 \lor \neg \left(x \le 7293.4466097268923\right):\\
\;\;\;\;0.25141790006653753 \cdot \frac{1}{{x}^{3}} + \left(0.1529819634592933 \cdot \frac{1}{{x}^{5}} + 0.5 \cdot \frac{1}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\left(\left(\left(\left(1 + 0.1049934947 \cdot \left(x \cdot x\right)\right) + 0.042406060400000001 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.00726441819999999999 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 5.0640340000000002 \cdot 10^{-4} \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) + 1.789971 \cdot 10^{-4} \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)}}{1} \cdot \frac{x \cdot \sqrt{\left(\left(\left(\left(1 + 0.1049934947 \cdot \left(x \cdot x\right)\right) + 0.042406060400000001 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.00726441819999999999 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 5.0640340000000002 \cdot 10^{-4} \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) + 1.789971 \cdot 10^{-4} \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({x}^{2} \cdot \left({x}^{2} \cdot \left({x}^{2} \cdot \left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right)\right)\right)\right) \cdot \left(2 \cdot 1.789971 \cdot 10^{-4}\right) + \left(\left(\left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right) \cdot 0.069455576099999999 + 1\right) + {x}^{2} \cdot \left(0.77154710189999998 + 0.29097386390000002 \cdot \left(x \cdot x\right)\right)\right)\right) + {x}^{2} \cdot \left(0.014000544199999999 \cdot \left(\left(\left(x \cdot x\right) \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)\right) + 8.32794500000000044 \cdot 10^{-4} \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)}\\
\end{array}double f(double x) {
double r254388 = 1.0;
double r254389 = 0.1049934947;
double r254390 = x;
double r254391 = r254390 * r254390;
double r254392 = r254389 * r254391;
double r254393 = r254388 + r254392;
double r254394 = 0.0424060604;
double r254395 = r254391 * r254391;
double r254396 = r254394 * r254395;
double r254397 = r254393 + r254396;
double r254398 = 0.0072644182;
double r254399 = r254395 * r254391;
double r254400 = r254398 * r254399;
double r254401 = r254397 + r254400;
double r254402 = 0.0005064034;
double r254403 = r254399 * r254391;
double r254404 = r254402 * r254403;
double r254405 = r254401 + r254404;
double r254406 = 0.0001789971;
double r254407 = r254403 * r254391;
double r254408 = r254406 * r254407;
double r254409 = r254405 + r254408;
double r254410 = 0.7715471019;
double r254411 = r254410 * r254391;
double r254412 = r254388 + r254411;
double r254413 = 0.2909738639;
double r254414 = r254413 * r254395;
double r254415 = r254412 + r254414;
double r254416 = 0.0694555761;
double r254417 = r254416 * r254399;
double r254418 = r254415 + r254417;
double r254419 = 0.0140005442;
double r254420 = r254419 * r254403;
double r254421 = r254418 + r254420;
double r254422 = 0.0008327945;
double r254423 = r254422 * r254407;
double r254424 = r254421 + r254423;
double r254425 = 2.0;
double r254426 = r254425 * r254406;
double r254427 = r254407 * r254391;
double r254428 = r254426 * r254427;
double r254429 = r254424 + r254428;
double r254430 = r254409 / r254429;
double r254431 = r254430 * r254390;
return r254431;
}
double f(double x) {
double r254432 = x;
double r254433 = -6753.504867123607;
bool r254434 = r254432 <= r254433;
double r254435 = 7293.446609726892;
bool r254436 = r254432 <= r254435;
double r254437 = !r254436;
bool r254438 = r254434 || r254437;
double r254439 = 0.2514179000665375;
double r254440 = 1.0;
double r254441 = 3.0;
double r254442 = pow(r254432, r254441);
double r254443 = r254440 / r254442;
double r254444 = r254439 * r254443;
double r254445 = 0.15298196345929327;
double r254446 = 5.0;
double r254447 = pow(r254432, r254446);
double r254448 = r254440 / r254447;
double r254449 = r254445 * r254448;
double r254450 = 0.5;
double r254451 = r254440 / r254432;
double r254452 = r254450 * r254451;
double r254453 = r254449 + r254452;
double r254454 = r254444 + r254453;
double r254455 = 1.0;
double r254456 = 0.1049934947;
double r254457 = r254432 * r254432;
double r254458 = r254456 * r254457;
double r254459 = r254455 + r254458;
double r254460 = 0.0424060604;
double r254461 = r254457 * r254457;
double r254462 = r254460 * r254461;
double r254463 = r254459 + r254462;
double r254464 = 0.0072644182;
double r254465 = r254461 * r254457;
double r254466 = r254464 * r254465;
double r254467 = r254463 + r254466;
double r254468 = 0.0005064034;
double r254469 = r254465 * r254457;
double r254470 = r254468 * r254469;
double r254471 = r254467 + r254470;
double r254472 = 0.0001789971;
double r254473 = r254469 * r254457;
double r254474 = r254472 * r254473;
double r254475 = r254471 + r254474;
double r254476 = sqrt(r254475);
double r254477 = r254476 / r254440;
double r254478 = r254432 * r254476;
double r254479 = 2.0;
double r254480 = pow(r254432, r254479);
double r254481 = r254432 * r254442;
double r254482 = r254480 * r254481;
double r254483 = r254480 * r254482;
double r254484 = r254480 * r254483;
double r254485 = r254480 * r254484;
double r254486 = 2.0;
double r254487 = r254486 * r254472;
double r254488 = r254485 * r254487;
double r254489 = 0.0694555761;
double r254490 = r254482 * r254489;
double r254491 = r254490 + r254455;
double r254492 = 0.7715471019;
double r254493 = 0.2909738639;
double r254494 = r254493 * r254457;
double r254495 = r254492 + r254494;
double r254496 = r254480 * r254495;
double r254497 = r254491 + r254496;
double r254498 = r254488 + r254497;
double r254499 = 0.0140005442;
double r254500 = r254457 * r254432;
double r254501 = r254500 * r254500;
double r254502 = r254499 * r254501;
double r254503 = 0.0008327945;
double r254504 = r254503 * r254469;
double r254505 = r254502 + r254504;
double r254506 = r254480 * r254505;
double r254507 = r254498 + r254506;
double r254508 = r254478 / r254507;
double r254509 = r254477 * r254508;
double r254510 = r254438 ? r254454 : r254509;
return r254510;
}



Bits error versus x
Results
if x < -6753.504867123607 or 7293.446609726892 < x Initial program 59.3
Taylor expanded around inf 0.0
if -6753.504867123607 < x < 7293.446609726892Initial program 0.0
rmApplied *-un-lft-identity0.0
Applied add-sqr-sqrt0.0
Applied times-frac0.0
Applied associate-*l*0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2020062
(FPCore (x)
:name "Jmat.Real.dawson"
:precision binary64
(* (/ (+ (+ (+ (+ (+ 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))