\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 -1.483471567911308015546022244990908340338 \cdot 10^{137}:\\
\;\;\;\;2 \cdot \left({\left(\frac{1}{{\left(\sqrt[3]{k} \cdot \sqrt[3]{k}\right)}^{2}}\right)}^{1} \cdot \left({\left(\frac{1}{{\left(\sqrt[3]{k}\right)}^{2}}\right)}^{1} \cdot \left(\frac{\sqrt[3]{\cos k} \cdot \sqrt[3]{\cos k}}{\frac{\sin k}{\ell}} \cdot \left({\left(\frac{1}{{t}^{1}}\right)}^{1} \cdot \frac{\sqrt[3]{\cos k}}{\frac{\sin k}{\ell}}\right)\right)\right)\right)\\
\mathbf{elif}\;k \le -9.730288504549286413960983657697353219034 \cdot 10^{-144}:\\
\;\;\;\;\left(\left({\left(\frac{1}{{t}^{1}}\right)}^{1} \cdot \frac{\sqrt[3]{\cos k}}{\frac{\sin k}{\ell}}\right) \cdot \left({\left(\frac{1}{{k}^{2}}\right)}^{1} \cdot \frac{\sqrt[3]{\cos k} \cdot \sqrt[3]{\cos k}}{\frac{\sin k}{\ell}}\right)\right) \cdot 2\\
\mathbf{elif}\;k \le 6.303690353095152979093645195617976968905 \cdot 10^{-135}:\\
\;\;\;\;\left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left(\frac{\sqrt[3]{\cos k} \cdot \sqrt[3]{\cos k}}{\frac{\sin k}{\ell}} \cdot \left({\left(\frac{1}{{t}^{1}}\right)}^{1} \cdot \frac{\sqrt[3]{\cos k}}{\frac{\sin k}{\ell}}\right)\right)\right) \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1}\right) \cdot 2\\
\mathbf{elif}\;k \le 1.488110455716832872059851006586872843391 \cdot 10^{172}:\\
\;\;\;\;\left(\left({\left(\frac{1}{{t}^{1}}\right)}^{1} \cdot \frac{\sqrt[3]{\cos k}}{\frac{\sin k}{\ell}}\right) \cdot \left({\left(\frac{1}{{k}^{2}}\right)}^{1} \cdot \frac{\sqrt[3]{\cos k} \cdot \sqrt[3]{\cos k}}{\frac{\sin k}{\ell}}\right)\right) \cdot 2\\
\mathbf{else}:\\
\;\;\;\;\left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left(\frac{\sqrt[3]{\cos k} \cdot \sqrt[3]{\cos k}}{\frac{\sin k}{\ell}} \cdot \left({\left(\frac{1}{{t}^{1}}\right)}^{1} \cdot \frac{\sqrt[3]{\cos k}}{\frac{\sin k}{\ell}}\right)\right)\right) \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1}\right) \cdot 2\\
\end{array}double f(double t, double l, double k) {
double r12786763 = 2.0;
double r12786764 = t;
double r12786765 = 3.0;
double r12786766 = pow(r12786764, r12786765);
double r12786767 = l;
double r12786768 = r12786767 * r12786767;
double r12786769 = r12786766 / r12786768;
double r12786770 = k;
double r12786771 = sin(r12786770);
double r12786772 = r12786769 * r12786771;
double r12786773 = tan(r12786770);
double r12786774 = r12786772 * r12786773;
double r12786775 = 1.0;
double r12786776 = r12786770 / r12786764;
double r12786777 = pow(r12786776, r12786763);
double r12786778 = r12786775 + r12786777;
double r12786779 = r12786778 - r12786775;
double r12786780 = r12786774 * r12786779;
double r12786781 = r12786763 / r12786780;
return r12786781;
}
double f(double t, double l, double k) {
double r12786782 = k;
double r12786783 = -1.483471567911308e+137;
bool r12786784 = r12786782 <= r12786783;
double r12786785 = 2.0;
double r12786786 = 1.0;
double r12786787 = cbrt(r12786782);
double r12786788 = r12786787 * r12786787;
double r12786789 = pow(r12786788, r12786785);
double r12786790 = r12786786 / r12786789;
double r12786791 = 1.0;
double r12786792 = pow(r12786790, r12786791);
double r12786793 = pow(r12786787, r12786785);
double r12786794 = r12786786 / r12786793;
double r12786795 = pow(r12786794, r12786791);
double r12786796 = cos(r12786782);
double r12786797 = cbrt(r12786796);
double r12786798 = r12786797 * r12786797;
double r12786799 = sin(r12786782);
double r12786800 = l;
double r12786801 = r12786799 / r12786800;
double r12786802 = r12786798 / r12786801;
double r12786803 = t;
double r12786804 = pow(r12786803, r12786791);
double r12786805 = r12786786 / r12786804;
double r12786806 = pow(r12786805, r12786791);
double r12786807 = r12786797 / r12786801;
double r12786808 = r12786806 * r12786807;
double r12786809 = r12786802 * r12786808;
double r12786810 = r12786795 * r12786809;
double r12786811 = r12786792 * r12786810;
double r12786812 = r12786785 * r12786811;
double r12786813 = -9.730288504549286e-144;
bool r12786814 = r12786782 <= r12786813;
double r12786815 = pow(r12786782, r12786785);
double r12786816 = r12786786 / r12786815;
double r12786817 = pow(r12786816, r12786791);
double r12786818 = r12786817 * r12786802;
double r12786819 = r12786808 * r12786818;
double r12786820 = r12786819 * r12786785;
double r12786821 = 6.303690353095153e-135;
bool r12786822 = r12786782 <= r12786821;
double r12786823 = 2.0;
double r12786824 = r12786785 / r12786823;
double r12786825 = pow(r12786782, r12786824);
double r12786826 = r12786786 / r12786825;
double r12786827 = pow(r12786826, r12786791);
double r12786828 = r12786827 * r12786809;
double r12786829 = r12786828 * r12786827;
double r12786830 = r12786829 * r12786785;
double r12786831 = 1.4881104557168329e+172;
bool r12786832 = r12786782 <= r12786831;
double r12786833 = r12786832 ? r12786820 : r12786830;
double r12786834 = r12786822 ? r12786830 : r12786833;
double r12786835 = r12786814 ? r12786820 : r12786834;
double r12786836 = r12786784 ? r12786812 : r12786835;
return r12786836;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -1.483471567911308e+137Initial program 40.6
Simplified34.1
Taylor expanded around inf 23.5
rmApplied *-un-lft-identity23.5
Applied times-frac23.5
Applied unpow-prod-down23.5
Applied associate-*l*23.7
Simplified23.7
rmApplied add-cube-cbrt23.8
Applied times-frac23.8
Applied associate-*l*22.1
rmApplied add-cube-cbrt22.1
Applied unpow-prod-down22.1
Applied *-un-lft-identity22.1
Applied times-frac21.8
Applied unpow-prod-down21.8
Applied associate-*l*16.9
if -1.483471567911308e+137 < k < -9.730288504549286e-144 or 6.303690353095153e-135 < k < 1.4881104557168329e+172Initial program 53.4
Simplified43.6
Taylor expanded around inf 18.4
rmApplied *-un-lft-identity18.4
Applied times-frac18.3
Applied unpow-prod-down18.3
Applied associate-*l*18.4
Simplified15.9
rmApplied add-cube-cbrt16.0
Applied times-frac15.8
Applied associate-*l*9.0
rmApplied associate-*r*4.4
if -9.730288504549286e-144 < k < 6.303690353095153e-135 or 1.4881104557168329e+172 < k Initial program 44.3
Simplified40.5
Taylor expanded around inf 31.0
rmApplied *-un-lft-identity31.0
Applied times-frac31.2
Applied unpow-prod-down31.2
Applied associate-*l*31.2
Simplified30.1
rmApplied add-cube-cbrt30.1
Applied times-frac30.0
Applied associate-*l*28.9
rmApplied sqr-pow28.9
Applied *-un-lft-identity28.9
Applied times-frac28.9
Applied unpow-prod-down28.9
Applied associate-*l*16.7
Final simplification9.9
herbie shell --seed 2019200
(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))))