\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}\;t \le -5.481507593861878460145816755067202944483 \cdot 10^{-263} \lor \neg \left(t \le 6.565726687888966078976606764192118335339 \cdot 10^{-111}\right):\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\sqrt[3]{\ell}}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}}\\
\end{array}double f(double t, double l, double k) {
double r117325 = 2.0;
double r117326 = t;
double r117327 = 3.0;
double r117328 = pow(r117326, r117327);
double r117329 = l;
double r117330 = r117329 * r117329;
double r117331 = r117328 / r117330;
double r117332 = k;
double r117333 = sin(r117332);
double r117334 = r117331 * r117333;
double r117335 = tan(r117332);
double r117336 = r117334 * r117335;
double r117337 = 1.0;
double r117338 = r117332 / r117326;
double r117339 = pow(r117338, r117325);
double r117340 = r117337 + r117339;
double r117341 = r117340 + r117337;
double r117342 = r117336 * r117341;
double r117343 = r117325 / r117342;
return r117343;
}
double f(double t, double l, double k) {
double r117344 = t;
double r117345 = -5.4815075938618785e-263;
bool r117346 = r117344 <= r117345;
double r117347 = 6.565726687888966e-111;
bool r117348 = r117344 <= r117347;
double r117349 = !r117348;
bool r117350 = r117346 || r117349;
double r117351 = 2.0;
double r117352 = cbrt(r117344);
double r117353 = 3.0;
double r117354 = pow(r117352, r117353);
double r117355 = l;
double r117356 = cbrt(r117355);
double r117357 = r117356 * r117356;
double r117358 = r117354 / r117357;
double r117359 = 0.3333333333333333;
double r117360 = r117359 * r117353;
double r117361 = pow(r117344, r117360);
double r117362 = r117361 / r117356;
double r117363 = r117358 * r117362;
double r117364 = r117354 / r117355;
double r117365 = k;
double r117366 = sin(r117365);
double r117367 = r117364 * r117366;
double r117368 = r117363 * r117367;
double r117369 = tan(r117365);
double r117370 = 1.0;
double r117371 = r117365 / r117344;
double r117372 = pow(r117371, r117351);
double r117373 = r117370 + r117372;
double r117374 = r117373 + r117370;
double r117375 = r117369 * r117374;
double r117376 = r117368 * r117375;
double r117377 = r117351 / r117376;
double r117378 = 3.0;
double r117379 = pow(r117344, r117378);
double r117380 = 2.0;
double r117381 = pow(r117366, r117380);
double r117382 = r117379 * r117381;
double r117383 = cos(r117365);
double r117384 = pow(r117355, r117380);
double r117385 = r117383 * r117384;
double r117386 = r117382 / r117385;
double r117387 = r117351 * r117386;
double r117388 = 1.0;
double r117389 = -1.0;
double r117390 = pow(r117389, r117353);
double r117391 = r117388 / r117390;
double r117392 = pow(r117391, r117370);
double r117393 = pow(r117365, r117380);
double r117394 = r117393 * r117381;
double r117395 = r117344 * r117394;
double r117396 = r117395 / r117385;
double r117397 = r117392 * r117396;
double r117398 = r117387 - r117397;
double r117399 = r117351 / r117398;
double r117400 = r117350 ? r117377 : r117399;
return r117400;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -5.4815075938618785e-263 or 6.565726687888966e-111 < t Initial program 27.6
rmApplied add-cube-cbrt27.8
Applied unpow-prod-down27.8
Applied times-frac21.3
Applied associate-*l*19.3
rmApplied add-cube-cbrt19.3
Applied unpow-prod-down19.3
Applied times-frac14.1
rmApplied associate-*l*14.0
rmApplied pow1/340.6
Applied pow-pow13.9
if -5.4815075938618785e-263 < t < 6.565726687888966e-111Initial program 64.0
Taylor expanded around -inf 39.7
Final simplification17.5
herbie shell --seed 2019362
(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))))