\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)}\frac{\left(\frac{1}{k} \cdot \frac{t}{t \cdot \frac{\sin k}{\ell}}\right) \cdot \frac{2 \cdot \frac{\ell}{k}}{t}}{\tan k}double f(double t, double l, double k) {
double r12019414 = 2.0;
double r12019415 = t;
double r12019416 = 3.0;
double r12019417 = pow(r12019415, r12019416);
double r12019418 = l;
double r12019419 = r12019418 * r12019418;
double r12019420 = r12019417 / r12019419;
double r12019421 = k;
double r12019422 = sin(r12019421);
double r12019423 = r12019420 * r12019422;
double r12019424 = tan(r12019421);
double r12019425 = r12019423 * r12019424;
double r12019426 = 1.0;
double r12019427 = r12019421 / r12019415;
double r12019428 = pow(r12019427, r12019414);
double r12019429 = r12019426 + r12019428;
double r12019430 = r12019429 - r12019426;
double r12019431 = r12019425 * r12019430;
double r12019432 = r12019414 / r12019431;
return r12019432;
}
double f(double t, double l, double k) {
double r12019433 = 1.0;
double r12019434 = k;
double r12019435 = r12019433 / r12019434;
double r12019436 = t;
double r12019437 = sin(r12019434);
double r12019438 = l;
double r12019439 = r12019437 / r12019438;
double r12019440 = r12019436 * r12019439;
double r12019441 = r12019436 / r12019440;
double r12019442 = r12019435 * r12019441;
double r12019443 = 2.0;
double r12019444 = r12019438 / r12019434;
double r12019445 = r12019443 * r12019444;
double r12019446 = r12019445 / r12019436;
double r12019447 = r12019442 * r12019446;
double r12019448 = tan(r12019434);
double r12019449 = r12019447 / r12019448;
return r12019449;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 46.8
Simplified29.4
rmApplied *-un-lft-identity29.4
Applied times-frac29.3
Applied times-frac17.7
rmApplied associate-/r/17.7
Applied associate-/r*10.9
rmApplied *-un-lft-identity10.9
Applied *-un-lft-identity10.9
Applied times-frac10.9
Applied *-un-lft-identity10.9
Applied times-frac10.9
Applied div-inv10.9
Applied add-sqr-sqrt10.9
Applied times-frac10.9
Applied times-frac10.9
Simplified10.9
Simplified10.8
Taylor expanded around inf 4.2
Final simplification4.2
herbie shell --seed 2019135
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))))