\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 -1.1607177078091258 \cdot 10^{142}:\\
\;\;\;\;\frac{2}{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \left(\left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right) \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;t \le -4.5548020353125647 \cdot 10^{-113}:\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\ell \cdot \cos k}}\\
\mathbf{elif}\;t \le 6.62256655119128422 \cdot 10^{-131}:\\
\;\;\;\;\frac{2}{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \left(\frac{k \cdot t}{\ell} - \frac{1}{6} \cdot \frac{{k}^{3} \cdot t}{\ell}\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \left(\left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right) \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\end{array}double f(double t, double l, double k) {
double r244 = 2.0;
double r245 = t;
double r246 = 3.0;
double r247 = pow(r245, r246);
double r248 = l;
double r249 = r248 * r248;
double r250 = r247 / r249;
double r251 = k;
double r252 = sin(r251);
double r253 = r250 * r252;
double r254 = tan(r251);
double r255 = r253 * r254;
double r256 = 1.0;
double r257 = r251 / r245;
double r258 = pow(r257, r244);
double r259 = r256 + r258;
double r260 = r259 + r256;
double r261 = r255 * r260;
double r262 = r244 / r261;
return r262;
}
double f(double t, double l, double k) {
double r263 = t;
double r264 = -1.1607177078091258e+142;
bool r265 = r263 <= r264;
double r266 = 2.0;
double r267 = cbrt(r263);
double r268 = r267 * r267;
double r269 = 3.0;
double r270 = 2.0;
double r271 = r269 / r270;
double r272 = pow(r268, r271);
double r273 = l;
double r274 = r272 / r273;
double r275 = pow(r267, r269);
double r276 = r275 / r273;
double r277 = k;
double r278 = sin(r277);
double r279 = r276 * r278;
double r280 = cbrt(r279);
double r281 = r280 * r280;
double r282 = r281 * r280;
double r283 = r274 * r282;
double r284 = r272 * r283;
double r285 = tan(r277);
double r286 = r284 * r285;
double r287 = 1.0;
double r288 = r277 / r263;
double r289 = pow(r288, r266);
double r290 = r287 + r289;
double r291 = r290 + r287;
double r292 = r286 * r291;
double r293 = r266 / r292;
double r294 = -4.5548020353125647e-113;
bool r295 = r263 <= r294;
double r296 = pow(r268, r269);
double r297 = r296 * r279;
double r298 = r297 * r278;
double r299 = r298 * r291;
double r300 = cos(r277);
double r301 = r273 * r300;
double r302 = r299 / r301;
double r303 = r266 / r302;
double r304 = 6.622566551191284e-131;
bool r305 = r263 <= r304;
double r306 = r277 * r263;
double r307 = r306 / r273;
double r308 = 0.16666666666666666;
double r309 = 3.0;
double r310 = pow(r277, r309);
double r311 = r310 * r263;
double r312 = r311 / r273;
double r313 = r308 * r312;
double r314 = r307 - r313;
double r315 = r274 * r314;
double r316 = r272 * r315;
double r317 = r316 * r285;
double r318 = r317 * r291;
double r319 = r266 / r318;
double r320 = r305 ? r319 : r293;
double r321 = r295 ? r303 : r320;
double r322 = r265 ? r293 : r321;
return r322;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -1.1607177078091258e+142 or 6.622566551191284e-131 < t Initial program 23.6
rmApplied add-cube-cbrt23.8
Applied unpow-prod-down23.8
Applied times-frac17.0
Applied associate-*l*15.7
rmApplied *-un-lft-identity15.7
Applied sqr-pow15.7
Applied times-frac9.3
Simplified9.3
rmApplied associate-*l*8.4
rmApplied add-cube-cbrt8.4
if -1.1607177078091258e+142 < t < -4.5548020353125647e-113Initial program 26.3
rmApplied add-cube-cbrt26.6
Applied unpow-prod-down26.6
Applied times-frac18.1
Applied associate-*l*13.9
rmApplied tan-quot13.9
Applied associate-*l/13.9
Applied frac-times13.9
Applied associate-*l/12.9
if -4.5548020353125647e-113 < t < 6.622566551191284e-131Initial program 64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac56.7
Applied associate-*l*56.7
rmApplied *-un-lft-identity56.7
Applied sqr-pow56.7
Applied times-frac46.0
Simplified46.0
rmApplied associate-*l*46.0
Taylor expanded around 0 40.8
Final simplification15.9
herbie shell --seed 2020025
(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))))