\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.910065564714718288643928397762008915017 \cdot 10^{183} \lor \neg \left(\ell \le 1.147874348448085091862458118841413831719 \cdot 10^{154}\right):\\
\;\;\;\;\frac{\frac{2}{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\frac{{\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \cos k}{\sin k} \cdot \frac{{\ell}^{2}}{\sin k}\right) \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1}\right)\\
\end{array}double f(double t, double l, double k) {
double r92677 = 2.0;
double r92678 = t;
double r92679 = 3.0;
double r92680 = pow(r92678, r92679);
double r92681 = l;
double r92682 = r92681 * r92681;
double r92683 = r92680 / r92682;
double r92684 = k;
double r92685 = sin(r92684);
double r92686 = r92683 * r92685;
double r92687 = tan(r92684);
double r92688 = r92686 * r92687;
double r92689 = 1.0;
double r92690 = r92684 / r92678;
double r92691 = pow(r92690, r92677);
double r92692 = r92689 + r92691;
double r92693 = r92692 - r92689;
double r92694 = r92688 * r92693;
double r92695 = r92677 / r92694;
return r92695;
}
double f(double t, double l, double k) {
double r92696 = l;
double r92697 = -1.9100655647147183e+183;
bool r92698 = r92696 <= r92697;
double r92699 = 1.1478743484480851e+154;
bool r92700 = r92696 <= r92699;
double r92701 = !r92700;
bool r92702 = r92698 || r92701;
double r92703 = 2.0;
double r92704 = t;
double r92705 = cbrt(r92704);
double r92706 = r92705 * r92705;
double r92707 = 3.0;
double r92708 = pow(r92706, r92707);
double r92709 = r92708 / r92696;
double r92710 = pow(r92705, r92707);
double r92711 = r92710 / r92696;
double r92712 = k;
double r92713 = sin(r92712);
double r92714 = r92711 * r92713;
double r92715 = r92709 * r92714;
double r92716 = tan(r92712);
double r92717 = r92715 * r92716;
double r92718 = r92703 / r92717;
double r92719 = r92712 / r92704;
double r92720 = pow(r92719, r92703);
double r92721 = r92718 / r92720;
double r92722 = 1.0;
double r92723 = 2.0;
double r92724 = r92703 / r92723;
double r92725 = pow(r92712, r92724);
double r92726 = 1.0;
double r92727 = pow(r92704, r92726);
double r92728 = r92725 * r92727;
double r92729 = r92722 / r92728;
double r92730 = pow(r92729, r92726);
double r92731 = cos(r92712);
double r92732 = r92730 * r92731;
double r92733 = r92732 / r92713;
double r92734 = pow(r92696, r92723);
double r92735 = r92734 / r92713;
double r92736 = r92733 * r92735;
double r92737 = r92722 / r92725;
double r92738 = pow(r92737, r92726);
double r92739 = r92736 * r92738;
double r92740 = r92703 * r92739;
double r92741 = r92702 ? r92721 : r92740;
return r92741;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -1.9100655647147183e+183 or 1.1478743484480851e+154 < l Initial program 64.0
Simplified64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac50.1
Applied associate-*l*50.1
if -1.9100655647147183e+183 < l < 1.1478743484480851e+154Initial program 45.4
Simplified36.8
Taylor expanded around inf 15.5
rmApplied sqr-pow15.5
Applied associate-*l*13.0
rmApplied add-sqr-sqrt13.0
Applied times-frac12.9
Applied unpow-prod-down12.9
Applied associate-*l*11.3
Simplified11.3
rmApplied sqr-pow11.3
Applied times-frac10.9
Applied associate-*r*10.9
Simplified10.8
Final simplification16.3
herbie shell --seed 2019235 +o rules:numerics
(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))))