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 = cbrt(l) * cbrt((1.0 / k));
double t_2 = pow(sin(k), 2.0);
double t_3 = cbrt((sin(k) * (tan(k) * (2.0 + pow((k / t), 2.0)))));
double t_4 = t * t_3;
double t_5 = pow(cbrt(l), 2.0);
double t_6 = t_3 * (t / t_5);
double tmp;
if (k <= -5.8e+48) {
tmp = 2.0 / (t_2 * ((k / l) * ((k / (l / t)) / cos(k))));
} else if (k <= -2.2e-155) {
tmp = (cbrt(l) * ((2.0 / t_4) / (1.0 / cbrt(l)))) / pow((1.0 / (t_5 / t_4)), 2.0);
} else if (k <= 4.3e-155) {
tmp = pow(((t_1 * t_1) / t), 3.0);
} else if (k <= 1.15e+43) {
tmp = (2.0 / t_6) / pow(t_6, 2.0);
} else {
tmp = 2.0 * (((l / k) * (l / k)) * (cos(k) / (t_2 * t)));
}
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.cbrt(l) * Math.cbrt((1.0 / k));
double t_2 = Math.pow(Math.sin(k), 2.0);
double t_3 = Math.cbrt((Math.sin(k) * (Math.tan(k) * (2.0 + Math.pow((k / t), 2.0)))));
double t_4 = t * t_3;
double t_5 = Math.pow(Math.cbrt(l), 2.0);
double t_6 = t_3 * (t / t_5);
double tmp;
if (k <= -5.8e+48) {
tmp = 2.0 / (t_2 * ((k / l) * ((k / (l / t)) / Math.cos(k))));
} else if (k <= -2.2e-155) {
tmp = (Math.cbrt(l) * ((2.0 / t_4) / (1.0 / Math.cbrt(l)))) / Math.pow((1.0 / (t_5 / t_4)), 2.0);
} else if (k <= 4.3e-155) {
tmp = Math.pow(((t_1 * t_1) / t), 3.0);
} else if (k <= 1.15e+43) {
tmp = (2.0 / t_6) / Math.pow(t_6, 2.0);
} else {
tmp = 2.0 * (((l / k) * (l / k)) * (Math.cos(k) / (t_2 * t)));
}
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[(N[Power[l, 1/3], $MachinePrecision] * N[Power[N[(1.0 / k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[Power[N[(N[Sin[k], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$4 = N[(t * t$95$3), $MachinePrecision]}, Block[{t$95$5 = N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$6 = N[(t$95$3 * N[(t / t$95$5), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, -5.8e+48], N[(2.0 / N[(t$95$2 * N[(N[(k / l), $MachinePrecision] * N[(N[(k / N[(l / t), $MachinePrecision]), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, -2.2e-155], N[(N[(N[Power[l, 1/3], $MachinePrecision] * N[(N[(2.0 / t$95$4), $MachinePrecision] / N[(1.0 / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[(1.0 / N[(t$95$5 / t$95$4), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 4.3e-155], N[Power[N[(N[(t$95$1 * t$95$1), $MachinePrecision] / t), $MachinePrecision], 3.0], $MachinePrecision], If[LessEqual[k, 1.15e+43], N[(N[(2.0 / t$95$6), $MachinePrecision] / N[Power[t$95$6, 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[(l / k), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k], $MachinePrecision] / N[(t$95$2 * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]