\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 -4.979647532007629 \cdot 10^{+167}:\\
\;\;\;\;\frac{1}{\frac{t}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{\sqrt[3]{k} \cdot \sqrt[3]{k}}}} \cdot \left(\frac{\frac{\sqrt{2} \cdot \ell}{\tan k}}{\sin k} \cdot \frac{\sqrt{2}}{\sqrt[3]{k} \cdot \frac{k}{\sqrt[3]{\ell}}}\right)\\
\mathbf{elif}\;k \le 1.6846040995988706 \cdot 10^{+107}:\\
\;\;\;\;\left(\frac{\sqrt{\sqrt{2}}}{\frac{k}{\frac{\ell}{k}}} \cdot \frac{\sqrt{\sqrt{2}}}{t}\right) \cdot \left(\frac{\sqrt{2}}{\tan k} \cdot \frac{\ell}{\sin k}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{t}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{\sqrt[3]{k} \cdot \sqrt[3]{k}}}} \cdot \left(\frac{\frac{\sqrt{2} \cdot \ell}{\tan k}}{\sin k} \cdot \frac{\sqrt{2}}{\sqrt[3]{k} \cdot \frac{k}{\sqrt[3]{\ell}}}\right)\\
\end{array}double f(double t, double l, double k) {
double r4266903 = 2.0;
double r4266904 = t;
double r4266905 = 3.0;
double r4266906 = pow(r4266904, r4266905);
double r4266907 = l;
double r4266908 = r4266907 * r4266907;
double r4266909 = r4266906 / r4266908;
double r4266910 = k;
double r4266911 = sin(r4266910);
double r4266912 = r4266909 * r4266911;
double r4266913 = tan(r4266910);
double r4266914 = r4266912 * r4266913;
double r4266915 = 1.0;
double r4266916 = r4266910 / r4266904;
double r4266917 = pow(r4266916, r4266903);
double r4266918 = r4266915 + r4266917;
double r4266919 = r4266918 - r4266915;
double r4266920 = r4266914 * r4266919;
double r4266921 = r4266903 / r4266920;
return r4266921;
}
double f(double t, double l, double k) {
double r4266922 = k;
double r4266923 = -4.979647532007629e+167;
bool r4266924 = r4266922 <= r4266923;
double r4266925 = 1.0;
double r4266926 = t;
double r4266927 = l;
double r4266928 = cbrt(r4266927);
double r4266929 = r4266928 * r4266928;
double r4266930 = cbrt(r4266922);
double r4266931 = r4266930 * r4266930;
double r4266932 = r4266929 / r4266931;
double r4266933 = r4266926 / r4266932;
double r4266934 = r4266925 / r4266933;
double r4266935 = 2.0;
double r4266936 = sqrt(r4266935);
double r4266937 = r4266936 * r4266927;
double r4266938 = tan(r4266922);
double r4266939 = r4266937 / r4266938;
double r4266940 = sin(r4266922);
double r4266941 = r4266939 / r4266940;
double r4266942 = r4266922 / r4266928;
double r4266943 = r4266930 * r4266942;
double r4266944 = r4266936 / r4266943;
double r4266945 = r4266941 * r4266944;
double r4266946 = r4266934 * r4266945;
double r4266947 = 1.6846040995988706e+107;
bool r4266948 = r4266922 <= r4266947;
double r4266949 = sqrt(r4266936);
double r4266950 = r4266927 / r4266922;
double r4266951 = r4266922 / r4266950;
double r4266952 = r4266949 / r4266951;
double r4266953 = r4266949 / r4266926;
double r4266954 = r4266952 * r4266953;
double r4266955 = r4266936 / r4266938;
double r4266956 = r4266927 / r4266940;
double r4266957 = r4266955 * r4266956;
double r4266958 = r4266954 * r4266957;
double r4266959 = r4266948 ? r4266958 : r4266946;
double r4266960 = r4266924 ? r4266946 : r4266959;
return r4266960;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -4.979647532007629e+167 or 1.6846040995988706e+107 < k Initial program 37.9
Simplified19.9
rmApplied div-inv19.9
Applied add-sqr-sqrt19.9
Applied times-frac19.9
Applied times-frac20.2
Simplified7.2
Simplified7.2
rmApplied add-sqr-sqrt7.2
Applied add-cube-cbrt7.3
Applied times-frac7.3
Applied add-cube-cbrt7.4
Applied times-frac7.4
Applied times-frac8.6
Applied *-un-lft-identity8.6
Applied times-frac8.2
Applied associate-*l*4.3
Simplified4.3
if -4.979647532007629e+167 < k < 1.6846040995988706e+107Initial program 53.4
Simplified13.9
rmApplied div-inv13.9
Applied add-sqr-sqrt14.1
Applied times-frac13.9
Applied times-frac10.9
Simplified6.9
Simplified3.6
rmApplied add-cube-cbrt3.6
Applied *-un-lft-identity3.6
Applied times-frac3.6
Applied *-un-lft-identity3.6
Applied times-frac3.6
Applied times-frac2.1
Applied add-sqr-sqrt2.1
Applied sqrt-prod2.0
Applied times-frac2.0
Simplified2.0
Simplified2.0
Final simplification3.0
herbie shell --seed 2019152
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))))