\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 -6.467661322206999577370568191648766431177 \cdot 10^{66}:\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \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)}{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \cos k}}\\
\mathbf{elif}\;k \le -1.923649962695752766452201338741225672413 \cdot 10^{-10}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{elif}\;k \le 4.398269809129575053279637547251081189272 \cdot 10^{-271}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{\sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{\sqrt{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}} \cdot \left(\frac{\sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}{\frac{\sqrt[3]{\ell}}{\sqrt{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}} \cdot \left(\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{else}:\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \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)}{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \cos k}}\\
\end{array}double f(double t, double l, double k) {
double r126057 = 2.0;
double r126058 = t;
double r126059 = 3.0;
double r126060 = pow(r126058, r126059);
double r126061 = l;
double r126062 = r126061 * r126061;
double r126063 = r126060 / r126062;
double r126064 = k;
double r126065 = sin(r126064);
double r126066 = r126063 * r126065;
double r126067 = tan(r126064);
double r126068 = r126066 * r126067;
double r126069 = 1.0;
double r126070 = r126064 / r126058;
double r126071 = pow(r126070, r126057);
double r126072 = r126069 + r126071;
double r126073 = r126072 + r126069;
double r126074 = r126068 * r126073;
double r126075 = r126057 / r126074;
return r126075;
}
double f(double t, double l, double k) {
double r126076 = k;
double r126077 = -6.4676613222069996e+66;
bool r126078 = r126076 <= r126077;
double r126079 = 2.0;
double r126080 = t;
double r126081 = cbrt(r126080);
double r126082 = r126081 * r126081;
double r126083 = 3.0;
double r126084 = 2.0;
double r126085 = r126083 / r126084;
double r126086 = pow(r126082, r126085);
double r126087 = pow(r126081, r126083);
double r126088 = l;
double r126089 = r126087 / r126088;
double r126090 = sin(r126076);
double r126091 = r126089 * r126090;
double r126092 = r126086 * r126091;
double r126093 = r126092 * r126090;
double r126094 = 1.0;
double r126095 = r126076 / r126080;
double r126096 = pow(r126095, r126079);
double r126097 = r126094 + r126096;
double r126098 = r126097 + r126094;
double r126099 = r126093 * r126098;
double r126100 = r126088 / r126086;
double r126101 = cos(r126076);
double r126102 = r126100 * r126101;
double r126103 = r126099 / r126102;
double r126104 = r126079 / r126103;
double r126105 = -1.9236499626957528e-10;
bool r126106 = r126076 <= r126105;
double r126107 = 3.0;
double r126108 = pow(r126080, r126107);
double r126109 = pow(r126090, r126084);
double r126110 = r126108 * r126109;
double r126111 = pow(r126088, r126084);
double r126112 = r126101 * r126111;
double r126113 = r126110 / r126112;
double r126114 = r126079 * r126113;
double r126115 = 1.0;
double r126116 = -1.0;
double r126117 = pow(r126116, r126083);
double r126118 = r126115 / r126117;
double r126119 = pow(r126118, r126094);
double r126120 = pow(r126076, r126084);
double r126121 = r126120 * r126109;
double r126122 = r126080 * r126121;
double r126123 = r126122 / r126112;
double r126124 = r126119 * r126123;
double r126125 = r126114 - r126124;
double r126126 = r126079 / r126125;
double r126127 = 4.398269809129575e-271;
bool r126128 = r126076 <= r126127;
double r126129 = cbrt(r126086);
double r126130 = r126129 * r126129;
double r126131 = cbrt(r126088);
double r126132 = r126131 * r126131;
double r126133 = sqrt(r126086);
double r126134 = r126132 / r126133;
double r126135 = r126130 / r126134;
double r126136 = r126131 / r126133;
double r126137 = r126129 / r126136;
double r126138 = r126137 * r126091;
double r126139 = r126135 * r126138;
double r126140 = tan(r126076);
double r126141 = r126139 * r126140;
double r126142 = r126141 * r126098;
double r126143 = r126079 / r126142;
double r126144 = r126128 ? r126143 : r126104;
double r126145 = r126106 ? r126126 : r126144;
double r126146 = r126078 ? r126104 : r126145;
return r126146;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -6.4676613222069996e+66 or 4.398269809129575e-271 < k Initial program 32.3
rmApplied add-cube-cbrt32.5
Applied unpow-prod-down32.5
Applied times-frac24.9
rmApplied sqr-pow24.9
Applied associate-/l*20.0
rmApplied associate-*l*18.7
rmApplied tan-quot18.7
Applied associate-*l/18.3
Applied frac-times17.4
Applied associate-*l/15.8
if -6.4676613222069996e+66 < k < -1.9236499626957528e-10Initial program 31.1
Taylor expanded around -inf 18.7
if -1.9236499626957528e-10 < k < 4.398269809129575e-271Initial program 35.5
rmApplied add-cube-cbrt35.7
Applied unpow-prod-down35.7
Applied times-frac28.5
rmApplied sqr-pow28.5
Applied associate-/l*24.2
rmApplied associate-*l*17.6
rmApplied add-sqr-sqrt17.6
Applied add-cube-cbrt17.6
Applied times-frac17.6
Applied add-cube-cbrt17.7
Applied times-frac17.7
Applied associate-*l*14.0
Final simplification15.7
herbie shell --seed 2019212
(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))))