\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}\;t \le -1.777700292473165462498818001107613277449 \cdot 10^{101}:\\
\;\;\;\;\frac{2}{\frac{\left(\sqrt[3]{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \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)} \cdot \sqrt[3]{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \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)}\right) \cdot \sqrt[3]{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \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)}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}}\\
\mathbf{elif}\;t \le -1.496151205830202533286098915686898280991 \cdot 10^{-95}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{1}^{3}}{\ell} \cdot \left(\frac{{t}^{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}\;t \le 1.065572922882695808628296268030592849626 \cdot 10^{-98}:\\
\;\;\;\;\frac{2}{\frac{\frac{{k}^{2} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot \ell} + 2 \cdot \frac{{t}^{2} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot \ell}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \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) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}}\\
\end{array}double f(double t, double l, double k) {
double r132844 = 2.0;
double r132845 = t;
double r132846 = 3.0;
double r132847 = pow(r132845, r132846);
double r132848 = l;
double r132849 = r132848 * r132848;
double r132850 = r132847 / r132849;
double r132851 = k;
double r132852 = sin(r132851);
double r132853 = r132850 * r132852;
double r132854 = tan(r132851);
double r132855 = r132853 * r132854;
double r132856 = 1.0;
double r132857 = r132851 / r132845;
double r132858 = pow(r132857, r132844);
double r132859 = r132856 + r132858;
double r132860 = r132859 + r132856;
double r132861 = r132855 * r132860;
double r132862 = r132844 / r132861;
return r132862;
}
double f(double t, double l, double k) {
double r132863 = t;
double r132864 = -1.7777002924731655e+101;
bool r132865 = r132863 <= r132864;
double r132866 = 2.0;
double r132867 = cbrt(r132863);
double r132868 = 3.0;
double r132869 = pow(r132867, r132868);
double r132870 = l;
double r132871 = r132869 / r132870;
double r132872 = k;
double r132873 = sin(r132872);
double r132874 = r132871 * r132873;
double r132875 = r132869 * r132874;
double r132876 = tan(r132872);
double r132877 = r132875 * r132876;
double r132878 = 1.0;
double r132879 = r132872 / r132863;
double r132880 = pow(r132879, r132866);
double r132881 = r132878 + r132880;
double r132882 = r132881 + r132878;
double r132883 = r132877 * r132882;
double r132884 = cbrt(r132883);
double r132885 = r132884 * r132884;
double r132886 = r132885 * r132884;
double r132887 = r132870 / r132869;
double r132888 = r132886 / r132887;
double r132889 = r132866 / r132888;
double r132890 = -1.4961512058302025e-95;
bool r132891 = r132863 <= r132890;
double r132892 = 1.0;
double r132893 = pow(r132892, r132868);
double r132894 = r132893 / r132870;
double r132895 = pow(r132863, r132868);
double r132896 = r132895 / r132870;
double r132897 = r132896 * r132873;
double r132898 = r132894 * r132897;
double r132899 = r132898 * r132876;
double r132900 = r132899 * r132882;
double r132901 = r132866 / r132900;
double r132902 = 1.0655729228826958e-98;
bool r132903 = r132863 <= r132902;
double r132904 = 2.0;
double r132905 = pow(r132872, r132904);
double r132906 = pow(r132873, r132904);
double r132907 = r132905 * r132906;
double r132908 = cos(r132872);
double r132909 = r132908 * r132870;
double r132910 = r132907 / r132909;
double r132911 = pow(r132863, r132904);
double r132912 = r132911 * r132906;
double r132913 = r132912 / r132909;
double r132914 = r132866 * r132913;
double r132915 = r132910 + r132914;
double r132916 = r132915 / r132887;
double r132917 = r132866 / r132916;
double r132918 = cbrt(r132873);
double r132919 = r132918 * r132918;
double r132920 = r132871 * r132919;
double r132921 = r132920 * r132918;
double r132922 = r132869 * r132921;
double r132923 = r132922 * r132876;
double r132924 = r132923 * r132882;
double r132925 = r132924 / r132887;
double r132926 = r132866 / r132925;
double r132927 = r132903 ? r132917 : r132926;
double r132928 = r132891 ? r132901 : r132927;
double r132929 = r132865 ? r132889 : r132928;
return r132929;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -1.7777002924731655e+101Initial program 23.2
rmApplied add-cube-cbrt23.2
Applied unpow-prod-down23.2
Applied times-frac15.4
Applied associate-*l*14.8
rmApplied unpow-prod-down14.8
Applied associate-/l*7.4
rmApplied associate-*l/5.4
Applied associate-*l/1.6
Applied associate-*l/1.8
rmApplied add-cube-cbrt1.8
if -1.7777002924731655e+101 < t < -1.4961512058302025e-95Initial program 21.6
rmApplied *-un-lft-identity21.6
Applied unpow-prod-down21.6
Applied times-frac17.1
Applied associate-*l*13.4
if -1.4961512058302025e-95 < t < 1.0655729228826958e-98Initial program 62.1
rmApplied add-cube-cbrt62.1
Applied unpow-prod-down62.1
Applied times-frac52.9
Applied associate-*l*52.9
rmApplied unpow-prod-down52.9
Applied associate-/l*43.8
rmApplied associate-*l/43.8
Applied associate-*l/45.1
Applied associate-*l/41.8
Taylor expanded around inf 23.2
if 1.0655729228826958e-98 < t Initial program 23.6
rmApplied add-cube-cbrt23.7
Applied unpow-prod-down23.7
Applied times-frac17.9
Applied associate-*l*15.6
rmApplied unpow-prod-down15.6
Applied associate-/l*10.5
rmApplied associate-*l/9.3
Applied associate-*l/7.3
Applied associate-*l/6.2
rmApplied add-cube-cbrt6.3
Applied associate-*r*6.3
Final simplification10.3
herbie shell --seed 2019297
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
:precision binary64
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (+ (+ 1 (pow (/ k t) 2)) 1))))