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