\frac{\left(\left(\left(\left(1 + 0.1049934946999999951788851149103720672429 \cdot \left(x \cdot x\right)\right) + 0.04240606040000000076517494562722276896238 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.007264418199999999985194687468492702464573 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 5.064034000000000243502107366566633572802 \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.789971000000000009994005623070734145585 \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.7715471018999999763821051601553335785866 \cdot \left(x \cdot x\right)\right) + 0.2909738639000000182122107617033179849386 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.06945557609999999937322456844412954524159 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.01400054419999999938406531896362139377743 \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.327945000000000442749725770852364803432 \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.789971000000000009994005623070734145585 \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 -246126550036688416 \lor \neg \left(x \le 3959.752675162330433522583916783332824707\right):\\
\;\;\;\;\frac{0.5}{x} + \left(\frac{0.2514179000665375252054900556686334311962}{{x}^{3}} + \frac{0.1529819634592932686700805788859724998474}{{x}^{5}}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{\left(\left(\left(\left(1 + 0.1049934946999999951788851149103720672429 \cdot \left(x \cdot x\right)\right) + 0.04240606040000000076517494562722276896238 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.007264418199999999985194687468492702464573 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 5.064034000000000243502107366566633572802 \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.789971000000000009994005623070734145585 \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.7715471018999999763821051601553335785866 \cdot \left(x \cdot x\right)\right) + 0.2909738639000000182122107617033179849386 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + 0.06945557609999999937322456844412954524159 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) + 0.01400054419999999938406531896362139377743 \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.327945000000000442749725770852364803432 \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.789971000000000009994005623070734145585 \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 \left(\sqrt{\frac{5.064034000000000243502107366566633572802 \cdot 10^{-4} \cdot {x}^{8} + \left(1.789971000000000009994005623070734145585 \cdot 10^{-4} \cdot {x}^{10} + \left(0.007264418199999999985194687468492702464573 \cdot {x}^{6} + \left({x}^{4} \cdot 0.04240606040000000076517494562722276896238 + \left(1 + 0.1049934946999999951788851149103720672429 \cdot \left(x \cdot x\right)\right)\right)\right)\right)}{\left({\left({x}^{4}\right)}^{3} \cdot \left(2 \cdot 1.789971000000000009994005623070734145585 \cdot 10^{-4}\right) + \left(0.06945557609999999937322456844412954524159 \cdot {x}^{6} + \left(\left(1 + 0.7715471018999999763821051601553335785866 \cdot \left(x \cdot x\right)\right) + {x}^{4} \cdot 0.2909738639000000182122107617033179849386\right)\right)\right) + \left(8.327945000000000442749725770852364803432 \cdot 10^{-4} \cdot {x}^{10} + {x}^{8} \cdot 0.01400054419999999938406531896362139377743\right)}} \cdot x\right)\\
\end{array}double f(double x) {
double r166330 = 1.0;
double r166331 = 0.1049934947;
double r166332 = x;
double r166333 = r166332 * r166332;
double r166334 = r166331 * r166333;
double r166335 = r166330 + r166334;
double r166336 = 0.0424060604;
double r166337 = r166333 * r166333;
double r166338 = r166336 * r166337;
double r166339 = r166335 + r166338;
double r166340 = 0.0072644182;
double r166341 = r166337 * r166333;
double r166342 = r166340 * r166341;
double r166343 = r166339 + r166342;
double r166344 = 0.0005064034;
double r166345 = r166341 * r166333;
double r166346 = r166344 * r166345;
double r166347 = r166343 + r166346;
double r166348 = 0.0001789971;
double r166349 = r166345 * r166333;
double r166350 = r166348 * r166349;
double r166351 = r166347 + r166350;
double r166352 = 0.7715471019;
double r166353 = r166352 * r166333;
double r166354 = r166330 + r166353;
double r166355 = 0.2909738639;
double r166356 = r166355 * r166337;
double r166357 = r166354 + r166356;
double r166358 = 0.0694555761;
double r166359 = r166358 * r166341;
double r166360 = r166357 + r166359;
double r166361 = 0.0140005442;
double r166362 = r166361 * r166345;
double r166363 = r166360 + r166362;
double r166364 = 0.0008327945;
double r166365 = r166364 * r166349;
double r166366 = r166363 + r166365;
double r166367 = 2.0;
double r166368 = r166367 * r166348;
double r166369 = r166349 * r166333;
double r166370 = r166368 * r166369;
double r166371 = r166366 + r166370;
double r166372 = r166351 / r166371;
double r166373 = r166372 * r166332;
return r166373;
}
double f(double x) {
double r166374 = x;
double r166375 = -2.4612655003668842e+17;
bool r166376 = r166374 <= r166375;
double r166377 = 3959.7526751623304;
bool r166378 = r166374 <= r166377;
double r166379 = !r166378;
bool r166380 = r166376 || r166379;
double r166381 = 0.5;
double r166382 = r166381 / r166374;
double r166383 = 0.2514179000665375;
double r166384 = 3.0;
double r166385 = pow(r166374, r166384);
double r166386 = r166383 / r166385;
double r166387 = 0.15298196345929327;
double r166388 = 5.0;
double r166389 = pow(r166374, r166388);
double r166390 = r166387 / r166389;
double r166391 = r166386 + r166390;
double r166392 = r166382 + r166391;
double r166393 = 1.0;
double r166394 = 0.1049934947;
double r166395 = r166374 * r166374;
double r166396 = r166394 * r166395;
double r166397 = r166393 + r166396;
double r166398 = 0.0424060604;
double r166399 = r166395 * r166395;
double r166400 = r166398 * r166399;
double r166401 = r166397 + r166400;
double r166402 = 0.0072644182;
double r166403 = r166399 * r166395;
double r166404 = r166402 * r166403;
double r166405 = r166401 + r166404;
double r166406 = 0.0005064034;
double r166407 = r166403 * r166395;
double r166408 = r166406 * r166407;
double r166409 = r166405 + r166408;
double r166410 = 0.0001789971;
double r166411 = r166407 * r166395;
double r166412 = r166410 * r166411;
double r166413 = r166409 + r166412;
double r166414 = 0.7715471019;
double r166415 = r166414 * r166395;
double r166416 = r166393 + r166415;
double r166417 = 0.2909738639;
double r166418 = r166417 * r166399;
double r166419 = r166416 + r166418;
double r166420 = 0.0694555761;
double r166421 = r166420 * r166403;
double r166422 = r166419 + r166421;
double r166423 = 0.0140005442;
double r166424 = r166423 * r166407;
double r166425 = r166422 + r166424;
double r166426 = 0.0008327945;
double r166427 = r166426 * r166411;
double r166428 = r166425 + r166427;
double r166429 = 2.0;
double r166430 = r166429 * r166410;
double r166431 = r166411 * r166395;
double r166432 = r166430 * r166431;
double r166433 = r166428 + r166432;
double r166434 = r166413 / r166433;
double r166435 = sqrt(r166434);
double r166436 = 8.0;
double r166437 = pow(r166374, r166436);
double r166438 = r166406 * r166437;
double r166439 = 10.0;
double r166440 = pow(r166374, r166439);
double r166441 = r166410 * r166440;
double r166442 = 6.0;
double r166443 = pow(r166374, r166442);
double r166444 = r166402 * r166443;
double r166445 = 4.0;
double r166446 = pow(r166374, r166445);
double r166447 = r166446 * r166398;
double r166448 = r166447 + r166397;
double r166449 = r166444 + r166448;
double r166450 = r166441 + r166449;
double r166451 = r166438 + r166450;
double r166452 = pow(r166446, r166384);
double r166453 = r166452 * r166430;
double r166454 = r166420 * r166443;
double r166455 = r166446 * r166417;
double r166456 = r166416 + r166455;
double r166457 = r166454 + r166456;
double r166458 = r166453 + r166457;
double r166459 = r166426 * r166440;
double r166460 = r166437 * r166423;
double r166461 = r166459 + r166460;
double r166462 = r166458 + r166461;
double r166463 = r166451 / r166462;
double r166464 = sqrt(r166463);
double r166465 = r166464 * r166374;
double r166466 = r166435 * r166465;
double r166467 = r166380 ? r166392 : r166466;
return r166467;
}



Bits error versus x
Results
if x < -2.4612655003668842e+17 or 3959.7526751623304 < x Initial program 60.5
Taylor expanded around inf 0.0
Simplified0.0
if -2.4612655003668842e+17 < x < 3959.7526751623304Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied associate-*l*0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019326
(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))