\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 -660441.7377912633:\\
\;\;\;\;\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 905.9497489091251:\\
\;\;\;\;\left(x \cdot \frac{\left(0.0005064034 \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) + \left(\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) + 0.0072644182 \cdot \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\right)\right) + 0.0001789971 \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)}{\sqrt{\left(\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.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(\left(1 + \left(x \cdot x\right) \cdot 0.7715471019\right) + 0.2909738639 \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.0140005442\right)\right) + \left(0.0001789971 \cdot 2\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)}}\right) \cdot \frac{1}{\sqrt{\left(\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.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(\left(1 + \left(x \cdot x\right) \cdot 0.7715471019\right) + 0.2909738639 \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.0140005442\right)\right) + \left(0.0001789971 \cdot 2\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)}}\\
\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 r3509627 = 1.0;
double r3509628 = 0.1049934947;
double r3509629 = x;
double r3509630 = r3509629 * r3509629;
double r3509631 = r3509628 * r3509630;
double r3509632 = r3509627 + r3509631;
double r3509633 = 0.0424060604;
double r3509634 = r3509630 * r3509630;
double r3509635 = r3509633 * r3509634;
double r3509636 = r3509632 + r3509635;
double r3509637 = 0.0072644182;
double r3509638 = r3509634 * r3509630;
double r3509639 = r3509637 * r3509638;
double r3509640 = r3509636 + r3509639;
double r3509641 = 0.0005064034;
double r3509642 = r3509638 * r3509630;
double r3509643 = r3509641 * r3509642;
double r3509644 = r3509640 + r3509643;
double r3509645 = 0.0001789971;
double r3509646 = r3509642 * r3509630;
double r3509647 = r3509645 * r3509646;
double r3509648 = r3509644 + r3509647;
double r3509649 = 0.7715471019;
double r3509650 = r3509649 * r3509630;
double r3509651 = r3509627 + r3509650;
double r3509652 = 0.2909738639;
double r3509653 = r3509652 * r3509634;
double r3509654 = r3509651 + r3509653;
double r3509655 = 0.0694555761;
double r3509656 = r3509655 * r3509638;
double r3509657 = r3509654 + r3509656;
double r3509658 = 0.0140005442;
double r3509659 = r3509658 * r3509642;
double r3509660 = r3509657 + r3509659;
double r3509661 = 0.0008327945;
double r3509662 = r3509661 * r3509646;
double r3509663 = r3509660 + r3509662;
double r3509664 = 2.0;
double r3509665 = r3509664 * r3509645;
double r3509666 = r3509646 * r3509630;
double r3509667 = r3509665 * r3509666;
double r3509668 = r3509663 + r3509667;
double r3509669 = r3509648 / r3509668;
double r3509670 = r3509669 * r3509629;
return r3509670;
}
double f(double x) {
double r3509671 = x;
double r3509672 = -660441.7377912633;
bool r3509673 = r3509671 <= r3509672;
double r3509674 = 0.2514179000665375;
double r3509675 = r3509671 * r3509671;
double r3509676 = r3509675 * r3509671;
double r3509677 = r3509674 / r3509676;
double r3509678 = 0.15298196345929327;
double r3509679 = 5.0;
double r3509680 = pow(r3509671, r3509679);
double r3509681 = r3509678 / r3509680;
double r3509682 = 0.5;
double r3509683 = r3509682 / r3509671;
double r3509684 = r3509681 + r3509683;
double r3509685 = r3509677 + r3509684;
double r3509686 = 905.9497489091251;
bool r3509687 = r3509671 <= r3509686;
double r3509688 = 0.0005064034;
double r3509689 = r3509675 * r3509675;
double r3509690 = r3509675 * r3509689;
double r3509691 = r3509675 * r3509690;
double r3509692 = r3509688 * r3509691;
double r3509693 = 0.0424060604;
double r3509694 = r3509693 * r3509689;
double r3509695 = 0.1049934947;
double r3509696 = r3509675 * r3509695;
double r3509697 = 1.0;
double r3509698 = r3509696 + r3509697;
double r3509699 = r3509694 + r3509698;
double r3509700 = 0.0072644182;
double r3509701 = r3509700 * r3509690;
double r3509702 = r3509699 + r3509701;
double r3509703 = r3509692 + r3509702;
double r3509704 = 0.0001789971;
double r3509705 = r3509691 * r3509675;
double r3509706 = r3509704 * r3509705;
double r3509707 = r3509703 + r3509706;
double r3509708 = 0.0008327945;
double r3509709 = r3509705 * r3509708;
double r3509710 = 0.0694555761;
double r3509711 = r3509710 * r3509690;
double r3509712 = 0.7715471019;
double r3509713 = r3509675 * r3509712;
double r3509714 = r3509697 + r3509713;
double r3509715 = 0.2909738639;
double r3509716 = r3509715 * r3509689;
double r3509717 = r3509714 + r3509716;
double r3509718 = r3509711 + r3509717;
double r3509719 = 0.0140005442;
double r3509720 = r3509691 * r3509719;
double r3509721 = r3509718 + r3509720;
double r3509722 = r3509709 + r3509721;
double r3509723 = 2.0;
double r3509724 = r3509704 * r3509723;
double r3509725 = r3509675 * r3509705;
double r3509726 = r3509724 * r3509725;
double r3509727 = r3509722 + r3509726;
double r3509728 = sqrt(r3509727);
double r3509729 = r3509707 / r3509728;
double r3509730 = r3509671 * r3509729;
double r3509731 = r3509697 / r3509728;
double r3509732 = r3509730 * r3509731;
double r3509733 = r3509687 ? r3509732 : r3509685;
double r3509734 = r3509673 ? r3509685 : r3509733;
return r3509734;
}



Bits error versus x
Results
if x < -660441.7377912633 or 905.9497489091251 < x Initial program 58.6
Taylor expanded around inf 0.0
Simplified0.0
if -660441.7377912633 < x < 905.9497489091251Initial 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 2019153
(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))