\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 -1.69166056710356561 \cdot 10^{109}:\\
\;\;\;\;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 -7.3151552301839124 \cdot 10^{-79}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k}{\sin k}\right) \cdot \ell}{\frac{\sin k}{\ell}}\\
\mathbf{elif}\;k \le 4.33536330724835694 \cdot 10^{-148}:\\
\;\;\;\;2 \cdot \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 \left(\left(\frac{\cos k}{\sin k} \cdot \ell\right) \cdot \frac{1}{\frac{\sin k}{\ell}}\right)\right)\\
\mathbf{elif}\;k \le 5.87932956341524449 \cdot 10^{115}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k}{\sin k}\right) \cdot \ell}{\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 \left(\frac{\cos k}{\sin k} \cdot \frac{\ell}{\frac{\sin k}{\ell}}\right)\right)\right)\\
\end{array}double code(double t, double l, double k) {
return ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) pow(t, 3.0)) / ((double) (l * l)))) * ((double) sin(k)))) * ((double) tan(k)))) * ((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) - 1.0))))));
}
double code(double t, double l, double k) {
double VAR;
if ((k <= -1.6916605671035656e+109)) {
VAR = ((double) (2.0 * ((double) (((double) pow(((double) (1.0 / ((double) pow(k, ((double) (2.0 / 2.0)))))), 1.0)) * ((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, ((double) (2.0 / 2.0)))) * ((double) pow(t, 1.0)))))), 1.0)) * ((double) (((double) (((double) cos(k)) / ((double) sin(k)))) * ((double) (l / ((double) (((double) sin(k)) / l))))))))))));
} else {
double VAR_1;
if ((k <= -7.315155230183912e-79)) {
VAR_1 = ((double) (2.0 * ((double) (((double) (((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, 2.0)) * ((double) pow(t, 1.0)))))), 1.0)) * ((double) (((double) cos(k)) / ((double) sin(k)))))) * l)) / ((double) (((double) sin(k)) / l))))));
} else {
double VAR_2;
if ((k <= 4.335363307248357e-148)) {
VAR_2 = ((double) (2.0 * ((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, ((double) (2.0 / 2.0)))) * ((double) (((double) pow(k, ((double) (2.0 / 2.0)))) * ((double) pow(t, 1.0)))))))), 1.0)) * ((double) (((double) (((double) (((double) cos(k)) / ((double) sin(k)))) * l)) * ((double) (1.0 / ((double) (((double) sin(k)) / l))))))))));
} else {
double VAR_3;
if ((k <= 5.8793295634152445e+115)) {
VAR_3 = ((double) (2.0 * ((double) (((double) (((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, 2.0)) * ((double) pow(t, 1.0)))))), 1.0)) * ((double) (((double) cos(k)) / ((double) sin(k)))))) * l)) / ((double) (((double) sin(k)) / l))))));
} else {
VAR_3 = ((double) (2.0 * ((double) (((double) pow(((double) (1.0 / ((double) pow(k, ((double) (2.0 / 2.0)))))), 1.0)) * ((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, ((double) (2.0 / 2.0)))) * ((double) pow(t, 1.0)))))), 1.0)) * ((double) (((double) (((double) cos(k)) / ((double) sin(k)))) * ((double) (l / ((double) (((double) sin(k)) / l))))))))))));
}
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 < -1.6916605671035656e+109 or 5.8793295634152445e+115 < k Initial program 40.7
Simplified34.6
Taylor expanded around inf 22.4
rmApplied sqr-pow22.4
Applied associate-*l*18.4
rmApplied add-sqr-sqrt41.3
Applied unpow-prod-down41.3
Applied times-frac41.3
Simplified41.3
Simplified18.4
rmApplied *-un-lft-identity18.4
Applied times-frac18.2
Applied unpow-prod-down18.2
Applied associate-*l*15.5
if -1.6916605671035656e+109 < k < -7.315155230183912e-79 or 4.335363307248357e-148 < k < 5.8793295634152445e+115Initial program 54.4
Simplified43.2
Taylor expanded around inf 16.4
rmApplied sqr-pow16.4
Applied associate-*l*16.4
rmApplied add-sqr-sqrt39.6
Applied unpow-prod-down39.6
Applied times-frac39.6
Simplified39.5
Simplified15.2
rmApplied associate-*r/14.0
Applied associate-*r/6.7
Simplified8.5
if -7.315155230183912e-79 < k < 4.335363307248357e-148Initial program 63.9
Simplified63.6
Taylor expanded around inf 47.6
rmApplied sqr-pow47.6
Applied associate-*l*47.6
rmApplied add-sqr-sqrt63.1
Applied unpow-prod-down63.1
Applied times-frac61.2
Simplified61.2
Simplified27.8
rmApplied div-inv27.8
Applied associate-*r*17.7
Final simplification12.8
herbie shell --seed 2020120
(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))))