\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 -1.15058253029693869 \cdot 10^{-125}:\\
\;\;\;\;\frac{2}{\left(\left(\sqrt[3]{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k} \cdot \sqrt[3]{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k}\right) \cdot \sqrt[3]{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;t \le 2.70311402462102548 \cdot 10^{-151}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}} \cdot \left(\frac{k \cdot t}{\ell} - \frac{1}{6} \cdot \frac{{k}^{3} \cdot t}{\ell}\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}}\right) \cdot \left(\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) \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\end{array}double f(double t, double l, double k) {
double r124276 = 2.0;
double r124277 = t;
double r124278 = 3.0;
double r124279 = pow(r124277, r124278);
double r124280 = l;
double r124281 = r124280 * r124280;
double r124282 = r124279 / r124281;
double r124283 = k;
double r124284 = sin(r124283);
double r124285 = r124282 * r124284;
double r124286 = tan(r124283);
double r124287 = r124285 * r124286;
double r124288 = 1.0;
double r124289 = r124283 / r124277;
double r124290 = pow(r124289, r124276);
double r124291 = r124288 + r124290;
double r124292 = r124291 + r124288;
double r124293 = r124287 * r124292;
double r124294 = r124276 / r124293;
return r124294;
}
double f(double t, double l, double k) {
double r124295 = t;
double r124296 = -1.1505825302969387e-125;
bool r124297 = r124295 <= r124296;
double r124298 = 2.0;
double r124299 = cbrt(r124295);
double r124300 = r124299 * r124299;
double r124301 = 3.0;
double r124302 = 2.0;
double r124303 = r124301 / r124302;
double r124304 = pow(r124300, r124303);
double r124305 = l;
double r124306 = cbrt(r124305);
double r124307 = r124306 * r124306;
double r124308 = r124304 / r124307;
double r124309 = r124304 / r124306;
double r124310 = pow(r124299, r124301);
double r124311 = r124310 / r124305;
double r124312 = k;
double r124313 = sin(r124312);
double r124314 = r124311 * r124313;
double r124315 = r124309 * r124314;
double r124316 = r124308 * r124315;
double r124317 = tan(r124312);
double r124318 = r124316 * r124317;
double r124319 = cbrt(r124318);
double r124320 = r124319 * r124319;
double r124321 = r124320 * r124319;
double r124322 = 1.0;
double r124323 = r124312 / r124295;
double r124324 = pow(r124323, r124298);
double r124325 = r124322 + r124324;
double r124326 = r124325 + r124322;
double r124327 = r124321 * r124326;
double r124328 = r124298 / r124327;
double r124329 = 2.7031140246210255e-151;
bool r124330 = r124295 <= r124329;
double r124331 = r124312 * r124295;
double r124332 = r124331 / r124305;
double r124333 = 0.16666666666666666;
double r124334 = 3.0;
double r124335 = pow(r124312, r124334);
double r124336 = r124335 * r124295;
double r124337 = r124336 / r124305;
double r124338 = r124333 * r124337;
double r124339 = r124332 - r124338;
double r124340 = r124309 * r124339;
double r124341 = r124308 * r124340;
double r124342 = r124341 * r124317;
double r124343 = r124342 * r124326;
double r124344 = r124298 / r124343;
double r124345 = r124308 * r124309;
double r124346 = cbrt(r124314);
double r124347 = r124346 * r124346;
double r124348 = r124347 * r124346;
double r124349 = r124345 * r124348;
double r124350 = r124349 * r124317;
double r124351 = r124350 * r124326;
double r124352 = r124298 / r124351;
double r124353 = r124330 ? r124344 : r124352;
double r124354 = r124297 ? r124328 : r124353;
return r124354;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -1.1505825302969387e-125Initial program 24.7
rmApplied add-cube-cbrt24.9
Applied unpow-prod-down24.9
Applied times-frac17.4
Applied associate-*l*15.2
rmApplied add-cube-cbrt15.2
Applied sqr-pow15.2
Applied times-frac10.7
rmApplied associate-*l*9.8
rmApplied add-cube-cbrt9.8
if -1.1505825302969387e-125 < t < 2.7031140246210255e-151Initial program 64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac59.2
Applied associate-*l*59.2
rmApplied add-cube-cbrt59.2
Applied sqr-pow59.2
Applied times-frac47.8
rmApplied associate-*l*47.8
Taylor expanded around 0 40.4
if 2.7031140246210255e-151 < t Initial program 27.5
rmApplied add-cube-cbrt27.7
Applied unpow-prod-down27.7
Applied times-frac19.2
Applied associate-*l*17.2
rmApplied add-cube-cbrt17.2
Applied sqr-pow17.2
Applied times-frac12.3
rmApplied add-cube-cbrt12.4
Final simplification16.5
herbie shell --seed 2020020
(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))))