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((k / t), 2.0);
double t_2 = 1.0 + (1.0 + t_1);
double t_3 = (((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * t_2;
double t_4 = cbrt((tan(k) * (2.0 + t_1))) * (t * cbrt(sin(k)));
double tmp;
if (t_3 <= -((double) INFINITY)) {
tmp = ((l * (l * 2.0)) / pow(t_4, 2.0)) / t_4;
} else if (t_3 <= 1e+163) {
tmp = 2.0 / (t_2 * (tan(k) * ((sin(k) * (pow(t, 3.0) / l)) / l)));
} else if (t_3 <= ((double) INFINITY)) {
tmp = l / pow(((t * cbrt(k)) / cbrt((l / k))), 3.0);
} else {
tmp = 2.0 / ((sin(k) * tan(k)) * ((k / l) * (t / (l / k))));
}
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((k / t), 2.0);
double t_2 = 1.0 + (1.0 + t_1);
double t_3 = (((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * t_2;
double t_4 = Math.cbrt((Math.tan(k) * (2.0 + t_1))) * (t * Math.cbrt(Math.sin(k)));
double tmp;
if (t_3 <= -Double.POSITIVE_INFINITY) {
tmp = ((l * (l * 2.0)) / Math.pow(t_4, 2.0)) / t_4;
} else if (t_3 <= 1e+163) {
tmp = 2.0 / (t_2 * (Math.tan(k) * ((Math.sin(k) * (Math.pow(t, 3.0) / l)) / l)));
} else if (t_3 <= Double.POSITIVE_INFINITY) {
tmp = l / Math.pow(((t * Math.cbrt(k)) / Math.cbrt((l / k))), 3.0);
} else {
tmp = 2.0 / ((Math.sin(k) * Math.tan(k)) * ((k / l) * (t / (l / k))));
}
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[(k / t), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(1.0 + N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = 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] * t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[(N[Power[N[(N[Tan[k], $MachinePrecision] * N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[(t * N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, (-Infinity)], N[(N[(N[(l * N[(l * 2.0), $MachinePrecision]), $MachinePrecision] / N[Power[t$95$4, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision], If[LessEqual[t$95$3, 1e+163], N[(2.0 / N[(t$95$2 * N[(N[Tan[k], $MachinePrecision] * N[(N[(N[Sin[k], $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, Infinity], N[(l / N[Power[N[(N[(t * N[Power[k, 1/3], $MachinePrecision]), $MachinePrecision] / N[Power[N[(l / k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(t / N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]