\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 -3.56310395914880465 \cdot 10^{-34} \lor \neg \left(t \le 6.3604377920357265 \cdot 10^{-52}\right):\\
\;\;\;\;\frac{2}{\frac{\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)\right) \cdot \sqrt[3]{\sin k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\frac{{k}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell} + 2 \cdot \frac{{t}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\
\end{array}double f(double t, double l, double k) {
double r107228 = 2.0;
double r107229 = t;
double r107230 = 3.0;
double r107231 = pow(r107229, r107230);
double r107232 = l;
double r107233 = r107232 * r107232;
double r107234 = r107231 / r107233;
double r107235 = k;
double r107236 = sin(r107235);
double r107237 = r107234 * r107236;
double r107238 = tan(r107235);
double r107239 = r107237 * r107238;
double r107240 = 1.0;
double r107241 = r107235 / r107229;
double r107242 = pow(r107241, r107228);
double r107243 = r107240 + r107242;
double r107244 = r107243 + r107240;
double r107245 = r107239 * r107244;
double r107246 = r107228 / r107245;
return r107246;
}
double f(double t, double l, double k) {
double r107247 = t;
double r107248 = -3.5631039591488047e-34;
bool r107249 = r107247 <= r107248;
double r107250 = 6.3604377920357265e-52;
bool r107251 = r107247 <= r107250;
double r107252 = !r107251;
bool r107253 = r107249 || r107252;
double r107254 = 2.0;
double r107255 = cbrt(r107247);
double r107256 = 3.0;
double r107257 = pow(r107255, r107256);
double r107258 = l;
double r107259 = r107257 / r107258;
double r107260 = k;
double r107261 = sin(r107260);
double r107262 = r107259 * r107261;
double r107263 = r107257 * r107262;
double r107264 = cbrt(r107261);
double r107265 = r107264 * r107264;
double r107266 = r107263 * r107265;
double r107267 = r107266 * r107264;
double r107268 = 1.0;
double r107269 = r107260 / r107247;
double r107270 = pow(r107269, r107254);
double r107271 = r107268 + r107270;
double r107272 = r107271 + r107268;
double r107273 = r107267 * r107272;
double r107274 = r107258 / r107257;
double r107275 = cos(r107260);
double r107276 = r107274 * r107275;
double r107277 = r107273 / r107276;
double r107278 = r107254 / r107277;
double r107279 = 2.0;
double r107280 = pow(r107260, r107279);
double r107281 = pow(r107261, r107279);
double r107282 = r107280 * r107281;
double r107283 = r107282 / r107258;
double r107284 = pow(r107247, r107279);
double r107285 = r107284 * r107281;
double r107286 = r107285 / r107258;
double r107287 = r107254 * r107286;
double r107288 = r107283 + r107287;
double r107289 = r107288 / r107276;
double r107290 = r107254 / r107289;
double r107291 = r107253 ? r107278 : r107290;
return r107291;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -3.5631039591488047e-34 or 6.3604377920357265e-52 < t Initial program 21.9
rmApplied add-cube-cbrt22.1
Applied unpow-prod-down22.1
Applied times-frac15.8
Applied associate-*l*13.7
rmApplied unpow-prod-down13.7
Applied associate-/l*8.6
rmApplied tan-quot8.6
Applied associate-*l/7.4
Applied frac-times5.1
Applied associate-*l/4.7
rmApplied add-cube-cbrt4.7
Applied associate-*r*4.7
if -3.5631039591488047e-34 < t < 6.3604377920357265e-52Initial program 54.8
rmApplied add-cube-cbrt54.8
Applied unpow-prod-down54.8
Applied times-frac46.5
Applied associate-*l*45.6
rmApplied unpow-prod-down45.6
Applied associate-/l*38.7
rmApplied tan-quot38.7
Applied associate-*l/38.7
Applied frac-times40.2
Applied associate-*l/36.4
Taylor expanded around inf 22.5
Final simplification10.1
herbie shell --seed 2020047
(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))))