\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}\;\ell \le -4.023562739479679874738129852385552406136 \cdot 10^{144}:\\
\;\;\;\;\frac{\frac{2}{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\frac{\cos k}{{\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)}^{2}} \cdot \left(\frac{{\ell}^{2}}{{\left(\sqrt[3]{\sin k}\right)}^{2}} \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1}\right)\right) \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1}\right) \cdot 2\\
\end{array}double f(double t, double l, double k) {
double r105163 = 2.0;
double r105164 = t;
double r105165 = 3.0;
double r105166 = pow(r105164, r105165);
double r105167 = l;
double r105168 = r105167 * r105167;
double r105169 = r105166 / r105168;
double r105170 = k;
double r105171 = sin(r105170);
double r105172 = r105169 * r105171;
double r105173 = tan(r105170);
double r105174 = r105172 * r105173;
double r105175 = 1.0;
double r105176 = r105170 / r105164;
double r105177 = pow(r105176, r105163);
double r105178 = r105175 + r105177;
double r105179 = r105178 - r105175;
double r105180 = r105174 * r105179;
double r105181 = r105163 / r105180;
return r105181;
}
double f(double t, double l, double k) {
double r105182 = l;
double r105183 = -4.02356273947968e+144;
bool r105184 = r105182 <= r105183;
double r105185 = 2.0;
double r105186 = t;
double r105187 = cbrt(r105186);
double r105188 = r105187 * r105187;
double r105189 = 3.0;
double r105190 = pow(r105188, r105189);
double r105191 = r105190 / r105182;
double r105192 = pow(r105187, r105189);
double r105193 = r105192 / r105182;
double r105194 = k;
double r105195 = sin(r105194);
double r105196 = r105193 * r105195;
double r105197 = r105191 * r105196;
double r105198 = tan(r105194);
double r105199 = r105197 * r105198;
double r105200 = r105185 / r105199;
double r105201 = r105194 / r105186;
double r105202 = pow(r105201, r105185);
double r105203 = r105200 / r105202;
double r105204 = cos(r105194);
double r105205 = cbrt(r105195);
double r105206 = r105205 * r105205;
double r105207 = 2.0;
double r105208 = pow(r105206, r105207);
double r105209 = r105204 / r105208;
double r105210 = pow(r105182, r105207);
double r105211 = pow(r105205, r105207);
double r105212 = r105210 / r105211;
double r105213 = 1.0;
double r105214 = r105185 / r105207;
double r105215 = pow(r105194, r105214);
double r105216 = 1.0;
double r105217 = pow(r105186, r105216);
double r105218 = r105215 * r105217;
double r105219 = r105213 / r105218;
double r105220 = pow(r105219, r105216);
double r105221 = r105212 * r105220;
double r105222 = r105209 * r105221;
double r105223 = r105213 / r105215;
double r105224 = pow(r105223, r105216);
double r105225 = r105222 * r105224;
double r105226 = r105225 * r105185;
double r105227 = r105184 ? r105203 : r105226;
return r105227;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -4.02356273947968e+144Initial program 63.1
Simplified62.6
rmApplied add-cube-cbrt62.6
Applied unpow-prod-down62.6
Applied times-frac49.4
Applied associate-*l*49.4
if -4.02356273947968e+144 < l Initial program 46.9
Simplified38.5
Taylor expanded around inf 18.5
rmApplied sqr-pow18.5
Applied associate-*l*16.2
rmApplied add-cube-cbrt16.2
Applied times-frac16.1
Applied unpow-prod-down16.1
Applied associate-*l*14.3
Simplified14.3
rmApplied add-cube-cbrt14.7
Applied unpow-prod-down14.7
Applied times-frac14.4
Applied associate-*l*14.3
Final simplification17.1
herbie shell --seed 2019347
(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))))