\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}\;\ell \le -1.0143648348712567 \cdot 10^{+22}:\\
\;\;\;\;\frac{\frac{2}{\frac{k}{t} \cdot \frac{k}{t} + 2}}{\frac{\frac{\sin k}{\frac{\ell}{t}}}{\cos k}} \cdot \frac{\frac{\ell}{\sin k \cdot t}}{t}\\
\mathbf{elif}\;\ell \le 8.796141127304192 \cdot 10^{+67}:\\
\;\;\;\;\frac{\frac{\frac{\sqrt{2}}{t}}{\frac{\sin k}{\ell} \cdot \left(t \cdot \sqrt[3]{\frac{k}{t} \cdot \frac{k}{t} + 2}\right)}}{\frac{\sin k}{\ell} \cdot \left(t \cdot \sqrt[3]{\frac{k}{t} \cdot \frac{k}{t} + 2}\right)} \cdot \left(\cos k \cdot \frac{\sqrt{2}}{\sqrt[3]{\frac{k}{t} \cdot \frac{k}{t} + 2}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\frac{k}{t} \cdot \frac{k}{t} + 2}}{\frac{\frac{\sin k}{\frac{\ell}{t}}}{\cos k}} \cdot \frac{\frac{\ell}{\sin k \cdot t}}{t}\\
\end{array}double f(double t, double l, double k) {
double r36794598 = 2.0;
double r36794599 = t;
double r36794600 = 3.0;
double r36794601 = pow(r36794599, r36794600);
double r36794602 = l;
double r36794603 = r36794602 * r36794602;
double r36794604 = r36794601 / r36794603;
double r36794605 = k;
double r36794606 = sin(r36794605);
double r36794607 = r36794604 * r36794606;
double r36794608 = tan(r36794605);
double r36794609 = r36794607 * r36794608;
double r36794610 = 1.0;
double r36794611 = r36794605 / r36794599;
double r36794612 = pow(r36794611, r36794598);
double r36794613 = r36794610 + r36794612;
double r36794614 = r36794613 + r36794610;
double r36794615 = r36794609 * r36794614;
double r36794616 = r36794598 / r36794615;
return r36794616;
}
double f(double t, double l, double k) {
double r36794617 = l;
double r36794618 = -1.0143648348712567e+22;
bool r36794619 = r36794617 <= r36794618;
double r36794620 = 2.0;
double r36794621 = k;
double r36794622 = t;
double r36794623 = r36794621 / r36794622;
double r36794624 = r36794623 * r36794623;
double r36794625 = r36794624 + r36794620;
double r36794626 = r36794620 / r36794625;
double r36794627 = sin(r36794621);
double r36794628 = r36794617 / r36794622;
double r36794629 = r36794627 / r36794628;
double r36794630 = cos(r36794621);
double r36794631 = r36794629 / r36794630;
double r36794632 = r36794626 / r36794631;
double r36794633 = r36794627 * r36794622;
double r36794634 = r36794617 / r36794633;
double r36794635 = r36794634 / r36794622;
double r36794636 = r36794632 * r36794635;
double r36794637 = 8.796141127304192e+67;
bool r36794638 = r36794617 <= r36794637;
double r36794639 = sqrt(r36794620);
double r36794640 = r36794639 / r36794622;
double r36794641 = r36794627 / r36794617;
double r36794642 = cbrt(r36794625);
double r36794643 = r36794622 * r36794642;
double r36794644 = r36794641 * r36794643;
double r36794645 = r36794640 / r36794644;
double r36794646 = r36794645 / r36794644;
double r36794647 = r36794639 / r36794642;
double r36794648 = r36794630 * r36794647;
double r36794649 = r36794646 * r36794648;
double r36794650 = r36794638 ? r36794649 : r36794636;
double r36794651 = r36794619 ? r36794636 : r36794650;
return r36794651;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -1.0143648348712567e+22 or 8.796141127304192e+67 < l Initial program 47.8
Simplified29.3
rmApplied tan-quot29.3
Applied associate-*r/29.3
Applied associate-*l/29.3
Simplified21.0
rmApplied *-un-lft-identity21.0
Applied times-frac21.0
Applied *-un-lft-identity21.0
Applied *-un-lft-identity21.0
Applied times-frac21.0
Applied times-frac19.0
Simplified18.9
if -1.0143648348712567e+22 < l < 8.796141127304192e+67Initial program 23.7
Simplified21.7
rmApplied tan-quot21.7
Applied associate-*r/21.7
Applied associate-*l/21.7
Simplified10.9
rmApplied div-inv10.9
Applied add-cube-cbrt11.1
Applied add-sqr-sqrt11.0
Applied times-frac11.1
Applied times-frac11.1
Simplified7.9
Simplified7.9
Final simplification11.5
herbie shell --seed 2019121
(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))))