\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 -660441.7377912633:\\
\;\;\;\;\frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x} + \left(\frac{0.15298196345929327}{{x}^{5}} + \frac{0.5}{x}\right)\\
\mathbf{elif}\;x \le 905.9497489091251:\\
\;\;\;\;\left(x \cdot \frac{\left(0.0005064034 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) + \left(\left(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(\left(x \cdot x\right) \cdot 0.1049934947 + 1\right)\right) + 0.0072644182 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) + 0.0001789971 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right)}{\sqrt{\left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0008327945 + \left(\left(0.0694555761 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + \left(\left(1 + \left(x \cdot x\right) \cdot 0.7715471019\right) + 0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot 0.0140005442\right)\right) + \left(0.0001789971 \cdot 2\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right)\right)}}\right) \cdot \frac{1}{\sqrt{\left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0008327945 + \left(\left(0.0694555761 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + \left(\left(1 + \left(x \cdot x\right) \cdot 0.7715471019\right) + 0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot 0.0140005442\right)\right) + \left(0.0001789971 \cdot 2\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot \left(x \cdot x\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x} + \left(\frac{0.15298196345929327}{{x}^{5}} + \frac{0.5}{x}\right)\\
\end{array}double f(double x) {
double r3692966 = 1.0;
double r3692967 = 0.1049934947;
double r3692968 = x;
double r3692969 = r3692968 * r3692968;
double r3692970 = r3692967 * r3692969;
double r3692971 = r3692966 + r3692970;
double r3692972 = 0.0424060604;
double r3692973 = r3692969 * r3692969;
double r3692974 = r3692972 * r3692973;
double r3692975 = r3692971 + r3692974;
double r3692976 = 0.0072644182;
double r3692977 = r3692973 * r3692969;
double r3692978 = r3692976 * r3692977;
double r3692979 = r3692975 + r3692978;
double r3692980 = 0.0005064034;
double r3692981 = r3692977 * r3692969;
double r3692982 = r3692980 * r3692981;
double r3692983 = r3692979 + r3692982;
double r3692984 = 0.0001789971;
double r3692985 = r3692981 * r3692969;
double r3692986 = r3692984 * r3692985;
double r3692987 = r3692983 + r3692986;
double r3692988 = 0.7715471019;
double r3692989 = r3692988 * r3692969;
double r3692990 = r3692966 + r3692989;
double r3692991 = 0.2909738639;
double r3692992 = r3692991 * r3692973;
double r3692993 = r3692990 + r3692992;
double r3692994 = 0.0694555761;
double r3692995 = r3692994 * r3692977;
double r3692996 = r3692993 + r3692995;
double r3692997 = 0.0140005442;
double r3692998 = r3692997 * r3692981;
double r3692999 = r3692996 + r3692998;
double r3693000 = 0.0008327945;
double r3693001 = r3693000 * r3692985;
double r3693002 = r3692999 + r3693001;
double r3693003 = 2.0;
double r3693004 = r3693003 * r3692984;
double r3693005 = r3692985 * r3692969;
double r3693006 = r3693004 * r3693005;
double r3693007 = r3693002 + r3693006;
double r3693008 = r3692987 / r3693007;
double r3693009 = r3693008 * r3692968;
return r3693009;
}
double f(double x) {
double r3693010 = x;
double r3693011 = -660441.7377912633;
bool r3693012 = r3693010 <= r3693011;
double r3693013 = 0.2514179000665375;
double r3693014 = r3693010 * r3693010;
double r3693015 = r3693014 * r3693010;
double r3693016 = r3693013 / r3693015;
double r3693017 = 0.15298196345929327;
double r3693018 = 5.0;
double r3693019 = pow(r3693010, r3693018);
double r3693020 = r3693017 / r3693019;
double r3693021 = 0.5;
double r3693022 = r3693021 / r3693010;
double r3693023 = r3693020 + r3693022;
double r3693024 = r3693016 + r3693023;
double r3693025 = 905.9497489091251;
bool r3693026 = r3693010 <= r3693025;
double r3693027 = 0.0005064034;
double r3693028 = r3693014 * r3693014;
double r3693029 = r3693014 * r3693028;
double r3693030 = r3693014 * r3693029;
double r3693031 = r3693027 * r3693030;
double r3693032 = 0.0424060604;
double r3693033 = r3693032 * r3693028;
double r3693034 = 0.1049934947;
double r3693035 = r3693014 * r3693034;
double r3693036 = 1.0;
double r3693037 = r3693035 + r3693036;
double r3693038 = r3693033 + r3693037;
double r3693039 = 0.0072644182;
double r3693040 = r3693039 * r3693029;
double r3693041 = r3693038 + r3693040;
double r3693042 = r3693031 + r3693041;
double r3693043 = 0.0001789971;
double r3693044 = r3693030 * r3693014;
double r3693045 = r3693043 * r3693044;
double r3693046 = r3693042 + r3693045;
double r3693047 = 0.0008327945;
double r3693048 = r3693044 * r3693047;
double r3693049 = 0.0694555761;
double r3693050 = r3693049 * r3693029;
double r3693051 = 0.7715471019;
double r3693052 = r3693014 * r3693051;
double r3693053 = r3693036 + r3693052;
double r3693054 = 0.2909738639;
double r3693055 = r3693054 * r3693028;
double r3693056 = r3693053 + r3693055;
double r3693057 = r3693050 + r3693056;
double r3693058 = 0.0140005442;
double r3693059 = r3693030 * r3693058;
double r3693060 = r3693057 + r3693059;
double r3693061 = r3693048 + r3693060;
double r3693062 = 2.0;
double r3693063 = r3693043 * r3693062;
double r3693064 = r3693014 * r3693044;
double r3693065 = r3693063 * r3693064;
double r3693066 = r3693061 + r3693065;
double r3693067 = sqrt(r3693066);
double r3693068 = r3693046 / r3693067;
double r3693069 = r3693010 * r3693068;
double r3693070 = r3693036 / r3693067;
double r3693071 = r3693069 * r3693070;
double r3693072 = r3693026 ? r3693071 : r3693024;
double r3693073 = r3693012 ? r3693024 : r3693072;
return r3693073;
}



Bits error versus x
Results
if x < -660441.7377912633 or 905.9497489091251 < x Initial program 58.6
Taylor expanded around -inf 0.0
Simplified0.0
if -660441.7377912633 < x < 905.9497489091251Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied *-un-lft-identity0.0
Applied times-frac0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2019153
(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))