\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 -1207.5387837908731:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\mathbf{elif}\;x \le 727.2141028040012:\\
\;\;\;\;x \cdot \frac{\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(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.1049934947 \cdot \left(x \cdot x\right) + 1\right)\right)\right) + \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 0.0005064034\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)}{\sqrt{\left(\left(x \cdot x\right) \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) \cdot \left(0.0001789971 \cdot 2\right) + \left(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) + \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 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}}}{\sqrt{\left(\left(x \cdot x\right) \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) \cdot \left(0.0001789971 \cdot 2\right) + \left(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) + \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 0.0140005442 + \left(\left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot x\right)\right) \cdot 0.0694555761 + \left(0.2909738639 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(0.7715471019 \cdot \left(x \cdot x\right) + 1\right)\right)\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5 + \frac{0.2514179000665375}{x \cdot x}}{x} + \frac{0.15298196345929327}{{x}^{5}}\\
\end{array}double f(double x) {
double r82746837 = 1.0;
double r82746838 = 0.1049934947;
double r82746839 = x;
double r82746840 = r82746839 * r82746839;
double r82746841 = r82746838 * r82746840;
double r82746842 = r82746837 + r82746841;
double r82746843 = 0.0424060604;
double r82746844 = r82746840 * r82746840;
double r82746845 = r82746843 * r82746844;
double r82746846 = r82746842 + r82746845;
double r82746847 = 0.0072644182;
double r82746848 = r82746844 * r82746840;
double r82746849 = r82746847 * r82746848;
double r82746850 = r82746846 + r82746849;
double r82746851 = 0.0005064034;
double r82746852 = r82746848 * r82746840;
double r82746853 = r82746851 * r82746852;
double r82746854 = r82746850 + r82746853;
double r82746855 = 0.0001789971;
double r82746856 = r82746852 * r82746840;
double r82746857 = r82746855 * r82746856;
double r82746858 = r82746854 + r82746857;
double r82746859 = 0.7715471019;
double r82746860 = r82746859 * r82746840;
double r82746861 = r82746837 + r82746860;
double r82746862 = 0.2909738639;
double r82746863 = r82746862 * r82746844;
double r82746864 = r82746861 + r82746863;
double r82746865 = 0.0694555761;
double r82746866 = r82746865 * r82746848;
double r82746867 = r82746864 + r82746866;
double r82746868 = 0.0140005442;
double r82746869 = r82746868 * r82746852;
double r82746870 = r82746867 + r82746869;
double r82746871 = 0.0008327945;
double r82746872 = r82746871 * r82746856;
double r82746873 = r82746870 + r82746872;
double r82746874 = 2.0;
double r82746875 = r82746874 * r82746855;
double r82746876 = r82746856 * r82746840;
double r82746877 = r82746875 * r82746876;
double r82746878 = r82746873 + r82746877;
double r82746879 = r82746858 / r82746878;
double r82746880 = r82746879 * r82746839;
return r82746880;
}
double f(double x) {
double r82746881 = x;
double r82746882 = -1207.5387837908731;
bool r82746883 = r82746881 <= r82746882;
double r82746884 = 0.5;
double r82746885 = 0.2514179000665375;
double r82746886 = r82746881 * r82746881;
double r82746887 = r82746885 / r82746886;
double r82746888 = r82746884 + r82746887;
double r82746889 = r82746888 / r82746881;
double r82746890 = 0.15298196345929327;
double r82746891 = 5.0;
double r82746892 = pow(r82746881, r82746891);
double r82746893 = r82746890 / r82746892;
double r82746894 = r82746889 + r82746893;
double r82746895 = 727.2141028040012;
bool r82746896 = r82746881 <= r82746895;
double r82746897 = 0.0072644182;
double r82746898 = r82746886 * r82746886;
double r82746899 = r82746898 * r82746886;
double r82746900 = r82746897 * r82746899;
double r82746901 = 0.0424060604;
double r82746902 = r82746901 * r82746898;
double r82746903 = 0.1049934947;
double r82746904 = r82746903 * r82746886;
double r82746905 = 1.0;
double r82746906 = r82746904 + r82746905;
double r82746907 = r82746902 + r82746906;
double r82746908 = r82746900 + r82746907;
double r82746909 = r82746899 * r82746886;
double r82746910 = 0.0005064034;
double r82746911 = r82746909 * r82746910;
double r82746912 = r82746908 + r82746911;
double r82746913 = 0.0001789971;
double r82746914 = r82746909 * r82746886;
double r82746915 = r82746913 * r82746914;
double r82746916 = r82746912 + r82746915;
double r82746917 = r82746886 * r82746914;
double r82746918 = 2.0;
double r82746919 = r82746913 * r82746918;
double r82746920 = r82746917 * r82746919;
double r82746921 = 0.0008327945;
double r82746922 = r82746921 * r82746914;
double r82746923 = 0.0140005442;
double r82746924 = r82746909 * r82746923;
double r82746925 = 0.0694555761;
double r82746926 = r82746899 * r82746925;
double r82746927 = 0.2909738639;
double r82746928 = r82746927 * r82746898;
double r82746929 = 0.7715471019;
double r82746930 = r82746929 * r82746886;
double r82746931 = r82746930 + r82746905;
double r82746932 = r82746928 + r82746931;
double r82746933 = r82746926 + r82746932;
double r82746934 = r82746924 + r82746933;
double r82746935 = r82746922 + r82746934;
double r82746936 = r82746920 + r82746935;
double r82746937 = sqrt(r82746936);
double r82746938 = r82746916 / r82746937;
double r82746939 = r82746938 / r82746937;
double r82746940 = r82746881 * r82746939;
double r82746941 = r82746896 ? r82746940 : r82746894;
double r82746942 = r82746883 ? r82746894 : r82746941;
return r82746942;
}



Bits error versus x
Results
if x < -1207.5387837908731 or 727.2141028040012 < x Initial program 58.0
rmApplied add-sqr-sqrt58.0
Taylor expanded around -inf 0.0
Simplified0.0
if -1207.5387837908731 < x < 727.2141028040012Initial program 0.0
rmApplied add-sqr-sqrt0.0
Applied associate-/r*0.0
Final simplification0.0
herbie shell --seed 2019107 +o rules:numerics
(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))