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(k), -2.0);
double t_2 = cbrt((l * 2.0)) / t;
double t_3 = 2.0 + pow((k / t), 2.0);
double tmp;
if (t <= -6.8e-46) {
tmp = pow((cbrt((l / (tan(k) * t_3))) * (cbrt((2.0 / (sin(k) / l))) / t)), 3.0);
} else if (t <= 0.00192) {
tmp = 2.0 * ((l * (cos(k) / k)) / (k * (t * (pow(sin(k), 2.0) / l))));
} else if (t <= 5e+151) {
tmp = (pow(t_2, 2.0) / (t_3 * sin(k))) * (t_2 / (tan(k) / l));
} else {
tmp = ((l / t) * (t_1 * (t_1 / t))) * (l / (t * pow(cbrt(k), 2.0)));
}
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(k), -2.0);
double t_2 = Math.cbrt((l * 2.0)) / t;
double t_3 = 2.0 + Math.pow((k / t), 2.0);
double tmp;
if (t <= -6.8e-46) {
tmp = Math.pow((Math.cbrt((l / (Math.tan(k) * t_3))) * (Math.cbrt((2.0 / (Math.sin(k) / l))) / t)), 3.0);
} else if (t <= 0.00192) {
tmp = 2.0 * ((l * (Math.cos(k) / k)) / (k * (t * (Math.pow(Math.sin(k), 2.0) / l))));
} else if (t <= 5e+151) {
tmp = (Math.pow(t_2, 2.0) / (t_3 * Math.sin(k))) * (t_2 / (Math.tan(k) / l));
} else {
tmp = ((l / t) * (t_1 * (t_1 / t))) * (l / (t * Math.pow(Math.cbrt(k), 2.0)));
}
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[Power[k, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, Block[{t$95$2 = N[(N[Power[N[(l * 2.0), $MachinePrecision], 1/3], $MachinePrecision] / t), $MachinePrecision]}, Block[{t$95$3 = N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -6.8e-46], N[Power[N[(N[Power[N[(l / N[(N[Tan[k], $MachinePrecision] * t$95$3), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[(N[Power[N[(2.0 / N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision], If[LessEqual[t, 0.00192], N[(2.0 * N[(N[(l * N[(N[Cos[k], $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / N[(k * N[(t * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5e+151], N[(N[(N[Power[t$95$2, 2.0], $MachinePrecision] / N[(t$95$3 * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$2 / N[(N[Tan[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / t), $MachinePrecision] * N[(t$95$1 * N[(t$95$1 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(t * N[Power[N[Power[k, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]