double code(double x) {
return ((((((1.0 + (0.1049934947 * (x * x))) + (0.0424060604 * ((x * x) * (x * x)))) + (0.0072644182 * (((x * x) * (x * x)) * (x * x)))) + (0.0005064034 * ((((x * x) * (x * x)) * (x * x)) * (x * x)))) + (0.0001789971 * (((((x * x) * (x * x)) * (x * x)) * (x * x)) * (x * x)))) / ((((((1.0 + (0.7715471019 * (x * x))) + (0.2909738639 * ((x * x) * (x * x)))) + (0.0694555761 * (((x * x) * (x * x)) * (x * x)))) + (0.0140005442 * ((((x * x) * (x * x)) * (x * x)) * (x * x)))) + (0.0008327945 * (((((x * x) * (x * x)) * (x * x)) * (x * x)) * (x * x)))) + ((2.0 * 0.0001789971) * ((((((x * x) * (x * x)) * (x * x)) * (x * x)) * (x * x)) * (x * x))))) * x;
}
↓
double code(double x) {
double t_0 = pow((x * x), 4.0);
double t_1 = pow((x * x), 5.0);
double tmp;
if (x <= -7000.0) {
tmp = ((1.0 / (x * x)) * (0.2514179000665374 / x)) + (0.5 / x);
} else if (x <= 50000000.0) {
tmp = x / ((1.0 + (fma(pow(x, 6.0), 0.0694555761, (t_0 * 0.0140005442)) + ((t_1 * (0.0008327945 + (x * (x * 0.0003579942)))) + fma((x * x), 0.7715471019, (0.2909738639 * pow(x, 4.0)))))) / (1.0 + (fma(0.0005064034, t_0, (t_1 * 0.0001789971)) + fma(x, (x * 0.1049934947), fma(0.0424060604, pow(x, 4.0), (pow(x, 6.0) * 0.0072644182))))));
} else {
tmp = 0.5 / x;
}
return tmp;
}
code[x_] := N[(N[(N[(N[(N[(N[(N[(1.0 + N[(0.1049934947 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0424060604 * N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0072644182 * N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0005064034 * N[(N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0001789971 * N[(N[(N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(1.0 + N[(0.7715471019 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.2909738639 * N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0694555761 * N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0140005442 * N[(N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.0008327945 * N[(N[(N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(2.0 * 0.0001789971), $MachinePrecision] * N[(N[(N[(N[(N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]
↓
code[x_] := Block[{t$95$0 = N[Power[N[(x * x), $MachinePrecision], 4.0], $MachinePrecision]}, Block[{t$95$1 = N[Power[N[(x * x), $MachinePrecision], 5.0], $MachinePrecision]}, If[LessEqual[x, -7000.0], N[(N[(N[(1.0 / N[(x * x), $MachinePrecision]), $MachinePrecision] * N[(0.2514179000665374 / x), $MachinePrecision]), $MachinePrecision] + N[(0.5 / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 50000000.0], N[(x / N[(N[(1.0 + N[(N[(N[Power[x, 6.0], $MachinePrecision] * 0.0694555761 + N[(t$95$0 * 0.0140005442), $MachinePrecision]), $MachinePrecision] + N[(N[(t$95$1 * N[(0.0008327945 + N[(x * N[(x * 0.0003579942), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(x * x), $MachinePrecision] * 0.7715471019 + N[(0.2909738639 * N[Power[x, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(N[(0.0005064034 * t$95$0 + N[(t$95$1 * 0.0001789971), $MachinePrecision]), $MachinePrecision] + N[(x * N[(x * 0.1049934947), $MachinePrecision] + N[(0.0424060604 * N[Power[x, 4.0], $MachinePrecision] + N[(N[Power[x, 6.0], $MachinePrecision] * 0.0072644182), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.5 / x), $MachinePrecision]]]]]