\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\begin{array}{l}
\mathbf{if}\;t \le -7.7186020624409743 \cdot 10^{211}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right)}\\
\mathbf{elif}\;t \le -4.4828637577847008 \cdot 10^{49}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right)}\\
\mathbf{elif}\;t \le 2.85659231871301485 \cdot 10^{41}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r209723 = 2.0;
double r209724 = n;
double r209725 = r209723 * r209724;
double r209726 = U;
double r209727 = r209725 * r209726;
double r209728 = t;
double r209729 = l;
double r209730 = r209729 * r209729;
double r209731 = Om;
double r209732 = r209730 / r209731;
double r209733 = r209723 * r209732;
double r209734 = r209728 - r209733;
double r209735 = r209729 / r209731;
double r209736 = pow(r209735, r209723);
double r209737 = r209724 * r209736;
double r209738 = U_;
double r209739 = r209726 - r209738;
double r209740 = r209737 * r209739;
double r209741 = r209734 - r209740;
double r209742 = r209727 * r209741;
double r209743 = sqrt(r209742);
return r209743;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r209744 = t;
double r209745 = -7.718602062440974e+211;
bool r209746 = r209744 <= r209745;
double r209747 = 2.0;
double r209748 = n;
double r209749 = r209747 * r209748;
double r209750 = U;
double r209751 = r209749 * r209750;
double r209752 = l;
double r209753 = Om;
double r209754 = r209752 / r209753;
double r209755 = r209752 * r209754;
double r209756 = r209747 * r209755;
double r209757 = r209744 - r209756;
double r209758 = r209751 * r209757;
double r209759 = sqrt(r209758);
double r209760 = -4.482863757784701e+49;
bool r209761 = r209744 <= r209760;
double r209762 = pow(r209754, r209747);
double r209763 = r209748 * r209762;
double r209764 = U_;
double r209765 = r209750 - r209764;
double r209766 = r209763 * r209765;
double r209767 = r209757 - r209766;
double r209768 = r209750 * r209767;
double r209769 = r209749 * r209768;
double r209770 = sqrt(r209769);
double r209771 = 2.856592318713015e+41;
bool r209772 = r209744 <= r209771;
double r209773 = 2.0;
double r209774 = r209747 / r209773;
double r209775 = pow(r209754, r209774);
double r209776 = r209748 * r209775;
double r209777 = r209775 * r209765;
double r209778 = r209776 * r209777;
double r209779 = r209757 - r209778;
double r209780 = r209751 * r209779;
double r209781 = sqrt(r209780);
double r209782 = sqrt(r209751);
double r209783 = r209773 * r209774;
double r209784 = pow(r209754, r209783);
double r209785 = r209748 * r209784;
double r209786 = r209765 * r209785;
double r209787 = r209757 - r209786;
double r209788 = sqrt(r209787);
double r209789 = r209782 * r209788;
double r209790 = r209772 ? r209781 : r209789;
double r209791 = r209761 ? r209770 : r209790;
double r209792 = r209746 ? r209759 : r209791;
return r209792;
}



Bits error versus n



Bits error versus U



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus U*
Results
if t < -7.718602062440974e+211Initial program 40.3
rmApplied *-un-lft-identity40.3
Applied times-frac37.1
Simplified37.1
Taylor expanded around 0 36.3
if -7.718602062440974e+211 < t < -4.482863757784701e+49Initial program 33.8
rmApplied *-un-lft-identity33.8
Applied times-frac30.7
Simplified30.7
rmApplied associate-*l*31.0
if -4.482863757784701e+49 < t < 2.856592318713015e+41Initial program 33.4
rmApplied *-un-lft-identity33.4
Applied times-frac30.7
Simplified30.7
rmApplied sqr-pow30.7
Applied associate-*r*29.4
rmApplied associate-*l*29.1
if 2.856592318713015e+41 < t Initial program 35.2
rmApplied *-un-lft-identity35.2
Applied times-frac32.3
Simplified32.3
rmApplied sqr-pow32.3
Applied associate-*r*31.9
rmApplied sqrt-prod24.4
Simplified24.6
Final simplification28.9
herbie shell --seed 2020046
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
:precision binary64
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))