\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 -1356525439.54039 \lor \neg \left(x \le 789.84959973016066\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}:\\
\;\;\;\;\left(\left(\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)} \cdot \sqrt{\frac{1}{\left(\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 {x}^{2}\right)\right) + \left(0.014000544199999999 \cdot {x}^{8} + {x}^{2} \cdot \left(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) + \left(2 \cdot 1.789971 \cdot 10^{-4}\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)\right)\right) \cdot 1}}\right) \cdot \left(\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)} \cdot \sqrt{\frac{1}{\left(\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 {x}^{2}\right)\right) + \left(0.014000544199999999 \cdot {x}^{8} + {x}^{2} \cdot \left(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) + \left(2 \cdot 1.789971 \cdot 10^{-4}\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)\right)\right) \cdot 1}}\right)\right) \cdot x\\
\end{array}double f(double x) {
double r236463 = 1.0;
double r236464 = 0.1049934947;
double r236465 = x;
double r236466 = r236465 * r236465;
double r236467 = r236464 * r236466;
double r236468 = r236463 + r236467;
double r236469 = 0.0424060604;
double r236470 = r236466 * r236466;
double r236471 = r236469 * r236470;
double r236472 = r236468 + r236471;
double r236473 = 0.0072644182;
double r236474 = r236470 * r236466;
double r236475 = r236473 * r236474;
double r236476 = r236472 + r236475;
double r236477 = 0.0005064034;
double r236478 = r236474 * r236466;
double r236479 = r236477 * r236478;
double r236480 = r236476 + r236479;
double r236481 = 0.0001789971;
double r236482 = r236478 * r236466;
double r236483 = r236481 * r236482;
double r236484 = r236480 + r236483;
double r236485 = 0.7715471019;
double r236486 = r236485 * r236466;
double r236487 = r236463 + r236486;
double r236488 = 0.2909738639;
double r236489 = r236488 * r236470;
double r236490 = r236487 + r236489;
double r236491 = 0.0694555761;
double r236492 = r236491 * r236474;
double r236493 = r236490 + r236492;
double r236494 = 0.0140005442;
double r236495 = r236494 * r236478;
double r236496 = r236493 + r236495;
double r236497 = 0.0008327945;
double r236498 = r236497 * r236482;
double r236499 = r236496 + r236498;
double r236500 = 2.0;
double r236501 = r236500 * r236481;
double r236502 = r236482 * r236466;
double r236503 = r236501 * r236502;
double r236504 = r236499 + r236503;
double r236505 = r236484 / r236504;
double r236506 = r236505 * r236465;
return r236506;
}
double f(double x) {
double r236507 = x;
double r236508 = -1356525439.54039;
bool r236509 = r236507 <= r236508;
double r236510 = 789.8495997301607;
bool r236511 = r236507 <= r236510;
double r236512 = !r236511;
bool r236513 = r236509 || r236512;
double r236514 = 0.2514179000665375;
double r236515 = 1.0;
double r236516 = 3.0;
double r236517 = pow(r236507, r236516);
double r236518 = r236515 / r236517;
double r236519 = r236514 * r236518;
double r236520 = 0.15298196345929327;
double r236521 = 5.0;
double r236522 = pow(r236507, r236521);
double r236523 = r236515 / r236522;
double r236524 = r236520 * r236523;
double r236525 = 0.5;
double r236526 = r236515 / r236507;
double r236527 = r236525 * r236526;
double r236528 = r236524 + r236527;
double r236529 = r236519 + r236528;
double r236530 = 1.0;
double r236531 = 0.1049934947;
double r236532 = r236507 * r236507;
double r236533 = r236531 * r236532;
double r236534 = r236530 + r236533;
double r236535 = 0.0424060604;
double r236536 = r236532 * r236532;
double r236537 = r236535 * r236536;
double r236538 = r236534 + r236537;
double r236539 = 0.0072644182;
double r236540 = r236536 * r236532;
double r236541 = r236539 * r236540;
double r236542 = r236538 + r236541;
double r236543 = 0.0005064034;
double r236544 = r236540 * r236532;
double r236545 = r236543 * r236544;
double r236546 = r236542 + r236545;
double r236547 = 0.0001789971;
double r236548 = r236544 * r236532;
double r236549 = r236547 * r236548;
double r236550 = r236546 + r236549;
double r236551 = sqrt(r236550);
double r236552 = 2.0;
double r236553 = pow(r236507, r236552);
double r236554 = r236507 * r236517;
double r236555 = r236553 * r236554;
double r236556 = 0.0694555761;
double r236557 = r236555 * r236556;
double r236558 = r236557 + r236530;
double r236559 = 0.7715471019;
double r236560 = 0.2909738639;
double r236561 = r236560 * r236553;
double r236562 = r236559 + r236561;
double r236563 = r236553 * r236562;
double r236564 = r236558 + r236563;
double r236565 = 0.0140005442;
double r236566 = 8.0;
double r236567 = pow(r236507, r236566);
double r236568 = r236565 * r236567;
double r236569 = 0.0008327945;
double r236570 = r236569 * r236544;
double r236571 = 2.0;
double r236572 = r236571 * r236547;
double r236573 = r236572 * r236548;
double r236574 = r236570 + r236573;
double r236575 = r236553 * r236574;
double r236576 = r236568 + r236575;
double r236577 = r236564 + r236576;
double r236578 = r236577 * r236515;
double r236579 = r236515 / r236578;
double r236580 = sqrt(r236579);
double r236581 = r236551 * r236580;
double r236582 = r236581 * r236581;
double r236583 = r236582 * r236507;
double r236584 = r236513 ? r236529 : r236583;
return r236584;
}



Bits error versus x
Results
if x < -1356525439.54039 or 789.8495997301607 < x Initial program 59.6
Taylor expanded around inf 0.0
if -1356525439.54039 < x < 789.8495997301607Initial program 0.0
rmApplied div-inv0.0
Simplified0.0
rmApplied add-sqr-sqrt0.0
Applied add-sqr-sqrt0.0
Applied unswap-sqr0.0
Final simplification0.0
herbie shell --seed 2020056
(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))