\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 -14623663576.376369 \lor \neg \left(x \le 702.997215207723116\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(\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)\right) \cdot \left(\frac{1}{\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\right)\\
\end{array}double f(double x) {
double r242279 = 1.0;
double r242280 = 0.1049934947;
double r242281 = x;
double r242282 = r242281 * r242281;
double r242283 = r242280 * r242282;
double r242284 = r242279 + r242283;
double r242285 = 0.0424060604;
double r242286 = r242282 * r242282;
double r242287 = r242285 * r242286;
double r242288 = r242284 + r242287;
double r242289 = 0.0072644182;
double r242290 = r242286 * r242282;
double r242291 = r242289 * r242290;
double r242292 = r242288 + r242291;
double r242293 = 0.0005064034;
double r242294 = r242290 * r242282;
double r242295 = r242293 * r242294;
double r242296 = r242292 + r242295;
double r242297 = 0.0001789971;
double r242298 = r242294 * r242282;
double r242299 = r242297 * r242298;
double r242300 = r242296 + r242299;
double r242301 = 0.7715471019;
double r242302 = r242301 * r242282;
double r242303 = r242279 + r242302;
double r242304 = 0.2909738639;
double r242305 = r242304 * r242286;
double r242306 = r242303 + r242305;
double r242307 = 0.0694555761;
double r242308 = r242307 * r242290;
double r242309 = r242306 + r242308;
double r242310 = 0.0140005442;
double r242311 = r242310 * r242294;
double r242312 = r242309 + r242311;
double r242313 = 0.0008327945;
double r242314 = r242313 * r242298;
double r242315 = r242312 + r242314;
double r242316 = 2.0;
double r242317 = r242316 * r242297;
double r242318 = r242298 * r242282;
double r242319 = r242317 * r242318;
double r242320 = r242315 + r242319;
double r242321 = r242300 / r242320;
double r242322 = r242321 * r242281;
return r242322;
}
double f(double x) {
double r242323 = x;
double r242324 = -14623663576.37637;
bool r242325 = r242323 <= r242324;
double r242326 = 702.9972152077231;
bool r242327 = r242323 <= r242326;
double r242328 = !r242327;
bool r242329 = r242325 || r242328;
double r242330 = 0.2514179000665375;
double r242331 = 1.0;
double r242332 = 3.0;
double r242333 = pow(r242323, r242332);
double r242334 = r242331 / r242333;
double r242335 = r242330 * r242334;
double r242336 = 0.15298196345929327;
double r242337 = 5.0;
double r242338 = pow(r242323, r242337);
double r242339 = r242331 / r242338;
double r242340 = r242336 * r242339;
double r242341 = 0.5;
double r242342 = r242331 / r242323;
double r242343 = r242341 * r242342;
double r242344 = r242340 + r242343;
double r242345 = r242335 + r242344;
double r242346 = 1.0;
double r242347 = 0.1049934947;
double r242348 = r242323 * r242323;
double r242349 = r242347 * r242348;
double r242350 = r242346 + r242349;
double r242351 = 0.0424060604;
double r242352 = r242348 * r242348;
double r242353 = r242351 * r242352;
double r242354 = r242350 + r242353;
double r242355 = 0.0072644182;
double r242356 = r242352 * r242348;
double r242357 = r242355 * r242356;
double r242358 = r242354 + r242357;
double r242359 = 0.0005064034;
double r242360 = r242356 * r242348;
double r242361 = r242359 * r242360;
double r242362 = r242358 + r242361;
double r242363 = 0.0001789971;
double r242364 = r242360 * r242348;
double r242365 = r242363 * r242364;
double r242366 = r242362 + r242365;
double r242367 = 0.7715471019;
double r242368 = r242367 * r242348;
double r242369 = r242346 + r242368;
double r242370 = 0.2909738639;
double r242371 = r242370 * r242352;
double r242372 = r242369 + r242371;
double r242373 = 0.0694555761;
double r242374 = r242373 * r242356;
double r242375 = r242372 + r242374;
double r242376 = 0.0140005442;
double r242377 = r242376 * r242360;
double r242378 = r242375 + r242377;
double r242379 = 0.0008327945;
double r242380 = r242379 * r242364;
double r242381 = r242378 + r242380;
double r242382 = 2.0;
double r242383 = r242382 * r242363;
double r242384 = r242364 * r242348;
double r242385 = r242383 * r242384;
double r242386 = r242381 + r242385;
double r242387 = r242331 / r242386;
double r242388 = r242387 * r242323;
double r242389 = r242366 * r242388;
double r242390 = r242329 ? r242345 : r242389;
return r242390;
}



Bits error versus x
Results
if x < -14623663576.37637 or 702.9972152077231 < x Initial program 59.8
Simplified59.7
Taylor expanded around inf 0.0
if -14623663576.37637 < x < 702.9972152077231Initial program 0.0
rmApplied div-inv0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2020027
(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))