\[\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\]
↓
\[\begin{array}{l}
\mathbf{if}\;\frac{angle}{180} \leq -0.0005:\\
\;\;\;\;\left(b \cdot b - a \cdot a\right) \cdot \sin \left(\sqrt[3]{angle} \cdot \left({\left(\sqrt[3]{angle}\right)}^{2} \cdot \left(\pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{elif}\;\frac{angle}{180} \leq 2 \cdot 10^{-23}:\\
\;\;\;\;\left(0.011111111111111112 \cdot \left(angle \cdot \left(b - a\right)\right)\right) \cdot \left(\pi \cdot \left(b + a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\sqrt[3]{\mathsf{fma}\left(b, b, -a \cdot a\right)}\right)}^{3} \cdot \sin \left(\mathsf{expm1}\left(\mathsf{log1p}\left(angle \cdot \left(\pi \cdot 0.011111111111111112\right)\right)\right)\right)\\
\end{array}
\]
double code(double a, double b, double angle) {
return ((2.0 * (pow(b, 2.0) - pow(a, 2.0))) * sin((((double) M_PI) * (angle / 180.0)))) * cos((((double) M_PI) * (angle / 180.0)));
}
↓
double code(double a, double b, double angle) {
double tmp;
if ((angle / 180.0) <= -0.0005) {
tmp = ((b * b) - (a * a)) * sin((cbrt(angle) * (pow(cbrt(angle), 2.0) * (((double) M_PI) * 0.011111111111111112))));
} else if ((angle / 180.0) <= 2e-23) {
tmp = (0.011111111111111112 * (angle * (b - a))) * (((double) M_PI) * (b + a));
} else {
tmp = pow(cbrt(fma(b, b, -(a * a))), 3.0) * sin(expm1(log1p((angle * (((double) M_PI) * 0.011111111111111112)))));
}
return tmp;
}
function code(a, b, angle)
return Float64(Float64(Float64(2.0 * Float64((b ^ 2.0) - (a ^ 2.0))) * sin(Float64(pi * Float64(angle / 180.0)))) * cos(Float64(pi * Float64(angle / 180.0))))
end
↓
function code(a, b, angle)
tmp = 0.0
if (Float64(angle / 180.0) <= -0.0005)
tmp = Float64(Float64(Float64(b * b) - Float64(a * a)) * sin(Float64(cbrt(angle) * Float64((cbrt(angle) ^ 2.0) * Float64(pi * 0.011111111111111112)))));
elseif (Float64(angle / 180.0) <= 2e-23)
tmp = Float64(Float64(0.011111111111111112 * Float64(angle * Float64(b - a))) * Float64(pi * Float64(b + a)));
else
tmp = Float64((cbrt(fma(b, b, Float64(-Float64(a * a)))) ^ 3.0) * sin(expm1(log1p(Float64(angle * Float64(pi * 0.011111111111111112))))));
end
return tmp
end
code[a_, b_, angle_] := N[(N[(N[(2.0 * N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
↓
code[a_, b_, angle_] := If[LessEqual[N[(angle / 180.0), $MachinePrecision], -0.0005], N[(N[(N[(b * b), $MachinePrecision] - N[(a * a), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(N[Power[angle, 1/3], $MachinePrecision] * N[(N[Power[N[Power[angle, 1/3], $MachinePrecision], 2.0], $MachinePrecision] * N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(angle / 180.0), $MachinePrecision], 2e-23], N[(N[(0.011111111111111112 * N[(angle * N[(b - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(Pi * N[(b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[N[Power[N[(b * b + (-N[(a * a), $MachinePrecision])), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision] * N[Sin[N[(Exp[N[Log[1 + N[(angle * N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
↓
\begin{array}{l}
\mathbf{if}\;\frac{angle}{180} \leq -0.0005:\\
\;\;\;\;\left(b \cdot b - a \cdot a\right) \cdot \sin \left(\sqrt[3]{angle} \cdot \left({\left(\sqrt[3]{angle}\right)}^{2} \cdot \left(\pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{elif}\;\frac{angle}{180} \leq 2 \cdot 10^{-23}:\\
\;\;\;\;\left(0.011111111111111112 \cdot \left(angle \cdot \left(b - a\right)\right)\right) \cdot \left(\pi \cdot \left(b + a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\sqrt[3]{\mathsf{fma}\left(b, b, -a \cdot a\right)}\right)}^{3} \cdot \sin \left(\mathsf{expm1}\left(\mathsf{log1p}\left(angle \cdot \left(\pi \cdot 0.011111111111111112\right)\right)\right)\right)\\
\end{array}