\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 -5.603375851643622422415683198718436485274 \cdot 10^{157} \lor \neg \left(k \le -3.459437993813782241454138708740781760984 \cdot 10^{-153}\right):\\
\;\;\;\;2 \cdot \frac{{\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \left(\frac{\cos k}{\left|\sin k\right|} \cdot \ell\right)}{\frac{\left|\sin k\right|}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\left(\frac{1}{{k}^{2}}\right)}^{1} \cdot \left({\left(\frac{1}{{t}^{1}}\right)}^{1} \cdot \left(\frac{\cos k}{\left|\sin k\right|} \cdot \ell\right)\right)}{\frac{\left|\sin k\right|}{\ell}}\\
\end{array}double f(double t, double l, double k) {
double r74930 = 2.0;
double r74931 = t;
double r74932 = 3.0;
double r74933 = pow(r74931, r74932);
double r74934 = l;
double r74935 = r74934 * r74934;
double r74936 = r74933 / r74935;
double r74937 = k;
double r74938 = sin(r74937);
double r74939 = r74936 * r74938;
double r74940 = tan(r74937);
double r74941 = r74939 * r74940;
double r74942 = 1.0;
double r74943 = r74937 / r74931;
double r74944 = pow(r74943, r74930);
double r74945 = r74942 + r74944;
double r74946 = r74945 - r74942;
double r74947 = r74941 * r74946;
double r74948 = r74930 / r74947;
return r74948;
}
double f(double t, double l, double k) {
double r74949 = k;
double r74950 = -5.603375851643622e+157;
bool r74951 = r74949 <= r74950;
double r74952 = -3.459437993813782e-153;
bool r74953 = r74949 <= r74952;
double r74954 = !r74953;
bool r74955 = r74951 || r74954;
double r74956 = 2.0;
double r74957 = 1.0;
double r74958 = 2.0;
double r74959 = r74956 / r74958;
double r74960 = pow(r74949, r74959);
double r74961 = t;
double r74962 = 1.0;
double r74963 = pow(r74961, r74962);
double r74964 = r74960 * r74963;
double r74965 = r74960 * r74964;
double r74966 = r74957 / r74965;
double r74967 = pow(r74966, r74962);
double r74968 = cos(r74949);
double r74969 = sin(r74949);
double r74970 = fabs(r74969);
double r74971 = r74968 / r74970;
double r74972 = l;
double r74973 = r74971 * r74972;
double r74974 = r74967 * r74973;
double r74975 = r74970 / r74972;
double r74976 = r74974 / r74975;
double r74977 = r74956 * r74976;
double r74978 = pow(r74949, r74956);
double r74979 = r74957 / r74978;
double r74980 = pow(r74979, r74962);
double r74981 = r74957 / r74963;
double r74982 = pow(r74981, r74962);
double r74983 = r74982 * r74973;
double r74984 = r74980 * r74983;
double r74985 = r74984 / r74975;
double r74986 = r74956 * r74985;
double r74987 = r74955 ? r74977 : r74986;
return r74987;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -5.603375851643622e+157 or -3.459437993813782e-153 < k Initial program 46.2
Simplified39.3
Taylor expanded around inf 23.7
rmApplied add-sqr-sqrt23.7
Applied times-frac23.7
Simplified23.7
Simplified22.9
rmApplied associate-*r/22.2
Applied associate-*r/18.7
rmApplied sqr-pow18.7
Applied associate-*l*12.3
if -5.603375851643622e+157 < k < -3.459437993813782e-153Initial program 53.5
Simplified43.6
Taylor expanded around inf 17.7
rmApplied add-sqr-sqrt17.7
Applied times-frac17.7
Simplified17.7
Simplified16.0
rmApplied associate-*r/14.1
Applied associate-*r/8.0
rmApplied *-un-lft-identity8.0
Applied times-frac7.7
Applied unpow-prod-down7.7
Applied associate-*l*4.1
Final simplification10.0
herbie shell --seed 2019291
(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))))