\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 -26828843184951967744:\\
\;\;\;\;\frac{0.1529819634592932686700805788859724998474}{{x}^{5}} + \left(\frac{0.2514179000665375252054900556686334311962}{x \cdot \left(x \cdot x\right)} + \frac{0.5}{x}\right)\\
\mathbf{elif}\;x \le 729.8020590425830960157327353954315185547:\\
\;\;\;\;x \cdot \frac{\frac{\left(\left(x \cdot x\right) \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)\right) \cdot 1.789971000000000009994005623070734145585 \cdot 10^{-4} + \left(\left(0.007264418199999999985194687468492702464573 \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 + 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)\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 5.064034000000000243502107366566633572802 \cdot 10^{-4}\right)}{\sqrt{\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \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)\right)\right) \cdot \left(2 \cdot 1.789971000000000009994005623070734145585 \cdot 10^{-4}\right) + \left(\left(\left(x \cdot x\right) \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)\right) \cdot 8.327945000000000442749725770852364803432 \cdot 10^{-4} + \left(\left(\left(\left(\left(x \cdot x\right) \cdot 0.7715471018999999763821051601553335785866 + 1\right) + 0.2909738639000000182122107617033179849386 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.06945557609999999937322456844412954524159\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.01400054419999999938406531896362139377743\right)\right)}}}{\sqrt{\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \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)\right)\right) \cdot \left(2 \cdot 1.789971000000000009994005623070734145585 \cdot 10^{-4}\right) + \left(\left(\left(x \cdot x\right) \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)\right) \cdot 8.327945000000000442749725770852364803432 \cdot 10^{-4} + \left(\left(\left(\left(\left(x \cdot x\right) \cdot 0.7715471018999999763821051601553335785866 + 1\right) + 0.2909738639000000182122107617033179849386 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.06945557609999999937322456844412954524159\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.01400054419999999938406531896362139377743\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.1529819634592932686700805788859724998474}{{x}^{5}} + \left(\frac{0.2514179000665375252054900556686334311962}{x \cdot \left(x \cdot x\right)} + \frac{0.5}{x}\right)\\
\end{array}double f(double x) {
double r7390095 = 1.0;
double r7390096 = 0.1049934947;
double r7390097 = x;
double r7390098 = r7390097 * r7390097;
double r7390099 = r7390096 * r7390098;
double r7390100 = r7390095 + r7390099;
double r7390101 = 0.0424060604;
double r7390102 = r7390098 * r7390098;
double r7390103 = r7390101 * r7390102;
double r7390104 = r7390100 + r7390103;
double r7390105 = 0.0072644182;
double r7390106 = r7390102 * r7390098;
double r7390107 = r7390105 * r7390106;
double r7390108 = r7390104 + r7390107;
double r7390109 = 0.0005064034;
double r7390110 = r7390106 * r7390098;
double r7390111 = r7390109 * r7390110;
double r7390112 = r7390108 + r7390111;
double r7390113 = 0.0001789971;
double r7390114 = r7390110 * r7390098;
double r7390115 = r7390113 * r7390114;
double r7390116 = r7390112 + r7390115;
double r7390117 = 0.7715471019;
double r7390118 = r7390117 * r7390098;
double r7390119 = r7390095 + r7390118;
double r7390120 = 0.2909738639;
double r7390121 = r7390120 * r7390102;
double r7390122 = r7390119 + r7390121;
double r7390123 = 0.0694555761;
double r7390124 = r7390123 * r7390106;
double r7390125 = r7390122 + r7390124;
double r7390126 = 0.0140005442;
double r7390127 = r7390126 * r7390110;
double r7390128 = r7390125 + r7390127;
double r7390129 = 0.0008327945;
double r7390130 = r7390129 * r7390114;
double r7390131 = r7390128 + r7390130;
double r7390132 = 2.0;
double r7390133 = r7390132 * r7390113;
double r7390134 = r7390114 * r7390098;
double r7390135 = r7390133 * r7390134;
double r7390136 = r7390131 + r7390135;
double r7390137 = r7390116 / r7390136;
double r7390138 = r7390137 * r7390097;
return r7390138;
}
double f(double x) {
double r7390139 = x;
double r7390140 = -2.6828843184951968e+19;
bool r7390141 = r7390139 <= r7390140;
double r7390142 = 0.15298196345929327;
double r7390143 = 5.0;
double r7390144 = pow(r7390139, r7390143);
double r7390145 = r7390142 / r7390144;
double r7390146 = 0.2514179000665375;
double r7390147 = r7390139 * r7390139;
double r7390148 = r7390139 * r7390147;
double r7390149 = r7390146 / r7390148;
double r7390150 = 0.5;
double r7390151 = r7390150 / r7390139;
double r7390152 = r7390149 + r7390151;
double r7390153 = r7390145 + r7390152;
double r7390154 = 729.8020590425831;
bool r7390155 = r7390139 <= r7390154;
double r7390156 = r7390147 * r7390147;
double r7390157 = r7390147 * r7390156;
double r7390158 = r7390147 * r7390157;
double r7390159 = r7390147 * r7390158;
double r7390160 = 0.0001789971;
double r7390161 = r7390159 * r7390160;
double r7390162 = 0.0072644182;
double r7390163 = r7390162 * r7390157;
double r7390164 = 1.0;
double r7390165 = 0.1049934947;
double r7390166 = r7390165 * r7390147;
double r7390167 = r7390164 + r7390166;
double r7390168 = 0.0424060604;
double r7390169 = r7390168 * r7390156;
double r7390170 = r7390167 + r7390169;
double r7390171 = r7390163 + r7390170;
double r7390172 = 0.0005064034;
double r7390173 = r7390158 * r7390172;
double r7390174 = r7390171 + r7390173;
double r7390175 = r7390161 + r7390174;
double r7390176 = r7390147 * r7390159;
double r7390177 = 2.0;
double r7390178 = r7390177 * r7390160;
double r7390179 = r7390176 * r7390178;
double r7390180 = 0.0008327945;
double r7390181 = r7390159 * r7390180;
double r7390182 = 0.7715471019;
double r7390183 = r7390147 * r7390182;
double r7390184 = r7390183 + r7390164;
double r7390185 = 0.2909738639;
double r7390186 = r7390185 * r7390156;
double r7390187 = r7390184 + r7390186;
double r7390188 = 0.0694555761;
double r7390189 = r7390157 * r7390188;
double r7390190 = r7390187 + r7390189;
double r7390191 = 0.0140005442;
double r7390192 = r7390158 * r7390191;
double r7390193 = r7390190 + r7390192;
double r7390194 = r7390181 + r7390193;
double r7390195 = r7390179 + r7390194;
double r7390196 = sqrt(r7390195);
double r7390197 = r7390175 / r7390196;
double r7390198 = r7390197 / r7390196;
double r7390199 = r7390139 * r7390198;
double r7390200 = r7390155 ? r7390199 : r7390153;
double r7390201 = r7390141 ? r7390153 : r7390200;
return r7390201;
}



Bits error versus x
Results
if x < -2.6828843184951968e+19 or 729.8020590425831 < x Initial program 60.6
Taylor expanded around inf 0.0
Simplified0.0
if -2.6828843184951968e+19 < x < 729.8020590425831Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied associate-/r*0.0
Final simplification0.0
herbie shell --seed 2019169
(FPCore (x)
:name "Jmat.Real.dawson"
(* (/ (+ (+ (+ (+ (+ 1.0 (* 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 (* 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 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x))