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((((cbrt(l) * cbrt((1.0 / sin(k)))) * cbrt((1.0 / tan(k)))) * (((cbrt(l) * cbrt(2.0)) / t) / cbrt((2.0 + pow((k / t), 2.0))))), 3.0);
double tmp;
if (t <= -2e-16) {
tmp = t_1;
} else if (t <= 1.0) {
tmp = (2.0 * ((l / k) / (t * k))) * ((l / sin(k)) / tan(k));
} else {
tmp = t_1;
}
return tmp;
}
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((((Math.cbrt(l) * Math.cbrt((1.0 / Math.sin(k)))) * Math.cbrt((1.0 / Math.tan(k)))) * (((Math.cbrt(l) * Math.cbrt(2.0)) / t) / Math.cbrt((2.0 + Math.pow((k / t), 2.0))))), 3.0);
double tmp;
if (t <= -2e-16) {
tmp = t_1;
} else if (t <= 1.0) {
tmp = (2.0 * ((l / k) / (t * k))) * ((l / Math.sin(k)) / Math.tan(k));
} else {
tmp = t_1;
}
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(Float64(Float64(cbrt(l) * cbrt(Float64(1.0 / sin(k)))) * cbrt(Float64(1.0 / tan(k)))) * Float64(Float64(Float64(cbrt(l) * cbrt(2.0)) / t) / cbrt(Float64(2.0 + (Float64(k / t) ^ 2.0))))) ^ 3.0
tmp = 0.0
if (t <= -2e-16)
tmp = t_1;
elseif (t <= 1.0)
tmp = Float64(Float64(2.0 * Float64(Float64(l / k) / Float64(t * k))) * Float64(Float64(l / sin(k)) / tan(k)));
else
tmp = t_1;
end
return 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[(N[(N[(N[Power[l, 1/3], $MachinePrecision] * N[Power[N[(1.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] * N[Power[N[(1.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[Power[l, 1/3], $MachinePrecision] * N[Power[2.0, 1/3], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] / N[Power[N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]}, If[LessEqual[t, -2e-16], t$95$1, If[LessEqual[t, 1.0], N[(N[(2.0 * N[(N[(l / k), $MachinePrecision] / N[(t * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
