\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}\;\ell \le -3.67292000376777888 \cdot 10^{146}:\\
\;\;\;\;\frac{2}{\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}} \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}}\right) \cdot \left(\sqrt[3]{\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)}\\
\mathbf{elif}\;\ell \le -1.54143404308992113 \cdot 10^{-91}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{elif}\;\ell \le 3.85026174332709824 \cdot 10^{-146}:\\
\;\;\;\;\frac{2}{\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 \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\
\mathbf{elif}\;\ell \le 2.14633943998772201 \cdot 10^{102}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)\right)}\\
\end{array}double f(double t, double l, double k) {
double r117952 = 2.0;
double r117953 = t;
double r117954 = 3.0;
double r117955 = pow(r117953, r117954);
double r117956 = l;
double r117957 = r117956 * r117956;
double r117958 = r117955 / r117957;
double r117959 = k;
double r117960 = sin(r117959);
double r117961 = r117958 * r117960;
double r117962 = tan(r117959);
double r117963 = r117961 * r117962;
double r117964 = 1.0;
double r117965 = r117959 / r117953;
double r117966 = pow(r117965, r117952);
double r117967 = r117964 + r117966;
double r117968 = r117967 + r117964;
double r117969 = r117963 * r117968;
double r117970 = r117952 / r117969;
return r117970;
}
double f(double t, double l, double k) {
double r117971 = l;
double r117972 = -3.672920003767779e+146;
bool r117973 = r117971 <= r117972;
double r117974 = 2.0;
double r117975 = t;
double r117976 = cbrt(r117975);
double r117977 = 3.0;
double r117978 = pow(r117976, r117977);
double r117979 = r117978 / r117971;
double r117980 = r117978 * r117979;
double r117981 = cbrt(r117979);
double r117982 = r117981 * r117981;
double r117983 = k;
double r117984 = sin(r117983);
double r117985 = r117981 * r117984;
double r117986 = r117982 * r117985;
double r117987 = r117980 * r117986;
double r117988 = tan(r117983);
double r117989 = r117987 * r117988;
double r117990 = 1.0;
double r117991 = r117983 / r117975;
double r117992 = pow(r117991, r117974);
double r117993 = r117990 + r117992;
double r117994 = r117993 + r117990;
double r117995 = r117989 * r117994;
double r117996 = r117974 / r117995;
double r117997 = -1.541434043089921e-91;
bool r117998 = r117971 <= r117997;
double r117999 = 3.0;
double r118000 = pow(r117975, r117999);
double r118001 = 2.0;
double r118002 = pow(r117984, r118001);
double r118003 = r118000 * r118002;
double r118004 = cos(r117983);
double r118005 = pow(r117971, r118001);
double r118006 = r118004 * r118005;
double r118007 = r118003 / r118006;
double r118008 = r117974 * r118007;
double r118009 = 1.0;
double r118010 = -1.0;
double r118011 = pow(r118010, r117977);
double r118012 = r118009 / r118011;
double r118013 = pow(r118012, r117990);
double r118014 = pow(r117983, r118001);
double r118015 = r118014 * r118002;
double r118016 = r117975 * r118015;
double r118017 = r118016 / r118006;
double r118018 = r118013 * r118017;
double r118019 = r118008 - r118018;
double r118020 = r117974 / r118019;
double r118021 = 3.8502617433270982e-146;
bool r118022 = r117971 <= r118021;
double r118023 = r117979 * r117984;
double r118024 = r117980 * r118023;
double r118025 = cbrt(r118024);
double r118026 = r118025 * r118025;
double r118027 = r118026 * r118025;
double r118028 = r117988 * r117994;
double r118029 = r118027 * r118028;
double r118030 = r117974 / r118029;
double r118031 = 2.146339439987722e+102;
bool r118032 = r117971 <= r118031;
double r118033 = r118023 * r118028;
double r118034 = r117980 * r118033;
double r118035 = r117974 / r118034;
double r118036 = r118032 ? r118020 : r118035;
double r118037 = r118022 ? r118030 : r118036;
double r118038 = r117998 ? r118020 : r118037;
double r118039 = r117973 ? r117996 : r118038;
return r118039;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -3.672920003767779e+146Initial program 62.8
rmApplied add-cube-cbrt62.8
Applied unpow-prod-down62.8
Applied times-frac42.4
Applied associate-*l*42.4
rmApplied *-un-lft-identity42.4
Applied unpow-prod-down42.4
Applied times-frac24.3
Simplified24.3
rmApplied add-cube-cbrt24.3
Applied associate-*l*24.3
if -3.672920003767779e+146 < l < -1.541434043089921e-91 or 3.8502617433270982e-146 < l < 2.146339439987722e+102Initial program 26.3
Taylor expanded around -inf 18.6
if -1.541434043089921e-91 < l < 3.8502617433270982e-146Initial program 24.2
rmApplied add-cube-cbrt24.3
Applied unpow-prod-down24.3
Applied times-frac17.9
Applied associate-*l*15.4
rmApplied *-un-lft-identity15.4
Applied unpow-prod-down15.4
Applied times-frac9.9
Simplified9.9
rmApplied associate-*l*9.4
rmApplied add-cube-cbrt9.4
if 2.146339439987722e+102 < l Initial program 55.4
rmApplied add-cube-cbrt55.5
Applied unpow-prod-down55.5
Applied times-frac39.2
Applied associate-*l*39.2
rmApplied *-un-lft-identity39.2
Applied unpow-prod-down39.2
Applied times-frac25.8
Simplified25.8
rmApplied associate-*l*25.8
rmApplied associate-*l*26.5
Final simplification16.6
herbie shell --seed 2020049
(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))))