\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 -8.573355553736766686780551655893902286167 \cdot 10^{70}:\\
\;\;\;\;\frac{2}{\frac{\left(\sin k \cdot \left(\left(\sqrt[3]{\sin k} \cdot \left(\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right) \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right)\right) \cdot {\left(\sqrt[3]{t}\right)}^{3}\right)\right) \cdot \left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}{\cos k \cdot \frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}}\\
\mathbf{elif}\;t \le 1.020078012341945963530160021779867815533 \cdot 10^{-104}:\\
\;\;\;\;\frac{2}{\frac{t}{\frac{\cos k}{\sin k \cdot \sin k} \cdot \left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right)} + \left(\frac{t \cdot \left(t \cdot t\right)}{\ell} \cdot \frac{\sin k \cdot \sin k}{\ell}\right) \cdot \frac{2}{\cos k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left(\sin k \cdot \left(\left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right)\right) \cdot {\left(\sqrt[3]{t}\right)}^{3}\right)\right) \cdot \left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}{\cos k \cdot \frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}}\\
\end{array}double f(double t, double l, double k) {
double r3956498 = 2.0;
double r3956499 = t;
double r3956500 = 3.0;
double r3956501 = pow(r3956499, r3956500);
double r3956502 = l;
double r3956503 = r3956502 * r3956502;
double r3956504 = r3956501 / r3956503;
double r3956505 = k;
double r3956506 = sin(r3956505);
double r3956507 = r3956504 * r3956506;
double r3956508 = tan(r3956505);
double r3956509 = r3956507 * r3956508;
double r3956510 = 1.0;
double r3956511 = r3956505 / r3956499;
double r3956512 = pow(r3956511, r3956498);
double r3956513 = r3956510 + r3956512;
double r3956514 = r3956513 + r3956510;
double r3956515 = r3956509 * r3956514;
double r3956516 = r3956498 / r3956515;
return r3956516;
}
double f(double t, double l, double k) {
double r3956517 = t;
double r3956518 = -8.573355553736767e+70;
bool r3956519 = r3956517 <= r3956518;
double r3956520 = 2.0;
double r3956521 = k;
double r3956522 = sin(r3956521);
double r3956523 = cbrt(r3956522);
double r3956524 = r3956523 * r3956523;
double r3956525 = cbrt(r3956517);
double r3956526 = 3.0;
double r3956527 = pow(r3956525, r3956526);
double r3956528 = l;
double r3956529 = r3956527 / r3956528;
double r3956530 = r3956524 * r3956529;
double r3956531 = r3956523 * r3956530;
double r3956532 = r3956531 * r3956527;
double r3956533 = r3956522 * r3956532;
double r3956534 = 1.0;
double r3956535 = r3956521 / r3956517;
double r3956536 = pow(r3956535, r3956520);
double r3956537 = r3956534 + r3956536;
double r3956538 = r3956534 + r3956537;
double r3956539 = r3956533 * r3956538;
double r3956540 = cos(r3956521);
double r3956541 = r3956528 / r3956527;
double r3956542 = r3956540 * r3956541;
double r3956543 = r3956539 / r3956542;
double r3956544 = r3956520 / r3956543;
double r3956545 = 1.020078012341946e-104;
bool r3956546 = r3956517 <= r3956545;
double r3956547 = r3956522 * r3956522;
double r3956548 = r3956540 / r3956547;
double r3956549 = r3956528 / r3956521;
double r3956550 = r3956549 * r3956549;
double r3956551 = r3956548 * r3956550;
double r3956552 = r3956517 / r3956551;
double r3956553 = r3956517 * r3956517;
double r3956554 = r3956517 * r3956553;
double r3956555 = r3956554 / r3956528;
double r3956556 = r3956547 / r3956528;
double r3956557 = r3956555 * r3956556;
double r3956558 = r3956520 / r3956540;
double r3956559 = r3956557 * r3956558;
double r3956560 = r3956552 + r3956559;
double r3956561 = r3956520 / r3956560;
double r3956562 = r3956529 * r3956522;
double r3956563 = cbrt(r3956562);
double r3956564 = r3956563 * r3956563;
double r3956565 = r3956563 * r3956564;
double r3956566 = r3956565 * r3956527;
double r3956567 = r3956522 * r3956566;
double r3956568 = r3956567 * r3956538;
double r3956569 = r3956568 / r3956542;
double r3956570 = r3956520 / r3956569;
double r3956571 = r3956546 ? r3956561 : r3956570;
double r3956572 = r3956519 ? r3956544 : r3956571;
return r3956572;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -8.573355553736767e+70Initial program 24.7
rmApplied add-cube-cbrt24.8
Applied unpow-prod-down24.8
Applied times-frac17.0
Applied associate-*l*15.4
rmApplied unpow-prod-down15.4
Applied associate-/l*7.4
rmApplied tan-quot7.4
Applied associate-*l/5.5
Applied frac-times2.2
Applied associate-*l/2.3
rmApplied add-cube-cbrt2.3
Applied associate-*r*2.3
if -8.573355553736767e+70 < t < 1.020078012341946e-104Initial program 48.6
rmApplied add-cube-cbrt48.7
Applied unpow-prod-down48.7
Applied times-frac41.5
Applied associate-*l*40.0
rmApplied unpow-prod-down40.0
Applied associate-/l*34.1
Taylor expanded around inf 34.0
Simplified11.8
if 1.020078012341946e-104 < t Initial program 23.6
rmApplied add-cube-cbrt23.8
Applied unpow-prod-down23.8
Applied times-frac17.6
Applied associate-*l*15.0
rmApplied unpow-prod-down15.0
Applied associate-/l*10.6
rmApplied tan-quot10.6
Applied associate-*l/9.6
Applied frac-times7.9
Applied associate-*l/6.8
rmApplied add-cube-cbrt6.8
Final simplification7.5
herbie shell --seed 2019174
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))