\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.5084231617352398 \cdot 10^{-212} \lor \neg \left(k \le 3.5586471885166657 \cdot 10^{-146}\right):\\
\;\;\;\;\frac{2}{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell} \cdot \sin k\right)\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\
\end{array}double f(double t, double l, double k) {
double r129152 = 2.0;
double r129153 = t;
double r129154 = 3.0;
double r129155 = pow(r129153, r129154);
double r129156 = l;
double r129157 = r129156 * r129156;
double r129158 = r129155 / r129157;
double r129159 = k;
double r129160 = sin(r129159);
double r129161 = r129158 * r129160;
double r129162 = tan(r129159);
double r129163 = r129161 * r129162;
double r129164 = 1.0;
double r129165 = r129159 / r129153;
double r129166 = pow(r129165, r129152);
double r129167 = r129164 + r129166;
double r129168 = r129167 + r129164;
double r129169 = r129163 * r129168;
double r129170 = r129152 / r129169;
return r129170;
}
double f(double t, double l, double k) {
double r129171 = k;
double r129172 = -1.5084231617352398e-212;
bool r129173 = r129171 <= r129172;
double r129174 = 3.5586471885166657e-146;
bool r129175 = r129171 <= r129174;
double r129176 = !r129175;
bool r129177 = r129173 || r129176;
double r129178 = 2.0;
double r129179 = t;
double r129180 = cbrt(r129179);
double r129181 = 3.0;
double r129182 = pow(r129180, r129181);
double r129183 = l;
double r129184 = r129182 / r129183;
double r129185 = r129182 * r129184;
double r129186 = sin(r129171);
double r129187 = r129184 * r129186;
double r129188 = tan(r129171);
double r129189 = 1.0;
double r129190 = r129171 / r129179;
double r129191 = pow(r129190, r129178);
double r129192 = r129189 + r129191;
double r129193 = r129192 + r129189;
double r129194 = r129188 * r129193;
double r129195 = r129187 * r129194;
double r129196 = r129185 * r129195;
double r129197 = r129178 / r129196;
double r129198 = 0.3333333333333333;
double r129199 = r129198 * r129181;
double r129200 = pow(r129179, r129199);
double r129201 = r129200 / r129183;
double r129202 = r129201 * r129186;
double r129203 = r129185 * r129202;
double r129204 = r129203 * r129194;
double r129205 = r129178 / r129204;
double r129206 = r129177 ? r129197 : r129205;
return r129206;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -1.5084231617352398e-212 or 3.5586471885166657e-146 < k Initial program 31.2
rmApplied add-cube-cbrt31.3
Applied unpow-prod-down31.3
Applied times-frac23.8
Applied associate-*l*23.1
rmApplied *-un-lft-identity23.1
Applied unpow-prod-down23.1
Applied times-frac17.7
Simplified17.7
rmApplied associate-*l*17.5
rmApplied associate-*l*15.6
if -1.5084231617352398e-212 < k < 3.5586471885166657e-146Initial program 41.3
rmApplied add-cube-cbrt41.4
Applied unpow-prod-down41.4
Applied times-frac36.7
Applied associate-*l*26.3
rmApplied *-un-lft-identity26.3
Applied unpow-prod-down26.3
Applied times-frac20.6
Simplified20.6
rmApplied associate-*l*20.6
rmApplied pow1/342.6
Applied pow-pow20.5
Final simplification16.2
herbie shell --seed 2020021
(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))))