\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 -3720.325581321098:\\
\;\;\;\;\frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x} + \left(\frac{0.15298196345929327}{{x}^{5}} + \frac{0.5}{x}\right)\\
\mathbf{elif}\;x \le 631.893299791396:\\
\;\;\;\;\left(x \cdot \frac{\sqrt{\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 0.0001789971 + \left(\left(0.0072644182 \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.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(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot 0.0005064034\right)}}{\sqrt{\left(2 \cdot 0.0001789971\right) \cdot \left(\left(x \cdot x\right) \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) + \left(0.0008327945 \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(x \cdot x\right)\right)\right) \cdot 0.0694555761 + \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.2909738639 + \left(1 + 0.7715471019 \cdot \left(x \cdot x\right)\right)\right)\right) + 0.0140005442 \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)}}\right) \cdot \frac{\sqrt{\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 0.0001789971 + \left(\left(0.0072644182 \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.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(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right) \cdot 0.0005064034\right)}}{\sqrt{\left(2 \cdot 0.0001789971\right) \cdot \left(\left(x \cdot x\right) \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) + \left(0.0008327945 \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(x \cdot x\right)\right)\right) \cdot 0.0694555761 + \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot 0.2909738639 + \left(1 + 0.7715471019 \cdot \left(x \cdot x\right)\right)\right)\right) + 0.0140005442 \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)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.2514179000665375}{\left(x \cdot x\right) \cdot x} + \left(\frac{0.15298196345929327}{{x}^{5}} + \frac{0.5}{x}\right)\\
\end{array}double f(double x) {
double r4536881 = 1.0;
double r4536882 = 0.1049934947;
double r4536883 = x;
double r4536884 = r4536883 * r4536883;
double r4536885 = r4536882 * r4536884;
double r4536886 = r4536881 + r4536885;
double r4536887 = 0.0424060604;
double r4536888 = r4536884 * r4536884;
double r4536889 = r4536887 * r4536888;
double r4536890 = r4536886 + r4536889;
double r4536891 = 0.0072644182;
double r4536892 = r4536888 * r4536884;
double r4536893 = r4536891 * r4536892;
double r4536894 = r4536890 + r4536893;
double r4536895 = 0.0005064034;
double r4536896 = r4536892 * r4536884;
double r4536897 = r4536895 * r4536896;
double r4536898 = r4536894 + r4536897;
double r4536899 = 0.0001789971;
double r4536900 = r4536896 * r4536884;
double r4536901 = r4536899 * r4536900;
double r4536902 = r4536898 + r4536901;
double r4536903 = 0.7715471019;
double r4536904 = r4536903 * r4536884;
double r4536905 = r4536881 + r4536904;
double r4536906 = 0.2909738639;
double r4536907 = r4536906 * r4536888;
double r4536908 = r4536905 + r4536907;
double r4536909 = 0.0694555761;
double r4536910 = r4536909 * r4536892;
double r4536911 = r4536908 + r4536910;
double r4536912 = 0.0140005442;
double r4536913 = r4536912 * r4536896;
double r4536914 = r4536911 + r4536913;
double r4536915 = 0.0008327945;
double r4536916 = r4536915 * r4536900;
double r4536917 = r4536914 + r4536916;
double r4536918 = 2.0;
double r4536919 = r4536918 * r4536899;
double r4536920 = r4536900 * r4536884;
double r4536921 = r4536919 * r4536920;
double r4536922 = r4536917 + r4536921;
double r4536923 = r4536902 / r4536922;
double r4536924 = r4536923 * r4536883;
return r4536924;
}
double f(double x) {
double r4536925 = x;
double r4536926 = -3720.325581321098;
bool r4536927 = r4536925 <= r4536926;
double r4536928 = 0.2514179000665375;
double r4536929 = r4536925 * r4536925;
double r4536930 = r4536929 * r4536925;
double r4536931 = r4536928 / r4536930;
double r4536932 = 0.15298196345929327;
double r4536933 = 5.0;
double r4536934 = pow(r4536925, r4536933);
double r4536935 = r4536932 / r4536934;
double r4536936 = 0.5;
double r4536937 = r4536936 / r4536925;
double r4536938 = r4536935 + r4536937;
double r4536939 = r4536931 + r4536938;
double r4536940 = 631.893299791396;
bool r4536941 = r4536925 <= r4536940;
double r4536942 = r4536929 * r4536929;
double r4536943 = r4536929 * r4536942;
double r4536944 = r4536929 * r4536943;
double r4536945 = r4536944 * r4536929;
double r4536946 = 0.0001789971;
double r4536947 = r4536945 * r4536946;
double r4536948 = 0.0072644182;
double r4536949 = r4536948 * r4536943;
double r4536950 = 1.0;
double r4536951 = 0.1049934947;
double r4536952 = r4536929 * r4536951;
double r4536953 = r4536950 + r4536952;
double r4536954 = 0.0424060604;
double r4536955 = r4536954 * r4536942;
double r4536956 = r4536953 + r4536955;
double r4536957 = r4536949 + r4536956;
double r4536958 = 0.0005064034;
double r4536959 = r4536944 * r4536958;
double r4536960 = r4536957 + r4536959;
double r4536961 = r4536947 + r4536960;
double r4536962 = sqrt(r4536961);
double r4536963 = 2.0;
double r4536964 = r4536963 * r4536946;
double r4536965 = r4536929 * r4536945;
double r4536966 = r4536964 * r4536965;
double r4536967 = 0.0008327945;
double r4536968 = r4536967 * r4536945;
double r4536969 = 0.0694555761;
double r4536970 = r4536943 * r4536969;
double r4536971 = 0.2909738639;
double r4536972 = r4536942 * r4536971;
double r4536973 = 0.7715471019;
double r4536974 = r4536973 * r4536929;
double r4536975 = r4536950 + r4536974;
double r4536976 = r4536972 + r4536975;
double r4536977 = r4536970 + r4536976;
double r4536978 = 0.0140005442;
double r4536979 = r4536978 * r4536944;
double r4536980 = r4536977 + r4536979;
double r4536981 = r4536968 + r4536980;
double r4536982 = r4536966 + r4536981;
double r4536983 = sqrt(r4536982);
double r4536984 = r4536962 / r4536983;
double r4536985 = r4536925 * r4536984;
double r4536986 = r4536985 * r4536984;
double r4536987 = r4536941 ? r4536986 : r4536939;
double r4536988 = r4536927 ? r4536939 : r4536987;
return r4536988;
}



Bits error versus x
Results
if x < -3720.325581321098 or 631.893299791396 < x Initial program 58.1
Taylor expanded around inf 0.0
Simplified0.0
if -3720.325581321098 < x < 631.893299791396Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied add-sqr-sqrt0.0
Applied times-frac0.0
Applied associate-*l*0.0
Final simplification0.0
herbie shell --seed 2019134
(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))