\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}\;k \le -1.932296528705668984693986702426810627746 \cdot 10^{146}:\\
\;\;\;\;\frac{2}{\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left({\left(\sqrt[3]{\sqrt[3]{t}} \cdot \sqrt[3]{\sqrt[3]{t}}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{\sqrt[3]{t}}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right)\right) \cdot \tan k\right)}\\
\mathbf{elif}\;k \le -2.465489873106864392898978142261571099549 \cdot 10^{69}:\\
\;\;\;\;\frac{2}{2 \cdot \left({\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({t}^{3} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}\right) + {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({\left(\sin k\right)}^{2} \cdot \left({k}^{2} \cdot t\right)\right)}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{elif}\;k \le -3.515615117987834914814812597504544001232 \cdot 10^{-226}:\\
\;\;\;\;\frac{2}{\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right)\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le 496413.6209539859555661678314208984375:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left(\frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\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}\;k \le 2.402994973162251042133283519939024529294 \cdot 10^{139}:\\
\;\;\;\;\frac{2}{2 \cdot \left({\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({t}^{3} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}\right) + {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({\left(\sin k\right)}^{2} \cdot \left({k}^{2} \cdot t\right)\right)}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{elif}\;k \le 3.506733478446760175749683528136847068451 \cdot 10^{202}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}}\right) \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right) \cdot \tan k\right)\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le 7.387345331331856574424040723290335999974 \cdot 10^{224}:\\
\;\;\;\;\frac{2}{2 \cdot \left({\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({t}^{3} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}\right) + {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({\left(\sin k\right)}^{2} \cdot \left({k}^{2} \cdot t\right)\right)}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k\right) \cdot \left(\sqrt[3]{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1} \cdot \sqrt[3]{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}\right)\right) \cdot \sqrt[3]{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}\\
\end{array}double f(double t, double l, double k) {
double r124160 = 2.0;
double r124161 = t;
double r124162 = 3.0;
double r124163 = pow(r124161, r124162);
double r124164 = l;
double r124165 = r124164 * r124164;
double r124166 = r124163 / r124165;
double r124167 = k;
double r124168 = sin(r124167);
double r124169 = r124166 * r124168;
double r124170 = tan(r124167);
double r124171 = r124169 * r124170;
double r124172 = 1.0;
double r124173 = r124167 / r124161;
double r124174 = pow(r124173, r124160);
double r124175 = r124172 + r124174;
double r124176 = r124175 + r124172;
double r124177 = r124171 * r124176;
double r124178 = r124160 / r124177;
return r124178;
}
double f(double t, double l, double k) {
double r124179 = k;
double r124180 = -1.932296528705669e+146;
bool r124181 = r124179 <= r124180;
double r124182 = 2.0;
double r124183 = 1.0;
double r124184 = t;
double r124185 = r124179 / r124184;
double r124186 = pow(r124185, r124182);
double r124187 = r124183 + r124186;
double r124188 = r124187 + r124183;
double r124189 = cbrt(r124184);
double r124190 = 3.0;
double r124191 = pow(r124189, r124190);
double r124192 = l;
double r124193 = cbrt(r124192);
double r124194 = r124193 * r124193;
double r124195 = r124191 / r124194;
double r124196 = r124191 / r124193;
double r124197 = cbrt(r124189);
double r124198 = r124197 * r124197;
double r124199 = pow(r124198, r124190);
double r124200 = pow(r124197, r124190);
double r124201 = r124200 / r124192;
double r124202 = sin(r124179);
double r124203 = r124201 * r124202;
double r124204 = r124199 * r124203;
double r124205 = r124196 * r124204;
double r124206 = r124195 * r124205;
double r124207 = tan(r124179);
double r124208 = r124206 * r124207;
double r124209 = r124188 * r124208;
double r124210 = r124182 / r124209;
double r124211 = -2.4654898731068644e+69;
bool r124212 = r124179 <= r124211;
double r124213 = 1.0;
double r124214 = -1.0;
double r124215 = pow(r124214, r124190);
double r124216 = r124213 / r124215;
double r124217 = pow(r124216, r124183);
double r124218 = cbrt(r124214);
double r124219 = 9.0;
double r124220 = pow(r124218, r124219);
double r124221 = 3.0;
double r124222 = pow(r124184, r124221);
double r124223 = 2.0;
double r124224 = pow(r124202, r124223);
double r124225 = r124222 * r124224;
double r124226 = r124220 * r124225;
double r124227 = cos(r124179);
double r124228 = pow(r124192, r124223);
double r124229 = r124227 * r124228;
double r124230 = r124226 / r124229;
double r124231 = r124217 * r124230;
double r124232 = r124182 * r124231;
double r124233 = pow(r124179, r124223);
double r124234 = r124233 * r124184;
double r124235 = r124224 * r124234;
double r124236 = r124220 * r124235;
double r124237 = r124236 / r124229;
double r124238 = r124217 * r124237;
double r124239 = r124232 + r124238;
double r124240 = r124182 / r124239;
double r124241 = -3.515615117987835e-226;
bool r124242 = r124179 <= r124241;
double r124243 = r124191 / r124192;
double r124244 = r124243 * r124202;
double r124245 = r124196 * r124244;
double r124246 = r124245 * r124207;
double r124247 = r124195 * r124246;
double r124248 = r124247 * r124188;
double r124249 = r124182 / r124248;
double r124250 = 496413.62095398596;
bool r124251 = r124179 <= r124250;
double r124252 = 0.3333333333333333;
double r124253 = r124252 * r124190;
double r124254 = pow(r124184, r124253);
double r124255 = r124254 / r124192;
double r124256 = r124255 * r124202;
double r124257 = r124196 * r124256;
double r124258 = r124195 * r124257;
double r124259 = r124258 * r124207;
double r124260 = r124259 * r124188;
double r124261 = r124182 / r124260;
double r124262 = 2.402994973162251e+139;
bool r124263 = r124179 <= r124262;
double r124264 = 3.50673347844676e+202;
bool r124265 = r124179 <= r124264;
double r124266 = r124195 * r124196;
double r124267 = r124244 * r124207;
double r124268 = r124266 * r124267;
double r124269 = r124268 * r124188;
double r124270 = r124182 / r124269;
double r124271 = 7.3873453313318566e+224;
bool r124272 = r124179 <= r124271;
double r124273 = r124195 * r124245;
double r124274 = r124273 * r124207;
double r124275 = cbrt(r124188);
double r124276 = r124275 * r124275;
double r124277 = r124274 * r124276;
double r124278 = r124277 * r124275;
double r124279 = r124182 / r124278;
double r124280 = r124272 ? r124240 : r124279;
double r124281 = r124265 ? r124270 : r124280;
double r124282 = r124263 ? r124240 : r124281;
double r124283 = r124251 ? r124261 : r124282;
double r124284 = r124242 ? r124249 : r124283;
double r124285 = r124212 ? r124240 : r124284;
double r124286 = r124181 ? r124210 : r124285;
return r124286;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -1.932296528705669e+146Initial program 34.9
rmApplied add-cube-cbrt34.9
Applied unpow-prod-down34.9
Applied times-frac27.8
Applied associate-*l*27.8
rmApplied add-cube-cbrt27.8
Applied unpow-prod-down27.8
Applied times-frac22.4
rmApplied associate-*l*22.4
rmApplied *-un-lft-identity22.4
Applied add-cube-cbrt22.4
Applied unpow-prod-down22.4
Applied times-frac22.4
Applied associate-*l*22.4
if -1.932296528705669e+146 < k < -2.4654898731068644e+69 or 496413.62095398596 < k < 2.402994973162251e+139 or 3.50673347844676e+202 < k < 7.3873453313318566e+224Initial program 32.7
rmApplied add-cube-cbrt32.8
Applied unpow-prod-down32.8
Applied times-frac25.2
Applied associate-*l*25.2
Taylor expanded around -inf 20.4
if -2.4654898731068644e+69 < k < -3.515615117987835e-226Initial program 30.4
rmApplied add-cube-cbrt30.6
Applied unpow-prod-down30.6
Applied times-frac23.0
Applied associate-*l*20.3
rmApplied add-cube-cbrt20.3
Applied unpow-prod-down20.3
Applied times-frac13.4
rmApplied associate-*l*12.8
rmApplied associate-*l*9.9
if -3.515615117987835e-226 < k < 496413.62095398596Initial program 34.2
rmApplied add-cube-cbrt34.4
Applied unpow-prod-down34.4
Applied times-frac27.6
Applied associate-*l*21.4
rmApplied add-cube-cbrt21.4
Applied unpow-prod-down21.4
Applied times-frac16.3
rmApplied associate-*l*13.0
rmApplied pow1/339.1
Applied pow-pow12.8
if 2.402994973162251e+139 < k < 3.50673347844676e+202Initial program 30.6
rmApplied add-cube-cbrt30.7
Applied unpow-prod-down30.7
Applied times-frac24.6
Applied associate-*l*24.6
rmApplied add-cube-cbrt24.6
Applied unpow-prod-down24.6
Applied times-frac19.4
rmApplied associate-*l*19.4
if 7.3873453313318566e+224 < k Initial program 32.0
rmApplied add-cube-cbrt32.0
Applied unpow-prod-down32.0
Applied times-frac26.2
Applied associate-*l*26.2
rmApplied add-cube-cbrt26.2
Applied unpow-prod-down26.2
Applied times-frac21.3
rmApplied associate-*l*21.3
rmApplied add-cube-cbrt21.3
Applied associate-*r*21.3
Final simplification16.9
herbie shell --seed 2019303
(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))))