\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}:\\
\;\;\;\;\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left(\frac{\cos k}{\sin k} \cdot \left({\left(\frac{\frac{1}{{k}^{\left(\frac{2}{2}\right)}}}{{t}^{1}}\right)}^{1} \cdot \frac{{\ell}^{2}}{\sin k}\right)\right)\right) \cdot 2\\
\end{array}double f(double t, double l, double k) {
double r82645 = 2.0;
double r82646 = t;
double r82647 = 3.0;
double r82648 = pow(r82646, r82647);
double r82649 = l;
double r82650 = r82649 * r82649;
double r82651 = r82648 / r82650;
double r82652 = k;
double r82653 = sin(r82652);
double r82654 = r82651 * r82653;
double r82655 = tan(r82652);
double r82656 = r82654 * r82655;
double r82657 = 1.0;
double r82658 = r82652 / r82646;
double r82659 = pow(r82658, r82645);
double r82660 = r82657 + r82659;
double r82661 = r82660 - r82657;
double r82662 = r82656 * r82661;
double r82663 = r82645 / r82662;
return r82663;
}
double f(double t, double l, double k) {
double r82664 = l;
double r82665 = -1.9100655647147183e+183;
bool r82666 = r82664 <= r82665;
double r82667 = 1.1478743484480851e+154;
bool r82668 = r82664 <= r82667;
double r82669 = !r82668;
bool r82670 = r82666 || r82669;
double r82671 = 2.0;
double r82672 = t;
double r82673 = cbrt(r82672);
double r82674 = r82673 * r82673;
double r82675 = 3.0;
double r82676 = pow(r82674, r82675);
double r82677 = r82676 / r82664;
double r82678 = pow(r82673, r82675);
double r82679 = r82678 / r82664;
double r82680 = k;
double r82681 = sin(r82680);
double r82682 = r82679 * r82681;
double r82683 = r82677 * r82682;
double r82684 = tan(r82680);
double r82685 = r82683 * r82684;
double r82686 = r82671 / r82685;
double r82687 = r82680 / r82672;
double r82688 = pow(r82687, r82671);
double r82689 = r82686 / r82688;
double r82690 = 1.0;
double r82691 = 2.0;
double r82692 = r82671 / r82691;
double r82693 = pow(r82680, r82692);
double r82694 = r82690 / r82693;
double r82695 = 1.0;
double r82696 = pow(r82694, r82695);
double r82697 = cos(r82680);
double r82698 = r82697 / r82681;
double r82699 = pow(r82672, r82695);
double r82700 = r82694 / r82699;
double r82701 = pow(r82700, r82695);
double r82702 = pow(r82664, r82691);
double r82703 = r82702 / r82681;
double r82704 = r82701 * r82703;
double r82705 = r82698 * r82704;
double r82706 = r82696 * r82705;
double r82707 = r82706 * r82671;
double r82708 = r82670 ? r82689 : r82707;
return r82708;
}



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 unpow211.3
Applied times-frac10.9
Applied associate-*l*10.8
Simplified10.7
Final simplification16.2
herbie shell --seed 2019235
(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))))