\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 -3.499363714053866281078841467858531447298 \cdot 10^{202}:\\
\;\;\;\;2 \cdot \left({\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{{1}^{2}}\right) \cdot \frac{\ell}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\\
\end{array}double f(double t, double l, double k) {
double r365897 = 2.0;
double r365898 = t;
double r365899 = 3.0;
double r365900 = pow(r365898, r365899);
double r365901 = l;
double r365902 = r365901 * r365901;
double r365903 = r365900 / r365902;
double r365904 = k;
double r365905 = sin(r365904);
double r365906 = r365903 * r365905;
double r365907 = tan(r365904);
double r365908 = r365906 * r365907;
double r365909 = 1.0;
double r365910 = r365904 / r365898;
double r365911 = pow(r365910, r365897);
double r365912 = r365909 + r365911;
double r365913 = r365912 - r365909;
double r365914 = r365908 * r365913;
double r365915 = r365897 / r365914;
return r365915;
}
double f(double t, double l, double k) {
double r365916 = k;
double r365917 = -3.4993637140538663e+202;
bool r365918 = r365916 <= r365917;
double r365919 = 2.0;
double r365920 = 1.0;
double r365921 = sqrt(r365920);
double r365922 = 2.0;
double r365923 = r365919 / r365922;
double r365924 = pow(r365916, r365923);
double r365925 = r365921 / r365924;
double r365926 = 1.0;
double r365927 = pow(r365925, r365926);
double r365928 = t;
double r365929 = pow(r365928, r365926);
double r365930 = r365924 * r365929;
double r365931 = r365920 / r365930;
double r365932 = pow(r365931, r365926);
double r365933 = cos(r365916);
double r365934 = l;
double r365935 = pow(r365934, r365922);
double r365936 = r365933 * r365935;
double r365937 = sin(r365916);
double r365938 = pow(r365937, r365922);
double r365939 = r365936 / r365938;
double r365940 = r365932 * r365939;
double r365941 = r365927 * r365940;
double r365942 = r365919 * r365941;
double r365943 = r365924 * r365930;
double r365944 = r365920 / r365943;
double r365945 = pow(r365944, r365926);
double r365946 = cbrt(r365937);
double r365947 = 4.0;
double r365948 = pow(r365946, r365947);
double r365949 = r365948 / r365934;
double r365950 = r365933 / r365949;
double r365951 = pow(r365920, r365922);
double r365952 = r365950 / r365951;
double r365953 = r365945 * r365952;
double r365954 = pow(r365946, r365922);
double r365955 = r365934 / r365954;
double r365956 = r365953 * r365955;
double r365957 = r365919 * r365956;
double r365958 = r365918 ? r365942 : r365957;
return r365958;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -3.4993637140538663e+202Initial program 36.6
Simplified33.0
Taylor expanded around inf 23.0
rmApplied sqr-pow23.0
Applied associate-*l*19.8
rmApplied add-sqr-sqrt19.8
Applied times-frac19.5
Applied unpow-prod-down19.5
Applied associate-*l*16.0
Simplified16.0
if -3.4993637140538663e+202 < k Initial program 50.4
Simplified41.8
Taylor expanded around inf 22.0
rmApplied sqr-pow22.0
Applied associate-*l*20.0
rmApplied add-cube-cbrt20.4
Applied unpow-prod-down20.4
Applied associate-/r*20.1
Simplified17.9
rmApplied *-un-lft-identity17.9
Applied unpow-prod-down17.9
Applied associate-/r/17.7
Applied times-frac16.8
Applied associate-*r*11.6
Final simplification12.2
herbie shell --seed 2019322
(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))))