\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\begin{array}{l}
\mathbf{if}\;k \le -7.601645767137617971798374690994089117302 \cdot 10^{149}:\\
\;\;\;\;\frac{2}{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)\right) \cdot \sqrt[3]{\sin k}\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le -5.044518950321330373642746174830737878481 \cdot 10^{47}:\\
\;\;\;\;\frac{2}{{\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \left(\frac{\sin k \cdot \sin k}{\frac{\left(\ell \cdot \ell\right) \cdot \cos k}{\left(k \cdot k\right) \cdot \left(t \cdot {\left(\sqrt[3]{-1}\right)}^{9}\right)}} + 2 \cdot \left(\frac{\left(t \cdot t\right) \cdot t}{\ell \cdot \ell} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left(\sin k \cdot \sin k\right)}{\cos k}\right)\right)}\\
\mathbf{elif}\;k \le 5.839305045058099752517812933614005683714 \cdot 10^{-206}:\\
\;\;\;\;\frac{2}{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right)\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le 1.136972446704892211624900722840210240971 \cdot 10^{-128}:\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\ell}}\\
\mathbf{elif}\;k \le 5.149546980068362951212817718174660270701 \cdot 10^{83}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\sqrt[3]{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)} \cdot \sqrt[3]{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)}\right) \cdot \sqrt[3]{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)}\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le 8.960378091421642333491931412091851428398 \cdot 10^{176}:\\
\;\;\;\;\frac{2}{\mathsf{fma}\left(\frac{t \cdot \left(t \cdot t\right)}{\frac{\left(\ell \cdot \ell\right) \cdot \cos k}{\sin k \cdot \sin k}}, 2, \frac{\frac{\left(t \cdot \left(\sin k \cdot \sin k\right)\right) \cdot \left(k \cdot k\right)}{\cos k}}{\ell \cdot \ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\frac{{\left(\sqrt[3]{\sqrt[3]{t} \cdot \sqrt[3]{t}}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{\sqrt[3]{t}}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \sin k\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\end{array}double f(double t, double l, double k) {
double r4908727 = 2.0;
double r4908728 = t;
double r4908729 = 3.0;
double r4908730 = pow(r4908728, r4908729);
double r4908731 = l;
double r4908732 = r4908731 * r4908731;
double r4908733 = r4908730 / r4908732;
double r4908734 = k;
double r4908735 = sin(r4908734);
double r4908736 = r4908733 * r4908735;
double r4908737 = tan(r4908734);
double r4908738 = r4908736 * r4908737;
double r4908739 = 1.0;
double r4908740 = r4908734 / r4908728;
double r4908741 = pow(r4908740, r4908727);
double r4908742 = r4908739 + r4908741;
double r4908743 = r4908742 + r4908739;
double r4908744 = r4908738 * r4908743;
double r4908745 = r4908727 / r4908744;
return r4908745;
}
double f(double t, double l, double k) {
double r4908746 = k;
double r4908747 = -7.601645767137618e+149;
bool r4908748 = r4908746 <= r4908747;
double r4908749 = 2.0;
double r4908750 = t;
double r4908751 = cbrt(r4908750);
double r4908752 = 3.0;
double r4908753 = pow(r4908751, r4908752);
double r4908754 = l;
double r4908755 = r4908753 / r4908754;
double r4908756 = sin(r4908746);
double r4908757 = cbrt(r4908756);
double r4908758 = r4908757 * r4908757;
double r4908759 = r4908755 * r4908758;
double r4908760 = r4908759 * r4908757;
double r4908761 = r4908755 * r4908760;
double r4908762 = r4908753 * r4908761;
double r4908763 = tan(r4908746);
double r4908764 = r4908762 * r4908763;
double r4908765 = 1.0;
double r4908766 = r4908746 / r4908750;
double r4908767 = pow(r4908766, r4908749);
double r4908768 = r4908765 + r4908767;
double r4908769 = r4908768 + r4908765;
double r4908770 = r4908764 * r4908769;
double r4908771 = r4908749 / r4908770;
double r4908772 = -5.04451895032133e+47;
bool r4908773 = r4908746 <= r4908772;
double r4908774 = 1.0;
double r4908775 = -1.0;
double r4908776 = pow(r4908775, r4908752);
double r4908777 = r4908774 / r4908776;
double r4908778 = pow(r4908777, r4908765);
double r4908779 = r4908756 * r4908756;
double r4908780 = r4908754 * r4908754;
double r4908781 = cos(r4908746);
double r4908782 = r4908780 * r4908781;
double r4908783 = r4908746 * r4908746;
double r4908784 = cbrt(r4908775);
double r4908785 = 9.0;
double r4908786 = pow(r4908784, r4908785);
double r4908787 = r4908750 * r4908786;
double r4908788 = r4908783 * r4908787;
double r4908789 = r4908782 / r4908788;
double r4908790 = r4908779 / r4908789;
double r4908791 = r4908750 * r4908750;
double r4908792 = r4908791 * r4908750;
double r4908793 = r4908792 / r4908780;
double r4908794 = r4908786 * r4908779;
double r4908795 = r4908794 / r4908781;
double r4908796 = r4908793 * r4908795;
double r4908797 = r4908749 * r4908796;
double r4908798 = r4908790 + r4908797;
double r4908799 = r4908778 * r4908798;
double r4908800 = r4908749 / r4908799;
double r4908801 = 5.8393050450581e-206;
bool r4908802 = r4908746 <= r4908801;
double r4908803 = r4908755 * r4908756;
double r4908804 = r4908755 * r4908803;
double r4908805 = r4908804 * r4908763;
double r4908806 = r4908753 * r4908805;
double r4908807 = r4908806 * r4908769;
double r4908808 = r4908749 / r4908807;
double r4908809 = 1.1369724467048922e-128;
bool r4908810 = r4908746 <= r4908809;
double r4908811 = r4908753 * r4908803;
double r4908812 = r4908753 * r4908811;
double r4908813 = r4908812 * r4908763;
double r4908814 = r4908813 * r4908769;
double r4908815 = r4908814 / r4908754;
double r4908816 = r4908749 / r4908815;
double r4908817 = 5.149546980068363e+83;
bool r4908818 = r4908746 <= r4908817;
double r4908819 = r4908753 * r4908755;
double r4908820 = r4908819 * r4908803;
double r4908821 = cbrt(r4908820);
double r4908822 = r4908821 * r4908821;
double r4908823 = r4908822 * r4908821;
double r4908824 = r4908823 * r4908763;
double r4908825 = r4908824 * r4908769;
double r4908826 = r4908749 / r4908825;
double r4908827 = 8.960378091421642e+176;
bool r4908828 = r4908746 <= r4908827;
double r4908829 = r4908750 * r4908791;
double r4908830 = r4908782 / r4908779;
double r4908831 = r4908829 / r4908830;
double r4908832 = r4908750 * r4908779;
double r4908833 = r4908832 * r4908783;
double r4908834 = r4908833 / r4908781;
double r4908835 = r4908834 / r4908780;
double r4908836 = fma(r4908831, r4908749, r4908835);
double r4908837 = r4908749 / r4908836;
double r4908838 = r4908751 * r4908751;
double r4908839 = cbrt(r4908838);
double r4908840 = pow(r4908839, r4908752);
double r4908841 = cbrt(r4908754);
double r4908842 = r4908841 * r4908841;
double r4908843 = r4908840 / r4908842;
double r4908844 = cbrt(r4908751);
double r4908845 = pow(r4908844, r4908752);
double r4908846 = r4908845 / r4908841;
double r4908847 = r4908846 * r4908756;
double r4908848 = r4908843 * r4908847;
double r4908849 = r4908819 * r4908848;
double r4908850 = r4908849 * r4908763;
double r4908851 = r4908850 * r4908769;
double r4908852 = r4908749 / r4908851;
double r4908853 = r4908828 ? r4908837 : r4908852;
double r4908854 = r4908818 ? r4908826 : r4908853;
double r4908855 = r4908810 ? r4908816 : r4908854;
double r4908856 = r4908802 ? r4908808 : r4908855;
double r4908857 = r4908773 ? r4908800 : r4908856;
double r4908858 = r4908748 ? r4908771 : r4908857;
return r4908858;
}



Bits error versus t



Bits error versus l



Bits error versus k
if k < -7.601645767137618e+149Initial program 34.5
rmApplied add-cube-cbrt34.5
Applied unpow-prod-down34.5
Applied times-frac28.6
Applied associate-*l*28.6
rmApplied *-un-lft-identity28.6
Applied unpow-prod-down28.6
Applied times-frac22.8
Simplified22.8
rmApplied associate-*l*23.0
rmApplied add-cube-cbrt23.0
Applied associate-*r*23.0
if -7.601645767137618e+149 < k < -5.04451895032133e+47Initial program 31.5
rmApplied add-cube-cbrt31.7
Applied unpow-prod-down31.6
Applied times-frac23.6
Applied associate-*l*23.6
Taylor expanded around -inf 19.3
Simplified19.3
if -5.04451895032133e+47 < k < 5.8393050450581e-206Initial program 33.4
rmApplied add-cube-cbrt33.7
Applied unpow-prod-down33.7
Applied times-frac27.7
Applied associate-*l*21.8
rmApplied *-un-lft-identity21.8
Applied unpow-prod-down21.8
Applied times-frac16.2
Simplified16.2
rmApplied associate-*l*13.8
rmApplied associate-*l*13.4
if 5.8393050450581e-206 < k < 1.1369724467048922e-128Initial program 33.1
rmApplied add-cube-cbrt33.3
Applied unpow-prod-down33.3
Applied times-frac25.8
Applied associate-*l*20.2
rmApplied *-un-lft-identity20.2
Applied unpow-prod-down20.2
Applied times-frac12.4
Simplified12.4
rmApplied associate-*l*10.8
rmApplied associate-*l/10.8
Applied associate-*r/16.0
Applied associate-*l/15.8
Applied associate-*l/15.8
if 1.1369724467048922e-128 < k < 5.149546980068363e+83Initial program 30.1
rmApplied add-cube-cbrt30.4
Applied unpow-prod-down30.4
Applied times-frac21.2
Applied associate-*l*20.5
rmApplied *-un-lft-identity20.5
Applied unpow-prod-down20.5
Applied times-frac14.7
Simplified14.7
rmApplied add-cube-cbrt14.7
if 5.149546980068363e+83 < k < 8.960378091421642e+176Initial program 33.4
Taylor expanded around inf 21.8
Simplified21.9
if 8.960378091421642e+176 < k Initial program 35.0
rmApplied add-cube-cbrt35.0
Applied unpow-prod-down35.0
Applied times-frac28.4
Applied associate-*l*28.4
rmApplied *-un-lft-identity28.4
Applied unpow-prod-down28.4
Applied times-frac22.7
Simplified22.7
rmApplied add-cube-cbrt22.7
Applied add-cube-cbrt22.7
Applied cbrt-prod22.7
Applied unpow-prod-down22.7
Applied times-frac22.7
Applied associate-*l*22.7
Final simplification18.2
herbie shell --seed 2019174 +o rules:numerics
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))