\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 -836.8739538644193:\\
\;\;\;\;\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 654.3238474206455:\\
\;\;\;\;x \cdot \sqrt[3]{\left(\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(\left(x \cdot x\right) \cdot 0.1049934947 + 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.0005064034\right) + 0.0001789971 \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)}{\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(0.0001789971 \cdot 2\right) + \left(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) + \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 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(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right)\right)\right)} \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(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(\left(x \cdot x\right) \cdot 0.1049934947 + 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.0005064034\right) + 0.0001789971 \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)}{\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(0.0001789971 \cdot 2\right) + \left(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) + \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 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(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right)\right)\right)}\right) \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(0.0424060604 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) + \left(\left(x \cdot x\right) \cdot 0.1049934947 + 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.0005064034\right) + 0.0001789971 \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)}{\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(0.0001789971 \cdot 2\right) + \left(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) + \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 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(\left(x \cdot x\right) \cdot 0.7715471019 + 1\right)\right)\right)\right)\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 r6729741 = 1.0;
double r6729742 = 0.1049934947;
double r6729743 = x;
double r6729744 = r6729743 * r6729743;
double r6729745 = r6729742 * r6729744;
double r6729746 = r6729741 + r6729745;
double r6729747 = 0.0424060604;
double r6729748 = r6729744 * r6729744;
double r6729749 = r6729747 * r6729748;
double r6729750 = r6729746 + r6729749;
double r6729751 = 0.0072644182;
double r6729752 = r6729748 * r6729744;
double r6729753 = r6729751 * r6729752;
double r6729754 = r6729750 + r6729753;
double r6729755 = 0.0005064034;
double r6729756 = r6729752 * r6729744;
double r6729757 = r6729755 * r6729756;
double r6729758 = r6729754 + r6729757;
double r6729759 = 0.0001789971;
double r6729760 = r6729756 * r6729744;
double r6729761 = r6729759 * r6729760;
double r6729762 = r6729758 + r6729761;
double r6729763 = 0.7715471019;
double r6729764 = r6729763 * r6729744;
double r6729765 = r6729741 + r6729764;
double r6729766 = 0.2909738639;
double r6729767 = r6729766 * r6729748;
double r6729768 = r6729765 + r6729767;
double r6729769 = 0.0694555761;
double r6729770 = r6729769 * r6729752;
double r6729771 = r6729768 + r6729770;
double r6729772 = 0.0140005442;
double r6729773 = r6729772 * r6729756;
double r6729774 = r6729771 + r6729773;
double r6729775 = 0.0008327945;
double r6729776 = r6729775 * r6729760;
double r6729777 = r6729774 + r6729776;
double r6729778 = 2.0;
double r6729779 = r6729778 * r6729759;
double r6729780 = r6729760 * r6729744;
double r6729781 = r6729779 * r6729780;
double r6729782 = r6729777 + r6729781;
double r6729783 = r6729762 / r6729782;
double r6729784 = r6729783 * r6729743;
return r6729784;
}
double f(double x) {
double r6729785 = x;
double r6729786 = -836.8739538644193;
bool r6729787 = r6729785 <= r6729786;
double r6729788 = 0.2514179000665375;
double r6729789 = r6729785 * r6729785;
double r6729790 = r6729789 * r6729785;
double r6729791 = r6729788 / r6729790;
double r6729792 = 0.5;
double r6729793 = r6729792 / r6729785;
double r6729794 = 0.15298196345929327;
double r6729795 = r6729789 * r6729790;
double r6729796 = r6729794 / r6729795;
double r6729797 = r6729793 + r6729796;
double r6729798 = r6729791 + r6729797;
double r6729799 = 654.3238474206455;
bool r6729800 = r6729785 <= r6729799;
double r6729801 = 0.0072644182;
double r6729802 = r6729789 * r6729789;
double r6729803 = r6729802 * r6729789;
double r6729804 = r6729801 * r6729803;
double r6729805 = 0.0424060604;
double r6729806 = r6729805 * r6729802;
double r6729807 = 0.1049934947;
double r6729808 = r6729789 * r6729807;
double r6729809 = 1.0;
double r6729810 = r6729808 + r6729809;
double r6729811 = r6729806 + r6729810;
double r6729812 = r6729804 + r6729811;
double r6729813 = r6729789 * r6729803;
double r6729814 = 0.0005064034;
double r6729815 = r6729813 * r6729814;
double r6729816 = r6729812 + r6729815;
double r6729817 = 0.0001789971;
double r6729818 = r6729813 * r6729789;
double r6729819 = r6729817 * r6729818;
double r6729820 = r6729816 + r6729819;
double r6729821 = r6729818 * r6729789;
double r6729822 = 2.0;
double r6729823 = r6729817 * r6729822;
double r6729824 = r6729821 * r6729823;
double r6729825 = 0.0008327945;
double r6729826 = r6729825 * r6729818;
double r6729827 = 0.0140005442;
double r6729828 = r6729813 * r6729827;
double r6729829 = 0.0694555761;
double r6729830 = r6729803 * r6729829;
double r6729831 = 0.2909738639;
double r6729832 = r6729831 * r6729802;
double r6729833 = 0.7715471019;
double r6729834 = r6729789 * r6729833;
double r6729835 = r6729834 + r6729809;
double r6729836 = r6729832 + r6729835;
double r6729837 = r6729830 + r6729836;
double r6729838 = r6729828 + r6729837;
double r6729839 = r6729826 + r6729838;
double r6729840 = r6729824 + r6729839;
double r6729841 = r6729820 / r6729840;
double r6729842 = r6729841 * r6729841;
double r6729843 = r6729842 * r6729841;
double r6729844 = cbrt(r6729843);
double r6729845 = r6729785 * r6729844;
double r6729846 = r6729800 ? r6729845 : r6729798;
double r6729847 = r6729787 ? r6729798 : r6729846;
return r6729847;
}



Bits error versus x
Results
if x < -836.8739538644193 or 654.3238474206455 < x Initial program 58.2
Taylor expanded around inf 0.0
Simplified0.0
if -836.8739538644193 < x < 654.3238474206455Initial program 0.0
rmApplied add-cbrt-cube0.0
Final simplification0.0
herbie shell --seed 2019163
(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))