\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 -783.4081442566906616775668226182460784912:\\
\;\;\;\;\left(\frac{0.5}{x} + \frac{0.2514179000665375252054900556686334311962}{\left(x \cdot x\right) \cdot x}\right) + \frac{0.1529819634592932686700805788859724998474}{\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot x}\\
\mathbf{elif}\;x \le 750.6664899017511061174445785582065582275:\\
\;\;\;\;\left(\frac{\sqrt[3]{\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 + \left(x \cdot x\right) \cdot 0.1049934946999999951788851149103720672429\right) + \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.04240606040000000076517494562722276896238\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) + 1.789971000000000009994005623070734145585 \cdot 10^{-4} \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)}}{\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 \left(x \cdot x\right)\right) \cdot \left(1.789971000000000009994005623070734145585 \cdot 10^{-4} \cdot 2\right) + \left(\left(0.01400054419999999938406531896362139377743 \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(\left(1 + \left(x \cdot x\right) \cdot 0.7715471018999999763821051601553335785866\right) + \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.2909738639000000182122107617033179849386\right) + 0.06945557609999999937322456844412954524159 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) + 8.327945000000000442749725770852364803432 \cdot 10^{-4} \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)} \cdot x\right) \cdot \left(\sqrt[3]{\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 + \left(x \cdot x\right) \cdot 0.1049934946999999951788851149103720672429\right) + \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.04240606040000000076517494562722276896238\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) + 1.789971000000000009994005623070734145585 \cdot 10^{-4} \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)} \cdot \sqrt[3]{\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 + \left(x \cdot x\right) \cdot 0.1049934946999999951788851149103720672429\right) + \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.04240606040000000076517494562722276896238\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) + 1.789971000000000009994005623070734145585 \cdot 10^{-4} \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}:\\
\;\;\;\;\left(\frac{0.5}{x} + \frac{0.2514179000665375252054900556686334311962}{\left(x \cdot x\right) \cdot x}\right) + \frac{0.1529819634592932686700805788859724998474}{\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot x}\\
\end{array}double f(double x) {
double r5673135 = 1.0;
double r5673136 = 0.1049934947;
double r5673137 = x;
double r5673138 = r5673137 * r5673137;
double r5673139 = r5673136 * r5673138;
double r5673140 = r5673135 + r5673139;
double r5673141 = 0.0424060604;
double r5673142 = r5673138 * r5673138;
double r5673143 = r5673141 * r5673142;
double r5673144 = r5673140 + r5673143;
double r5673145 = 0.0072644182;
double r5673146 = r5673142 * r5673138;
double r5673147 = r5673145 * r5673146;
double r5673148 = r5673144 + r5673147;
double r5673149 = 0.0005064034;
double r5673150 = r5673146 * r5673138;
double r5673151 = r5673149 * r5673150;
double r5673152 = r5673148 + r5673151;
double r5673153 = 0.0001789971;
double r5673154 = r5673150 * r5673138;
double r5673155 = r5673153 * r5673154;
double r5673156 = r5673152 + r5673155;
double r5673157 = 0.7715471019;
double r5673158 = r5673157 * r5673138;
double r5673159 = r5673135 + r5673158;
double r5673160 = 0.2909738639;
double r5673161 = r5673160 * r5673142;
double r5673162 = r5673159 + r5673161;
double r5673163 = 0.0694555761;
double r5673164 = r5673163 * r5673146;
double r5673165 = r5673162 + r5673164;
double r5673166 = 0.0140005442;
double r5673167 = r5673166 * r5673150;
double r5673168 = r5673165 + r5673167;
double r5673169 = 0.0008327945;
double r5673170 = r5673169 * r5673154;
double r5673171 = r5673168 + r5673170;
double r5673172 = 2.0;
double r5673173 = r5673172 * r5673153;
double r5673174 = r5673154 * r5673138;
double r5673175 = r5673173 * r5673174;
double r5673176 = r5673171 + r5673175;
double r5673177 = r5673156 / r5673176;
double r5673178 = r5673177 * r5673137;
return r5673178;
}
double f(double x) {
double r5673179 = x;
double r5673180 = -783.4081442566907;
bool r5673181 = r5673179 <= r5673180;
double r5673182 = 0.5;
double r5673183 = r5673182 / r5673179;
double r5673184 = 0.2514179000665375;
double r5673185 = r5673179 * r5673179;
double r5673186 = r5673185 * r5673179;
double r5673187 = r5673184 / r5673186;
double r5673188 = r5673183 + r5673187;
double r5673189 = 0.15298196345929327;
double r5673190 = r5673185 * r5673185;
double r5673191 = r5673190 * r5673179;
double r5673192 = r5673189 / r5673191;
double r5673193 = r5673188 + r5673192;
double r5673194 = 750.6664899017511;
bool r5673195 = r5673179 <= r5673194;
double r5673196 = 0.0072644182;
double r5673197 = r5673185 * r5673190;
double r5673198 = r5673196 * r5673197;
double r5673199 = 1.0;
double r5673200 = 0.1049934947;
double r5673201 = r5673185 * r5673200;
double r5673202 = r5673199 + r5673201;
double r5673203 = 0.0424060604;
double r5673204 = r5673190 * r5673203;
double r5673205 = r5673202 + r5673204;
double r5673206 = r5673198 + r5673205;
double r5673207 = r5673185 * r5673197;
double r5673208 = 0.0005064034;
double r5673209 = r5673207 * r5673208;
double r5673210 = r5673206 + r5673209;
double r5673211 = 0.0001789971;
double r5673212 = r5673207 * r5673185;
double r5673213 = r5673211 * r5673212;
double r5673214 = r5673210 + r5673213;
double r5673215 = cbrt(r5673214);
double r5673216 = r5673212 * r5673185;
double r5673217 = 2.0;
double r5673218 = r5673211 * r5673217;
double r5673219 = r5673216 * r5673218;
double r5673220 = 0.0140005442;
double r5673221 = r5673220 * r5673207;
double r5673222 = 0.7715471019;
double r5673223 = r5673185 * r5673222;
double r5673224 = r5673199 + r5673223;
double r5673225 = 0.2909738639;
double r5673226 = r5673190 * r5673225;
double r5673227 = r5673224 + r5673226;
double r5673228 = 0.0694555761;
double r5673229 = r5673228 * r5673197;
double r5673230 = r5673227 + r5673229;
double r5673231 = r5673221 + r5673230;
double r5673232 = 0.0008327945;
double r5673233 = r5673232 * r5673212;
double r5673234 = r5673231 + r5673233;
double r5673235 = r5673219 + r5673234;
double r5673236 = r5673215 / r5673235;
double r5673237 = r5673236 * r5673179;
double r5673238 = r5673215 * r5673215;
double r5673239 = r5673237 * r5673238;
double r5673240 = r5673195 ? r5673239 : r5673193;
double r5673241 = r5673181 ? r5673193 : r5673240;
return r5673241;
}



Bits error versus x
Results
if x < -783.4081442566907 or 750.6664899017511 < x Initial program 59.1
Taylor expanded around inf 0.0
Simplified0.0
if -783.4081442566907 < x < 750.6664899017511Initial program 0.0
rmApplied *-un-lft-identity0.0
Applied add-cube-cbrt0.0
Applied times-frac0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2019192 +o rules:numerics
(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))