\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 -6753.5048671236073 \lor \neg \left(x \le 7293.4466097268923\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}:\\
\;\;\;\;\frac{\sqrt{\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)}}{1} \cdot \frac{x \cdot \sqrt{\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({x}^{2} \cdot \left({x}^{2} \cdot \left({x}^{2} \cdot \left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right)\right)\right)\right) \cdot \left(2 \cdot 1.789971 \cdot 10^{-4}\right) + \left(\left(\left({x}^{2} \cdot \left(x \cdot {x}^{3}\right)\right) \cdot 0.069455576099999999 + 1\right) + {x}^{2} \cdot \left(0.77154710189999998 + 0.29097386390000002 \cdot \left(x \cdot x\right)\right)\right)\right) + {x}^{2} \cdot \left(0.014000544199999999 \cdot \left(\left(\left(x \cdot x\right) \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)\right) + 8.32794500000000044 \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)}\\
\end{array}double f(double x) {
double r243116 = 1.0;
double r243117 = 0.1049934947;
double r243118 = x;
double r243119 = r243118 * r243118;
double r243120 = r243117 * r243119;
double r243121 = r243116 + r243120;
double r243122 = 0.0424060604;
double r243123 = r243119 * r243119;
double r243124 = r243122 * r243123;
double r243125 = r243121 + r243124;
double r243126 = 0.0072644182;
double r243127 = r243123 * r243119;
double r243128 = r243126 * r243127;
double r243129 = r243125 + r243128;
double r243130 = 0.0005064034;
double r243131 = r243127 * r243119;
double r243132 = r243130 * r243131;
double r243133 = r243129 + r243132;
double r243134 = 0.0001789971;
double r243135 = r243131 * r243119;
double r243136 = r243134 * r243135;
double r243137 = r243133 + r243136;
double r243138 = 0.7715471019;
double r243139 = r243138 * r243119;
double r243140 = r243116 + r243139;
double r243141 = 0.2909738639;
double r243142 = r243141 * r243123;
double r243143 = r243140 + r243142;
double r243144 = 0.0694555761;
double r243145 = r243144 * r243127;
double r243146 = r243143 + r243145;
double r243147 = 0.0140005442;
double r243148 = r243147 * r243131;
double r243149 = r243146 + r243148;
double r243150 = 0.0008327945;
double r243151 = r243150 * r243135;
double r243152 = r243149 + r243151;
double r243153 = 2.0;
double r243154 = r243153 * r243134;
double r243155 = r243135 * r243119;
double r243156 = r243154 * r243155;
double r243157 = r243152 + r243156;
double r243158 = r243137 / r243157;
double r243159 = r243158 * r243118;
return r243159;
}
double f(double x) {
double r243160 = x;
double r243161 = -6753.504867123607;
bool r243162 = r243160 <= r243161;
double r243163 = 7293.446609726892;
bool r243164 = r243160 <= r243163;
double r243165 = !r243164;
bool r243166 = r243162 || r243165;
double r243167 = 0.2514179000665375;
double r243168 = 1.0;
double r243169 = 3.0;
double r243170 = pow(r243160, r243169);
double r243171 = r243168 / r243170;
double r243172 = r243167 * r243171;
double r243173 = 0.15298196345929327;
double r243174 = 5.0;
double r243175 = pow(r243160, r243174);
double r243176 = r243168 / r243175;
double r243177 = r243173 * r243176;
double r243178 = 0.5;
double r243179 = r243168 / r243160;
double r243180 = r243178 * r243179;
double r243181 = r243177 + r243180;
double r243182 = r243172 + r243181;
double r243183 = 1.0;
double r243184 = 0.1049934947;
double r243185 = r243160 * r243160;
double r243186 = r243184 * r243185;
double r243187 = r243183 + r243186;
double r243188 = 0.0424060604;
double r243189 = r243185 * r243185;
double r243190 = r243188 * r243189;
double r243191 = r243187 + r243190;
double r243192 = 0.0072644182;
double r243193 = r243189 * r243185;
double r243194 = r243192 * r243193;
double r243195 = r243191 + r243194;
double r243196 = 0.0005064034;
double r243197 = r243193 * r243185;
double r243198 = r243196 * r243197;
double r243199 = r243195 + r243198;
double r243200 = 0.0001789971;
double r243201 = r243197 * r243185;
double r243202 = r243200 * r243201;
double r243203 = r243199 + r243202;
double r243204 = sqrt(r243203);
double r243205 = r243204 / r243168;
double r243206 = r243160 * r243204;
double r243207 = 2.0;
double r243208 = pow(r243160, r243207);
double r243209 = r243160 * r243170;
double r243210 = r243208 * r243209;
double r243211 = r243208 * r243210;
double r243212 = r243208 * r243211;
double r243213 = r243208 * r243212;
double r243214 = 2.0;
double r243215 = r243214 * r243200;
double r243216 = r243213 * r243215;
double r243217 = 0.0694555761;
double r243218 = r243210 * r243217;
double r243219 = r243218 + r243183;
double r243220 = 0.7715471019;
double r243221 = 0.2909738639;
double r243222 = r243221 * r243185;
double r243223 = r243220 + r243222;
double r243224 = r243208 * r243223;
double r243225 = r243219 + r243224;
double r243226 = r243216 + r243225;
double r243227 = 0.0140005442;
double r243228 = r243185 * r243160;
double r243229 = r243228 * r243228;
double r243230 = r243227 * r243229;
double r243231 = 0.0008327945;
double r243232 = r243231 * r243197;
double r243233 = r243230 + r243232;
double r243234 = r243208 * r243233;
double r243235 = r243226 + r243234;
double r243236 = r243206 / r243235;
double r243237 = r243205 * r243236;
double r243238 = r243166 ? r243182 : r243237;
return r243238;
}



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