\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 -201055.117688729166 \lor \neg \left(x \le 566.32829468601153\right):\\
\;\;\;\;\mathsf{fma}\left(0.25141790006653753, \frac{1}{{x}^{3}}, \mathsf{fma}\left(0.1529819634592933, \frac{1}{{x}^{5}}, \frac{0.5}{x}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\sqrt[3]{\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 \sqrt[3]{\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)}}\right) \cdot \left(\sqrt[3]{\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\right)\\
\end{array}double f(double x) {
double r207913 = 1.0;
double r207914 = 0.1049934947;
double r207915 = x;
double r207916 = r207915 * r207915;
double r207917 = r207914 * r207916;
double r207918 = r207913 + r207917;
double r207919 = 0.0424060604;
double r207920 = r207916 * r207916;
double r207921 = r207919 * r207920;
double r207922 = r207918 + r207921;
double r207923 = 0.0072644182;
double r207924 = r207920 * r207916;
double r207925 = r207923 * r207924;
double r207926 = r207922 + r207925;
double r207927 = 0.0005064034;
double r207928 = r207924 * r207916;
double r207929 = r207927 * r207928;
double r207930 = r207926 + r207929;
double r207931 = 0.0001789971;
double r207932 = r207928 * r207916;
double r207933 = r207931 * r207932;
double r207934 = r207930 + r207933;
double r207935 = 0.7715471019;
double r207936 = r207935 * r207916;
double r207937 = r207913 + r207936;
double r207938 = 0.2909738639;
double r207939 = r207938 * r207920;
double r207940 = r207937 + r207939;
double r207941 = 0.0694555761;
double r207942 = r207941 * r207924;
double r207943 = r207940 + r207942;
double r207944 = 0.0140005442;
double r207945 = r207944 * r207928;
double r207946 = r207943 + r207945;
double r207947 = 0.0008327945;
double r207948 = r207947 * r207932;
double r207949 = r207946 + r207948;
double r207950 = 2.0;
double r207951 = r207950 * r207931;
double r207952 = r207932 * r207916;
double r207953 = r207951 * r207952;
double r207954 = r207949 + r207953;
double r207955 = r207934 / r207954;
double r207956 = r207955 * r207915;
return r207956;
}
double f(double x) {
double r207957 = x;
double r207958 = -201055.11768872917;
bool r207959 = r207957 <= r207958;
double r207960 = 566.3282946860115;
bool r207961 = r207957 <= r207960;
double r207962 = !r207961;
bool r207963 = r207959 || r207962;
double r207964 = 0.2514179000665375;
double r207965 = 1.0;
double r207966 = 3.0;
double r207967 = pow(r207957, r207966);
double r207968 = r207965 / r207967;
double r207969 = 0.15298196345929327;
double r207970 = 5.0;
double r207971 = pow(r207957, r207970);
double r207972 = r207965 / r207971;
double r207973 = 0.5;
double r207974 = r207973 / r207957;
double r207975 = fma(r207969, r207972, r207974);
double r207976 = fma(r207964, r207968, r207975);
double r207977 = 1.0;
double r207978 = 0.1049934947;
double r207979 = r207957 * r207957;
double r207980 = r207978 * r207979;
double r207981 = r207977 + r207980;
double r207982 = 0.0424060604;
double r207983 = r207979 * r207979;
double r207984 = r207982 * r207983;
double r207985 = r207981 + r207984;
double r207986 = 0.0072644182;
double r207987 = r207983 * r207979;
double r207988 = r207986 * r207987;
double r207989 = r207985 + r207988;
double r207990 = 0.0005064034;
double r207991 = r207987 * r207979;
double r207992 = r207990 * r207991;
double r207993 = r207989 + r207992;
double r207994 = 0.0001789971;
double r207995 = r207991 * r207979;
double r207996 = r207994 * r207995;
double r207997 = r207993 + r207996;
double r207998 = 0.7715471019;
double r207999 = r207998 * r207979;
double r208000 = r207977 + r207999;
double r208001 = 0.2909738639;
double r208002 = r208001 * r207983;
double r208003 = r208000 + r208002;
double r208004 = 0.0694555761;
double r208005 = r208004 * r207987;
double r208006 = r208003 + r208005;
double r208007 = 0.0140005442;
double r208008 = r208007 * r207991;
double r208009 = r208006 + r208008;
double r208010 = 0.0008327945;
double r208011 = r208010 * r207995;
double r208012 = r208009 + r208011;
double r208013 = 2.0;
double r208014 = r208013 * r207994;
double r208015 = r207995 * r207979;
double r208016 = r208014 * r208015;
double r208017 = r208012 + r208016;
double r208018 = r207997 / r208017;
double r208019 = cbrt(r208018);
double r208020 = r208019 * r208019;
double r208021 = r208019 * r207957;
double r208022 = r208020 * r208021;
double r208023 = r207963 ? r207976 : r208022;
return r208023;
}



Bits error versus x
if x < -201055.11768872917 or 566.3282946860115 < x Initial program 59.3
Taylor expanded around inf 0.0
Simplified0.0
if -201055.11768872917 < x < 566.3282946860115Initial program 0.0
rmApplied add-cube-cbrt0.1
Applied associate-*l*0.1
Final simplification0.0
herbie shell --seed 2020065 +o rules:numerics
(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))