\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{\frac{\frac{2}{\frac{\frac{k}{\ell}}{\frac{1}{t}}}}{\sin k}}{\frac{k}{\ell} \cdot \tan k}double f(double t, double l, double k) {
double r8416127 = 2.0;
double r8416128 = t;
double r8416129 = 3.0;
double r8416130 = pow(r8416128, r8416129);
double r8416131 = l;
double r8416132 = r8416131 * r8416131;
double r8416133 = r8416130 / r8416132;
double r8416134 = k;
double r8416135 = sin(r8416134);
double r8416136 = r8416133 * r8416135;
double r8416137 = tan(r8416134);
double r8416138 = r8416136 * r8416137;
double r8416139 = 1.0;
double r8416140 = r8416134 / r8416128;
double r8416141 = pow(r8416140, r8416127);
double r8416142 = r8416139 + r8416141;
double r8416143 = r8416142 - r8416139;
double r8416144 = r8416138 * r8416143;
double r8416145 = r8416127 / r8416144;
return r8416145;
}
double f(double t, double l, double k) {
double r8416146 = 2.0;
double r8416147 = k;
double r8416148 = l;
double r8416149 = r8416147 / r8416148;
double r8416150 = 1.0;
double r8416151 = t;
double r8416152 = r8416150 / r8416151;
double r8416153 = r8416149 / r8416152;
double r8416154 = r8416146 / r8416153;
double r8416155 = sin(r8416147);
double r8416156 = r8416154 / r8416155;
double r8416157 = tan(r8416147);
double r8416158 = r8416149 * r8416157;
double r8416159 = r8416156 / r8416158;
return r8416159;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 46.9
Simplified26.5
Taylor expanded around inf 14.5
Simplified14.5
rmApplied *-un-lft-identity14.5
Applied div-inv14.6
Applied *-un-lft-identity14.6
Applied times-frac10.6
Applied times-frac3.0
Applied *-un-lft-identity3.0
Applied times-frac2.6
Applied times-frac1.6
Simplified1.6
rmApplied associate-*l/1.6
Applied associate-/l/1.3
Final simplification1.3
herbie shell --seed 2019132
(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))))