\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 -6.2391434554590811 \cdot 10^{149}:\\
\;\;\;\;2 \cdot \left({\left(\frac{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{elif}\;k \le -5.15729782667803904 \cdot 10^{-154}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\
\mathbf{elif}\;k \le 3.3090433695573168 \cdot 10^{-155}:\\
\;\;\;\;2 \cdot \left({\left(\frac{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 \left(\frac{\cos k}{\sin k} \cdot \frac{\ell}{\frac{\sin k}{\ell}}\right)\right)\right)\\
\mathbf{elif}\;k \le 1.22181510276113542 \cdot 10^{154}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left({\left(\frac{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)\\
\end{array}double code(double t, double l, double k) {
return (2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0)));
}
double code(double t, double l, double k) {
double VAR;
if ((k <= -6.239143455459081e+149)) {
VAR = (2.0 * (pow((1.0 / pow(k, (2.0 / 2.0))), 1.0) * (pow((1.0 / (pow(k, (2.0 / 2.0)) * pow(t, 1.0))), 1.0) * ((cos(k) * pow(l, 2.0)) / pow(sin(k), 2.0)))));
} else {
double VAR_1;
if ((k <= -5.157297826678039e-154)) {
VAR_1 = (2.0 * (((pow((1.0 / (pow(k, 2.0) * pow(t, 1.0))), 1.0) * l) * cos(k)) / (sin(k) * (sin(k) / l))));
} else {
double VAR_2;
if ((k <= 3.3090433695573168e-155)) {
VAR_2 = (2.0 * (pow((1.0 / pow(k, (2.0 / 2.0))), 1.0) * (pow((1.0 / (pow(k, (2.0 / 2.0)) * pow(t, 1.0))), 1.0) * ((cos(k) / sin(k)) * (l / (sin(k) / l))))));
} else {
double VAR_3;
if ((k <= 1.2218151027611354e+154)) {
VAR_3 = (2.0 * (((pow((1.0 / (pow(k, 2.0) * pow(t, 1.0))), 1.0) * l) * cos(k)) / (sin(k) * (sin(k) / l))));
} else {
VAR_3 = (2.0 * (pow((1.0 / pow(k, (2.0 / 2.0))), 1.0) * (pow((1.0 / (pow(k, (2.0 / 2.0)) * pow(t, 1.0))), 1.0) * ((cos(k) * pow(l, 2.0)) / pow(sin(k), 2.0)))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -6.239143455459081e+149 or 1.2218151027611354e+154 < k Initial program 38.4
Simplified33.8
Taylor expanded around inf 23.9
rmApplied sqr-pow23.9
Applied associate-*l*18.5
rmApplied *-un-lft-identity18.5
Applied times-frac18.2
Applied unpow-prod-down18.2
Applied associate-*l*14.9
if -6.239143455459081e+149 < k < -5.157297826678039e-154 or 3.3090433695573168e-155 < k < 1.2218151027611354e+154Initial program 53.4
Simplified43.8
Taylor expanded around inf 17.9
rmApplied sqr-pow17.9
Applied associate-*l*17.9
rmApplied add-sqr-sqrt40.6
Applied unpow-prod-down40.6
Applied times-frac40.6
Simplified40.6
Simplified16.5
rmApplied frac-times14.7
Applied associate-*r/8.2
Simplified8.2
if -5.157297826678039e-154 < k < 3.3090433695573168e-155Initial program 64.0
Simplified64.0
Taylor expanded around inf 62.0
rmApplied sqr-pow62.0
Applied associate-*l*61.9
rmApplied add-sqr-sqrt63.0
Applied unpow-prod-down63.0
Applied times-frac59.3
Simplified59.3
Simplified29.0
rmApplied *-un-lft-identity29.0
Applied times-frac29.1
Applied unpow-prod-down29.1
Applied associate-*l*23.0
Final simplification11.4
herbie shell --seed 2020100
(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))))