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((((1.0 / t) * cbrt((l / (tan(k) * (2.0 + pow((k / t), 2.0)))))) * cbrt(((l * 2.0) / sin(k)))), 3.0);
double t_2 = pow(cbrt(k), 2.0);
double tmp;
if (t <= -6.5e-6) {
tmp = t_1;
} else if (t <= 6.9e-124) {
tmp = (2.0 * (((l / k) * (1.0 / k)) / t)) * ((l / sin(k)) / tan(k));
} else if (t <= 7.2e+229) {
tmp = t_1;
} else {
tmp = pow((t / (cbrt(l) / t_2)), -2.0) * (cbrt(l) * (l / (t * t_2)));
}
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((((1.0 / t) * Math.cbrt((l / (Math.tan(k) * (2.0 + Math.pow((k / t), 2.0)))))) * Math.cbrt(((l * 2.0) / Math.sin(k)))), 3.0);
double t_2 = Math.pow(Math.cbrt(k), 2.0);
double tmp;
if (t <= -6.5e-6) {
tmp = t_1;
} else if (t <= 6.9e-124) {
tmp = (2.0 * (((l / k) * (1.0 / k)) / t)) * ((l / Math.sin(k)) / Math.tan(k));
} else if (t <= 7.2e+229) {
tmp = t_1;
} else {
tmp = Math.pow((t / (Math.cbrt(l) / t_2)), -2.0) * (Math.cbrt(l) * (l / (t * t_2)));
}
return tmp;
}
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[(1.0 / t), $MachinePrecision] * N[Power[N[(l / N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[(l * 2.0), $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]}, Block[{t$95$2 = N[Power[N[Power[k, 1/3], $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[t, -6.5e-6], t$95$1, If[LessEqual[t, 6.9e-124], N[(N[(2.0 * N[(N[(N[(l / k), $MachinePrecision] * N[(1.0 / k), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] * N[(N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 7.2e+229], t$95$1, N[(N[Power[N[(t / N[(N[Power[l, 1/3], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision] * N[(N[Power[l, 1/3], $MachinePrecision] * N[(l / N[(t * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]