\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{-2}{-\sin k} \cdot \frac{\frac{1}{\frac{\sqrt[3]{t} \cdot \sqrt[3]{t}}{\frac{\ell}{k}}}}{\frac{\tan k}{\frac{1}{\frac{\sqrt[3]{t}}{\frac{\ell}{k}}}}}double f(double t, double l, double k) {
double r2661163 = 2.0;
double r2661164 = t;
double r2661165 = 3.0;
double r2661166 = pow(r2661164, r2661165);
double r2661167 = l;
double r2661168 = r2661167 * r2661167;
double r2661169 = r2661166 / r2661168;
double r2661170 = k;
double r2661171 = sin(r2661170);
double r2661172 = r2661169 * r2661171;
double r2661173 = tan(r2661170);
double r2661174 = r2661172 * r2661173;
double r2661175 = 1.0;
double r2661176 = r2661170 / r2661164;
double r2661177 = pow(r2661176, r2661163);
double r2661178 = r2661175 + r2661177;
double r2661179 = r2661178 - r2661175;
double r2661180 = r2661174 * r2661179;
double r2661181 = r2661163 / r2661180;
return r2661181;
}
double f(double t, double l, double k) {
double r2661182 = -2.0;
double r2661183 = k;
double r2661184 = sin(r2661183);
double r2661185 = -r2661184;
double r2661186 = r2661182 / r2661185;
double r2661187 = 1.0;
double r2661188 = t;
double r2661189 = cbrt(r2661188);
double r2661190 = r2661189 * r2661189;
double r2661191 = l;
double r2661192 = r2661191 / r2661183;
double r2661193 = r2661190 / r2661192;
double r2661194 = r2661187 / r2661193;
double r2661195 = tan(r2661183);
double r2661196 = r2661189 / r2661192;
double r2661197 = r2661187 / r2661196;
double r2661198 = r2661195 / r2661197;
double r2661199 = r2661194 / r2661198;
double r2661200 = r2661186 * r2661199;
return r2661200;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
Initial program 47.2
Simplified21.1
rmApplied frac-2neg21.1
Simplified16.7
rmApplied distribute-lft-neg-in16.7
Applied div-inv16.7
Applied times-frac16.5
Simplified14.8
Taylor expanded around inf 21.6
Simplified8.3
rmApplied add-cube-cbrt8.6
Applied times-frac3.2
Applied *-un-lft-identity3.2
Applied times-frac2.8
Applied associate-/l*1.4
Final simplification1.4
herbie shell --seed 2019143 +o rules:numerics
(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))))