\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)}2 \cdot \left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\right)double f(double t, double l, double k) {
double r89369 = 2.0;
double r89370 = t;
double r89371 = 3.0;
double r89372 = pow(r89370, r89371);
double r89373 = l;
double r89374 = r89373 * r89373;
double r89375 = r89372 / r89374;
double r89376 = k;
double r89377 = sin(r89376);
double r89378 = r89375 * r89377;
double r89379 = tan(r89376);
double r89380 = r89378 * r89379;
double r89381 = 1.0;
double r89382 = r89376 / r89370;
double r89383 = pow(r89382, r89369);
double r89384 = r89381 + r89383;
double r89385 = r89384 - r89381;
double r89386 = r89380 * r89385;
double r89387 = r89369 / r89386;
return r89387;
}
double f(double t, double l, double k) {
double r89388 = 2.0;
double r89389 = 1.0;
double r89390 = cbrt(r89389);
double r89391 = r89390 * r89390;
double r89392 = k;
double r89393 = 2.0;
double r89394 = r89388 / r89393;
double r89395 = pow(r89392, r89394);
double r89396 = r89391 / r89395;
double r89397 = 1.0;
double r89398 = pow(r89396, r89397);
double r89399 = t;
double r89400 = pow(r89399, r89397);
double r89401 = r89395 * r89400;
double r89402 = r89390 / r89401;
double r89403 = pow(r89402, r89397);
double r89404 = cos(r89392);
double r89405 = sin(r89392);
double r89406 = fabs(r89405);
double r89407 = r89404 / r89406;
double r89408 = l;
double r89409 = r89406 / r89408;
double r89410 = r89408 / r89409;
double r89411 = r89407 * r89410;
double r89412 = r89403 * r89411;
double r89413 = r89398 * r89412;
double r89414 = r89388 * r89413;
return r89414;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 48.1
Simplified40.3
Taylor expanded around inf 22.5
rmApplied sqr-pow22.5
Applied associate-*l*20.3
rmApplied add-cube-cbrt20.3
Applied times-frac20.1
Applied unpow-prod-down20.1
Applied associate-*l*18.6
rmApplied add-sqr-sqrt18.6
Applied times-frac18.6
Simplified18.6
Simplified16.2
Final simplification16.2
herbie shell --seed 2020064
(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))))