| Alternative 1 | |
|---|---|
| Accuracy | 78.9% |
| Cost | 53828 |

(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (/ k t) 2.0)))
(if (<=
(/
2.0
(*
(* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k))
(+ 1.0 (+ 1.0 t_1))))
5e+130)
(* (* l l) (/ 2.0 (* (tan k) (* (+ 2.0 t_1) (* (pow t 3.0) (sin k))))))
(* 2.0 (* (/ (cos k) t) (/ (pow (/ l k) 2.0) (pow (sin k) 2.0)))))))double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) + 1.0));
}
double code(double t, double l, double k) {
double t_1 = pow((k / t), 2.0);
double tmp;
if ((2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * (1.0 + (1.0 + t_1)))) <= 5e+130) {
tmp = (l * l) * (2.0 / (tan(k) * ((2.0 + t_1) * (pow(t, 3.0) * sin(k)))));
} else {
tmp = 2.0 * ((cos(k) / t) * (pow((l / k), 2.0) / pow(sin(k), 2.0)));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) + 1.0d0))
end function
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = (k / t) ** 2.0d0
if ((2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * (1.0d0 + (1.0d0 + t_1)))) <= 5d+130) then
tmp = (l * l) * (2.0d0 / (tan(k) * ((2.0d0 + t_1) * ((t ** 3.0d0) * sin(k)))))
else
tmp = 2.0d0 * ((cos(k) / t) * (((l / k) ** 2.0d0) / (sin(k) ** 2.0d0)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) + 1.0));
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow((k / t), 2.0);
double tmp;
if ((2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * (1.0 + (1.0 + t_1)))) <= 5e+130) {
tmp = (l * l) * (2.0 / (Math.tan(k) * ((2.0 + t_1) * (Math.pow(t, 3.0) * Math.sin(k)))));
} else {
tmp = 2.0 * ((Math.cos(k) / t) * (Math.pow((l / k), 2.0) / Math.pow(Math.sin(k), 2.0)));
}
return tmp;
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) + 1.0))
def code(t, l, k): t_1 = math.pow((k / t), 2.0) tmp = 0 if (2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * (1.0 + (1.0 + t_1)))) <= 5e+130: tmp = (l * l) * (2.0 / (math.tan(k) * ((2.0 + t_1) * (math.pow(t, 3.0) * math.sin(k))))) else: tmp = 2.0 * ((math.cos(k) / t) * (math.pow((l / k), 2.0) / math.pow(math.sin(k), 2.0))) return tmp
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) + 1.0))) end
function code(t, l, k) t_1 = Float64(k / t) ^ 2.0 tmp = 0.0 if (Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(1.0 + Float64(1.0 + t_1)))) <= 5e+130) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(tan(k) * Float64(Float64(2.0 + t_1) * Float64((t ^ 3.0) * sin(k)))))); else tmp = Float64(2.0 * Float64(Float64(cos(k) / t) * Float64((Float64(l / k) ^ 2.0) / (sin(k) ^ 2.0)))); end return tmp end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) + 1.0)); end
function tmp_2 = code(t, l, k) t_1 = (k / t) ^ 2.0; tmp = 0.0; if ((2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * (1.0 + (1.0 + t_1)))) <= 5e+130) tmp = (l * l) * (2.0 / (tan(k) * ((2.0 + t_1) * ((t ^ 3.0) * sin(k))))); else tmp = 2.0 * ((cos(k) / t) * (((l / k) ^ 2.0) / (sin(k) ^ 2.0))); end tmp_2 = tmp; end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 5e+130], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Tan[k], $MachinePrecision] * N[(N[(2.0 + t$95$1), $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] / t), $MachinePrecision] * N[(N[Power[N[(l / k), $MachinePrecision], 2.0], $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\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}
t_1 := {\left(\frac{k}{t}\right)}^{2}\\
\mathbf{if}\;\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(1 + \left(1 + t_1\right)\right)} \leq 5 \cdot 10^{+130}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{\tan k \cdot \left(\left(2 + t_1\right) \cdot \left({t}^{3} \cdot \sin k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k}{t} \cdot \frac{{\left(\frac{\ell}{k}\right)}^{2}}{{\sin k}^{2}}\right)\\
\end{array}
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
Results
if (/.f64 2 (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t 3) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 1 (pow.f64 (/.f64 k t) 2)) 1))) < 4.9999999999999996e130Initial program 85.2%
Simplified86.4%
[Start]85.2% | \[ \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)}
\] |
|---|---|
associate-/l/ [<=]85.1% | \[ \color{blue}{\frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k}}
\] |
associate-*l/ [=>]86.4% | \[ \frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\color{blue}{\frac{{t}^{3} \cdot \sin k}{\ell \cdot \ell}} \cdot \tan k}
\] |
associate-*l/ [=>]86.4% | \[ \frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\color{blue}{\frac{\left({t}^{3} \cdot \sin k\right) \cdot \tan k}{\ell \cdot \ell}}}
\] |
associate-/r/ [=>]86.4% | \[ \color{blue}{\frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\left({t}^{3} \cdot \sin k\right) \cdot \tan k} \cdot \left(\ell \cdot \ell\right)}
\] |
*-commutative [=>]86.4% | \[ \color{blue}{\left(\ell \cdot \ell\right) \cdot \frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\left({t}^{3} \cdot \sin k\right) \cdot \tan k}}
\] |
associate-/l/ [=>]86.4% | \[ \left(\ell \cdot \ell\right) \cdot \color{blue}{\frac{2}{\left(\left({t}^{3} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}}
\] |
associate-*r* [<=]86.4% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\color{blue}{\left({t}^{3} \cdot \sin k\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}}
\] |
*-commutative [=>]86.4% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\left({t}^{3} \cdot \sin k\right) \cdot \color{blue}{\left(\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \tan k\right)}}
\] |
associate-*r* [=>]86.4% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\color{blue}{\left(\left({t}^{3} \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right) \cdot \tan k}}
\] |
*-commutative [=>]86.4% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\color{blue}{\tan k \cdot \left(\left({t}^{3} \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}}
\] |
if 4.9999999999999996e130 < (/.f64 2 (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t 3) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 1 (pow.f64 (/.f64 k t) 2)) 1))) Initial program 21.2%
Simplified21.3%
[Start]21.2% | \[ \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)}
\] |
|---|---|
associate-/l/ [<=]21.2% | \[ \color{blue}{\frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k}}
\] |
associate-*l/ [=>]21.2% | \[ \frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\color{blue}{\frac{{t}^{3} \cdot \sin k}{\ell \cdot \ell}} \cdot \tan k}
\] |
associate-*l/ [=>]21.2% | \[ \frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\color{blue}{\frac{\left({t}^{3} \cdot \sin k\right) \cdot \tan k}{\ell \cdot \ell}}}
\] |
associate-/r/ [=>]21.3% | \[ \color{blue}{\frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\left({t}^{3} \cdot \sin k\right) \cdot \tan k} \cdot \left(\ell \cdot \ell\right)}
\] |
*-commutative [=>]21.3% | \[ \color{blue}{\left(\ell \cdot \ell\right) \cdot \frac{\frac{2}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\left({t}^{3} \cdot \sin k\right) \cdot \tan k}}
\] |
associate-/l/ [=>]21.3% | \[ \left(\ell \cdot \ell\right) \cdot \color{blue}{\frac{2}{\left(\left({t}^{3} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}}
\] |
associate-*r* [<=]21.3% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\color{blue}{\left({t}^{3} \cdot \sin k\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}}
\] |
*-commutative [=>]21.3% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\left({t}^{3} \cdot \sin k\right) \cdot \color{blue}{\left(\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \tan k\right)}}
\] |
associate-*r* [=>]21.3% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\color{blue}{\left(\left({t}^{3} \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right) \cdot \tan k}}
\] |
*-commutative [=>]21.3% | \[ \left(\ell \cdot \ell\right) \cdot \frac{2}{\color{blue}{\tan k \cdot \left(\left({t}^{3} \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}}
\] |
Taylor expanded in k around inf 59.7%
Simplified80.0%
[Start]59.7% | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
*-commutative [=>]59.7% | \[ 2 \cdot \frac{\color{blue}{{\ell}^{2} \cdot \cos k}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
times-frac [=>]60.1% | \[ 2 \cdot \color{blue}{\left(\frac{{\ell}^{2}}{{k}^{2}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)}
\] |
unpow2 [=>]60.1% | \[ 2 \cdot \left(\frac{\color{blue}{\ell \cdot \ell}}{{k}^{2}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
unpow2 [=>]60.1% | \[ 2 \cdot \left(\frac{\ell \cdot \ell}{\color{blue}{k \cdot k}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
times-frac [=>]80.0% | \[ 2 \cdot \left(\color{blue}{\left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right)} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
*-commutative [=>]80.0% | \[ 2 \cdot \left(\left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right) \cdot \frac{\cos k}{\color{blue}{t \cdot {\sin k}^{2}}}\right)
\] |
Taylor expanded in l around 0 59.7%
Simplified81.6%
[Start]59.7% | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
*-commutative [=>]59.7% | \[ 2 \cdot \frac{\color{blue}{{\ell}^{2} \cdot \cos k}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
*-commutative [<=]59.7% | \[ 2 \cdot \frac{{\ell}^{2} \cdot \cos k}{{k}^{2} \cdot \color{blue}{\left(t \cdot {\sin k}^{2}\right)}}
\] |
times-frac [=>]60.1% | \[ 2 \cdot \color{blue}{\left(\frac{{\ell}^{2}}{{k}^{2}} \cdot \frac{\cos k}{t \cdot {\sin k}^{2}}\right)}
\] |
unpow2 [=>]60.1% | \[ 2 \cdot \left(\frac{\color{blue}{\ell \cdot \ell}}{{k}^{2}} \cdot \frac{\cos k}{t \cdot {\sin k}^{2}}\right)
\] |
unpow2 [=>]60.1% | \[ 2 \cdot \left(\frac{\ell \cdot \ell}{\color{blue}{k \cdot k}} \cdot \frac{\cos k}{t \cdot {\sin k}^{2}}\right)
\] |
times-frac [=>]80.0% | \[ 2 \cdot \left(\color{blue}{\left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right)} \cdot \frac{\cos k}{t \cdot {\sin k}^{2}}\right)
\] |
unpow2 [<=]80.0% | \[ 2 \cdot \left(\color{blue}{{\left(\frac{\ell}{k}\right)}^{2}} \cdot \frac{\cos k}{t \cdot {\sin k}^{2}}\right)
\] |
associate-*r/ [=>]80.0% | \[ 2 \cdot \color{blue}{\frac{{\left(\frac{\ell}{k}\right)}^{2} \cdot \cos k}{t \cdot {\sin k}^{2}}}
\] |
times-frac [=>]81.6% | \[ 2 \cdot \color{blue}{\left(\frac{{\left(\frac{\ell}{k}\right)}^{2}}{t} \cdot \frac{\cos k}{{\sin k}^{2}}\right)}
\] |
Applied egg-rr81.6%
[Start]81.6% | \[ 2 \cdot \left(\frac{{\left(\frac{\ell}{k}\right)}^{2}}{t} \cdot \frac{\cos k}{{\sin k}^{2}}\right)
\] |
|---|---|
unpow2 [=>]81.6% | \[ 2 \cdot \left(\frac{\color{blue}{\frac{\ell}{k} \cdot \frac{\ell}{k}}}{t} \cdot \frac{\cos k}{{\sin k}^{2}}\right)
\] |
Taylor expanded in l around 0 59.7%
Simplified81.7%
[Start]59.7% | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
associate-/r* [=>]60.2% | \[ 2 \cdot \color{blue}{\frac{\frac{\cos k \cdot {\ell}^{2}}{{k}^{2}}}{{\sin k}^{2} \cdot t}}
\] |
associate-*r/ [<=]60.2% | \[ 2 \cdot \frac{\color{blue}{\cos k \cdot \frac{{\ell}^{2}}{{k}^{2}}}}{{\sin k}^{2} \cdot t}
\] |
unpow2 [=>]60.2% | \[ 2 \cdot \frac{\cos k \cdot \frac{\color{blue}{\ell \cdot \ell}}{{k}^{2}}}{{\sin k}^{2} \cdot t}
\] |
unpow2 [=>]60.2% | \[ 2 \cdot \frac{\cos k \cdot \frac{\ell \cdot \ell}{\color{blue}{k \cdot k}}}{{\sin k}^{2} \cdot t}
\] |
times-frac [=>]80.0% | \[ 2 \cdot \frac{\cos k \cdot \color{blue}{\left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right)}}{{\sin k}^{2} \cdot t}
\] |
unpow2 [<=]80.0% | \[ 2 \cdot \frac{\cos k \cdot \color{blue}{{\left(\frac{\ell}{k}\right)}^{2}}}{{\sin k}^{2} \cdot t}
\] |
*-commutative [<=]80.0% | \[ 2 \cdot \frac{\cos k \cdot {\left(\frac{\ell}{k}\right)}^{2}}{\color{blue}{t \cdot {\sin k}^{2}}}
\] |
times-frac [=>]81.7% | \[ 2 \cdot \color{blue}{\left(\frac{\cos k}{t} \cdot \frac{{\left(\frac{\ell}{k}\right)}^{2}}{{\sin k}^{2}}\right)}
\] |
Final simplification84.4%
| Alternative 1 | |
|---|---|
| Accuracy | 78.9% |
| Cost | 53828 |
| Alternative 2 | |
|---|---|
| Accuracy | 76.7% |
| Cost | 20489 |
| Alternative 3 | |
|---|---|
| Accuracy | 77.1% |
| Cost | 20489 |
| Alternative 4 | |
|---|---|
| Accuracy | 79.4% |
| Cost | 14540 |
| Alternative 5 | |
|---|---|
| Accuracy | 79.4% |
| Cost | 14540 |
| Alternative 6 | |
|---|---|
| Accuracy | 71.9% |
| Cost | 14025 |
| Alternative 7 | |
|---|---|
| Accuracy | 66.0% |
| Cost | 7305 |
| Alternative 8 | |
|---|---|
| Accuracy | 63.2% |
| Cost | 1481 |
| Alternative 9 | |
|---|---|
| Accuracy | 60.5% |
| Cost | 1353 |
| Alternative 10 | |
|---|---|
| Accuracy | 60.5% |
| Cost | 1353 |
| Alternative 11 | |
|---|---|
| Accuracy | 61.9% |
| Cost | 1353 |
| Alternative 12 | |
|---|---|
| Accuracy | 54.8% |
| Cost | 832 |
herbie shell --seed 2023165
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
:precision binary64
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))