\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.2391434554590811 \cdot 10^{149}:\\
\;\;\;\;2 \cdot \left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\right)\\
\mathbf{elif}\;k \le -5.15729782667803904 \cdot 10^{-154}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\left|\sin k\right| \cdot \frac{\left|\sin k\right|}{\ell}}\\
\mathbf{elif}\;k \le 3.3090433695573168 \cdot 10^{-155}:\\
\;\;\;\;2 \cdot \left({\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\right)\\
\mathbf{elif}\;k \le 1.22181510276113542 \cdot 10^{154}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\left|\sin k\right| \cdot \frac{\left|\sin k\right|}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\right)\\
\end{array}double f(double t, double l, double k) {
double r87910 = 2.0;
double r87911 = t;
double r87912 = 3.0;
double r87913 = pow(r87911, r87912);
double r87914 = l;
double r87915 = r87914 * r87914;
double r87916 = r87913 / r87915;
double r87917 = k;
double r87918 = sin(r87917);
double r87919 = r87916 * r87918;
double r87920 = tan(r87917);
double r87921 = r87919 * r87920;
double r87922 = 1.0;
double r87923 = r87917 / r87911;
double r87924 = pow(r87923, r87910);
double r87925 = r87922 + r87924;
double r87926 = r87925 - r87922;
double r87927 = r87921 * r87926;
double r87928 = r87910 / r87927;
return r87928;
}
double f(double t, double l, double k) {
double r87929 = k;
double r87930 = -6.239143455459081e+149;
bool r87931 = r87929 <= r87930;
double r87932 = 2.0;
double r87933 = 1.0;
double r87934 = cbrt(r87933);
double r87935 = r87934 * r87934;
double r87936 = 2.0;
double r87937 = r87932 / r87936;
double r87938 = pow(r87929, r87937);
double r87939 = r87935 / r87938;
double r87940 = 1.0;
double r87941 = pow(r87939, r87940);
double r87942 = t;
double r87943 = pow(r87942, r87940);
double r87944 = r87938 * r87943;
double r87945 = r87934 / r87944;
double r87946 = pow(r87945, r87940);
double r87947 = cos(r87929);
double r87948 = l;
double r87949 = pow(r87948, r87936);
double r87950 = r87947 * r87949;
double r87951 = sin(r87929);
double r87952 = pow(r87951, r87936);
double r87953 = r87950 / r87952;
double r87954 = r87946 * r87953;
double r87955 = r87941 * r87954;
double r87956 = r87932 * r87955;
double r87957 = -5.157297826678039e-154;
bool r87958 = r87929 <= r87957;
double r87959 = pow(r87929, r87932);
double r87960 = r87959 * r87943;
double r87961 = r87933 / r87960;
double r87962 = pow(r87961, r87940);
double r87963 = r87962 * r87948;
double r87964 = r87963 * r87947;
double r87965 = fabs(r87951);
double r87966 = r87965 / r87948;
double r87967 = r87965 * r87966;
double r87968 = r87964 / r87967;
double r87969 = r87932 * r87968;
double r87970 = 3.3090433695573168e-155;
bool r87971 = r87929 <= r87970;
double r87972 = sqrt(r87933);
double r87973 = r87972 / r87938;
double r87974 = pow(r87973, r87940);
double r87975 = r87933 / r87944;
double r87976 = pow(r87975, r87940);
double r87977 = r87947 / r87965;
double r87978 = r87948 / r87966;
double r87979 = r87977 * r87978;
double r87980 = r87976 * r87979;
double r87981 = r87974 * r87980;
double r87982 = r87932 * r87981;
double r87983 = 1.2218151027611354e+154;
bool r87984 = r87929 <= r87983;
double r87985 = r87984 ? r87969 : r87956;
double r87986 = r87971 ? r87982 : r87985;
double r87987 = r87958 ? r87969 : r87986;
double r87988 = r87931 ? r87956 : r87987;
return r87988;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -6.239143455459081e+149 or 1.2218151027611354e+154 < k Initial program 38.4
Simplified33.8
Taylor expanded around inf 23.9
rmApplied sqr-pow23.9
Applied associate-*l*18.5
rmApplied add-cube-cbrt18.5
Applied times-frac18.2
Applied unpow-prod-down18.2
Applied associate-*l*14.9
if -6.239143455459081e+149 < k < -5.157297826678039e-154 or 3.3090433695573168e-155 < k < 1.2218151027611354e+154Initial program 53.4
Simplified43.8
Taylor expanded around inf 17.9
rmApplied sqr-pow17.9
Applied associate-*l*17.9
rmApplied add-sqr-sqrt17.9
Applied times-frac17.9
Simplified17.9
Simplified16.5
rmApplied frac-times14.7
Applied associate-*r/8.2
Simplified8.2
if -5.157297826678039e-154 < k < 3.3090433695573168e-155Initial program 64.0
Simplified64.0
Taylor expanded around inf 62.0
rmApplied sqr-pow62.0
Applied associate-*l*61.9
rmApplied add-sqr-sqrt61.9
Applied times-frac61.9
Simplified61.9
Simplified29.0
rmApplied add-sqr-sqrt29.0
Applied times-frac29.1
Applied unpow-prod-down29.1
Applied associate-*l*23.0
Simplified23.0
Final simplification11.4
herbie shell --seed 2020100
(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))))