\frac{\left(\left(\left(\left(1.0 + 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 + 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.0 \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 -1254.9677244857505:\\
\;\;\;\;\frac{0.5}{x} + \left(\frac{0.15298196345929327}{{x}^{5}} + \frac{0.2514179000665375}{x \cdot \left(x \cdot x\right)}\right)\\
\mathbf{elif}\;x \le 759.899268558491:\\
\;\;\;\;\frac{\sqrt{\left(\left(x \cdot x\right) \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) \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.0 + 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)\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)} \cdot \left(x \cdot \sqrt{\left(\left(x \cdot x\right) \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) \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.0 + 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)\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)}\right)}{\left(\left(\left(x \cdot x\right) \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) \cdot 0.0008327945 + \left(\left(0.0694555761 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right) + \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.0\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.0140005442\right)\right) + \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \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) \cdot \left(2.0 \cdot 0.0001789971\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5}{x} + \left(\frac{0.15298196345929327}{{x}^{5}} + \frac{0.2514179000665375}{x \cdot \left(x \cdot x\right)}\right)\\
\end{array}double f(double x) {
double r6702759 = 1.0;
double r6702760 = 0.1049934947;
double r6702761 = x;
double r6702762 = r6702761 * r6702761;
double r6702763 = r6702760 * r6702762;
double r6702764 = r6702759 + r6702763;
double r6702765 = 0.0424060604;
double r6702766 = r6702762 * r6702762;
double r6702767 = r6702765 * r6702766;
double r6702768 = r6702764 + r6702767;
double r6702769 = 0.0072644182;
double r6702770 = r6702766 * r6702762;
double r6702771 = r6702769 * r6702770;
double r6702772 = r6702768 + r6702771;
double r6702773 = 0.0005064034;
double r6702774 = r6702770 * r6702762;
double r6702775 = r6702773 * r6702774;
double r6702776 = r6702772 + r6702775;
double r6702777 = 0.0001789971;
double r6702778 = r6702774 * r6702762;
double r6702779 = r6702777 * r6702778;
double r6702780 = r6702776 + r6702779;
double r6702781 = 0.7715471019;
double r6702782 = r6702781 * r6702762;
double r6702783 = r6702759 + r6702782;
double r6702784 = 0.2909738639;
double r6702785 = r6702784 * r6702766;
double r6702786 = r6702783 + r6702785;
double r6702787 = 0.0694555761;
double r6702788 = r6702787 * r6702770;
double r6702789 = r6702786 + r6702788;
double r6702790 = 0.0140005442;
double r6702791 = r6702790 * r6702774;
double r6702792 = r6702789 + r6702791;
double r6702793 = 0.0008327945;
double r6702794 = r6702793 * r6702778;
double r6702795 = r6702792 + r6702794;
double r6702796 = 2.0;
double r6702797 = r6702796 * r6702777;
double r6702798 = r6702778 * r6702762;
double r6702799 = r6702797 * r6702798;
double r6702800 = r6702795 + r6702799;
double r6702801 = r6702780 / r6702800;
double r6702802 = r6702801 * r6702761;
return r6702802;
}
double f(double x) {
double r6702803 = x;
double r6702804 = -1254.9677244857505;
bool r6702805 = r6702803 <= r6702804;
double r6702806 = 0.5;
double r6702807 = r6702806 / r6702803;
double r6702808 = 0.15298196345929327;
double r6702809 = 5.0;
double r6702810 = pow(r6702803, r6702809);
double r6702811 = r6702808 / r6702810;
double r6702812 = 0.2514179000665375;
double r6702813 = r6702803 * r6702803;
double r6702814 = r6702803 * r6702813;
double r6702815 = r6702812 / r6702814;
double r6702816 = r6702811 + r6702815;
double r6702817 = r6702807 + r6702816;
double r6702818 = 759.899268558491;
bool r6702819 = r6702803 <= r6702818;
double r6702820 = r6702813 * r6702813;
double r6702821 = r6702813 * r6702820;
double r6702822 = r6702813 * r6702821;
double r6702823 = r6702813 * r6702822;
double r6702824 = 0.0001789971;
double r6702825 = r6702823 * r6702824;
double r6702826 = 0.0072644182;
double r6702827 = r6702826 * r6702821;
double r6702828 = 1.0;
double r6702829 = 0.1049934947;
double r6702830 = r6702829 * r6702813;
double r6702831 = r6702828 + r6702830;
double r6702832 = 0.0424060604;
double r6702833 = r6702832 * r6702820;
double r6702834 = r6702831 + r6702833;
double r6702835 = r6702827 + r6702834;
double r6702836 = 0.0005064034;
double r6702837 = r6702822 * r6702836;
double r6702838 = r6702835 + r6702837;
double r6702839 = r6702825 + r6702838;
double r6702840 = sqrt(r6702839);
double r6702841 = r6702803 * r6702840;
double r6702842 = r6702840 * r6702841;
double r6702843 = 0.0008327945;
double r6702844 = r6702823 * r6702843;
double r6702845 = 0.0694555761;
double r6702846 = r6702845 * r6702821;
double r6702847 = 0.2909738639;
double r6702848 = r6702847 * r6702820;
double r6702849 = 0.7715471019;
double r6702850 = r6702849 * r6702813;
double r6702851 = r6702850 + r6702828;
double r6702852 = r6702848 + r6702851;
double r6702853 = r6702846 + r6702852;
double r6702854 = 0.0140005442;
double r6702855 = r6702822 * r6702854;
double r6702856 = r6702853 + r6702855;
double r6702857 = r6702844 + r6702856;
double r6702858 = r6702813 * r6702823;
double r6702859 = 2.0;
double r6702860 = r6702859 * r6702824;
double r6702861 = r6702858 * r6702860;
double r6702862 = r6702857 + r6702861;
double r6702863 = r6702842 / r6702862;
double r6702864 = r6702819 ? r6702863 : r6702817;
double r6702865 = r6702805 ? r6702817 : r6702864;
return r6702865;
}



Bits error versus x
Results
if x < -1254.9677244857505 or 759.899268558491 < x Initial program 59.2
Taylor expanded around inf 0.0
Simplified0.0
if -1254.9677244857505 < x < 759.899268558491Initial program 0.0
rmApplied *-un-lft-identity0.0
Applied add-sqr-sqrt0.0
Applied times-frac0.0
Applied associate-*l*0.0
rmApplied associate-*l/0.0
Applied associate-*r/0.0
Final simplification0.0
herbie shell --seed 2019165
(FPCore (x)
:name "Jmat.Real.dawson"
(* (/ (+ (+ (+ (+ (+ 1.0 (* 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 (* 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 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x))