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