\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 -1972.5019371328153:\\
\;\;\;\;\frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x} + \left(\frac{0.5}{x} + \frac{0.15298196345929327}{\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)}\right)\\
\mathbf{elif}\;x \le 1119.150968415439:\\
\;\;\;\;\left(x \cdot \frac{\left(\left(0.0072644182 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(\left(1 + \left(x \cdot x\right) \cdot 0.1049934947\right) + 0.0424060604 \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(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.0005064034\right) + \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0001789971}{\sqrt{\left(\left(\left(0.0694555761 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.0140005442\right) + 0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(2 \cdot 0.0001789971\right)}}\right) \cdot \frac{1}{\sqrt{\left(\left(\left(0.0694555761 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot 0.0140005442\right) + 0.0008327945 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right)\right) + \left(\left(\left(\left(x \cdot x\right) \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot \left(2 \cdot 0.0001789971\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x} + \left(\frac{0.5}{x} + \frac{0.15298196345929327}{\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot x\right)}\right)\\
\end{array}double f(double x) {
double r6874874 = 1.0;
double r6874875 = 0.1049934947;
double r6874876 = x;
double r6874877 = r6874876 * r6874876;
double r6874878 = r6874875 * r6874877;
double r6874879 = r6874874 + r6874878;
double r6874880 = 0.0424060604;
double r6874881 = r6874877 * r6874877;
double r6874882 = r6874880 * r6874881;
double r6874883 = r6874879 + r6874882;
double r6874884 = 0.0072644182;
double r6874885 = r6874881 * r6874877;
double r6874886 = r6874884 * r6874885;
double r6874887 = r6874883 + r6874886;
double r6874888 = 0.0005064034;
double r6874889 = r6874885 * r6874877;
double r6874890 = r6874888 * r6874889;
double r6874891 = r6874887 + r6874890;
double r6874892 = 0.0001789971;
double r6874893 = r6874889 * r6874877;
double r6874894 = r6874892 * r6874893;
double r6874895 = r6874891 + r6874894;
double r6874896 = 0.7715471019;
double r6874897 = r6874896 * r6874877;
double r6874898 = r6874874 + r6874897;
double r6874899 = 0.2909738639;
double r6874900 = r6874899 * r6874881;
double r6874901 = r6874898 + r6874900;
double r6874902 = 0.0694555761;
double r6874903 = r6874902 * r6874885;
double r6874904 = r6874901 + r6874903;
double r6874905 = 0.0140005442;
double r6874906 = r6874905 * r6874889;
double r6874907 = r6874904 + r6874906;
double r6874908 = 0.0008327945;
double r6874909 = r6874908 * r6874893;
double r6874910 = r6874907 + r6874909;
double r6874911 = 2.0;
double r6874912 = r6874911 * r6874892;
double r6874913 = r6874893 * r6874877;
double r6874914 = r6874912 * r6874913;
double r6874915 = r6874910 + r6874914;
double r6874916 = r6874895 / r6874915;
double r6874917 = r6874916 * r6874876;
return r6874917;
}
double f(double x) {
double r6874918 = x;
double r6874919 = -1972.5019371328153;
bool r6874920 = r6874918 <= r6874919;
double r6874921 = 0.2514179000665375;
double r6874922 = r6874918 * r6874918;
double r6874923 = r6874922 * r6874918;
double r6874924 = r6874921 / r6874923;
double r6874925 = 0.5;
double r6874926 = r6874925 / r6874918;
double r6874927 = 0.15298196345929327;
double r6874928 = r6874922 * r6874923;
double r6874929 = r6874927 / r6874928;
double r6874930 = r6874926 + r6874929;
double r6874931 = r6874924 + r6874930;
double r6874932 = 1119.150968415439;
bool r6874933 = r6874918 <= r6874932;
double r6874934 = 0.0072644182;
double r6874935 = r6874922 * r6874922;
double r6874936 = r6874935 * r6874922;
double r6874937 = r6874934 * r6874936;
double r6874938 = 1.0;
double r6874939 = 0.1049934947;
double r6874940 = r6874922 * r6874939;
double r6874941 = r6874938 + r6874940;
double r6874942 = 0.0424060604;
double r6874943 = r6874942 * r6874935;
double r6874944 = r6874941 + r6874943;
double r6874945 = r6874937 + r6874944;
double r6874946 = r6874922 * r6874936;
double r6874947 = 0.0005064034;
double r6874948 = r6874946 * r6874947;
double r6874949 = r6874945 + r6874948;
double r6874950 = r6874946 * r6874922;
double r6874951 = 0.0001789971;
double r6874952 = r6874950 * r6874951;
double r6874953 = r6874949 + r6874952;
double r6874954 = 0.0694555761;
double r6874955 = r6874954 * r6874936;
double r6874956 = 0.2909738639;
double r6874957 = r6874956 * r6874935;
double r6874958 = 0.7715471019;
double r6874959 = r6874922 * r6874958;
double r6874960 = r6874959 + r6874938;
double r6874961 = r6874957 + r6874960;
double r6874962 = r6874955 + r6874961;
double r6874963 = 0.0140005442;
double r6874964 = r6874946 * r6874963;
double r6874965 = r6874962 + r6874964;
double r6874966 = 0.0008327945;
double r6874967 = r6874966 * r6874950;
double r6874968 = r6874965 + r6874967;
double r6874969 = r6874950 * r6874922;
double r6874970 = 2.0;
double r6874971 = r6874970 * r6874951;
double r6874972 = r6874969 * r6874971;
double r6874973 = r6874968 + r6874972;
double r6874974 = sqrt(r6874973);
double r6874975 = r6874953 / r6874974;
double r6874976 = r6874918 * r6874975;
double r6874977 = r6874938 / r6874974;
double r6874978 = r6874976 * r6874977;
double r6874979 = r6874933 ? r6874978 : r6874931;
double r6874980 = r6874920 ? r6874931 : r6874979;
return r6874980;
}



Bits error versus x
Results
if x < -1972.5019371328153 or 1119.150968415439 < x Initial program 57.8
Taylor expanded around inf 0.0
Simplified0.0
if -1972.5019371328153 < x < 1119.150968415439Initial 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 2019152
(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))