function code(NN, MM, N1, N)
return Float64(sqrt(Float64(Float64(Float64(4.0 * NN) * NN) / Float64(Float64(Float64(Float64(Float64(Float64(NN + MM) * Float64(NN - MM)) * Float64(Float64(N1 - N) + NN)) * Float64(Float64(N - N1) + NN)) * Float64(Float64(Float64(N + N1) - NN) + 1.0)) * Float64(Float64(Float64(N + N1) + NN) + 1.0)))) * Float64(Float64(Float64(2.0 * NN) + 1.0) * Float64(Float64(2.0 * NN) - 1.0)))
end
↓
function code(NN, MM, N1, N)
t_0 = Float64(Float64(NN + Float64(N - N1)) * Float64(N1 - Float64(-1.0 + Float64(NN - N))))
t_1 = Float64(Float64(Float64(N1 + N) - NN) + 1.0)
t_2 = Float64(sqrt(Float64(Float64(NN * Float64(4.0 * NN)) / Float64(Float64(Float64(Float64(Float64(N1 - N) - NN) * Float64(Float64(Float64(NN + MM) * Float64(NN - MM)) * Float64(Float64(N - N1) - NN))) * t_1) * Float64(1.0 + Float64(NN + Float64(N1 + N)))))) * Float64(Float64(1.0 + Float64(NN * 2.0)) * Float64(Float64(NN * 2.0) + -1.0)))
t_3 = Float64(NN / Float64(NN + Float64(Float64(N1 + N) + 1.0)))
t_4 = Float64(NN + Float64(N1 - N))
t_5 = Float64(fma(2.0, NN, 1.0) * fma(2.0, NN, -1.0))
tmp = 0.0
if (t_2 <= -2e-62)
tmp = Float64(sqrt(Float64(Float64(Float64(NN / Float64(t_4 / Float64(NN / Float64(Float64(NN + N) - N1)))) / Float64(NN + MM)) * Float64(Float64(4.0 / Float64(t_1 * Float64(NN - Float64(Float64(-1.0 - N) - N1)))) / Float64(NN - MM)))) * Float64(fma(NN, 2.0, 1.0) * fma(NN, 2.0, -1.0)));
elseif (t_2 <= 0.0)
tmp = Float64(Float64(sqrt(Float64(Float64(4.0 * NN) / Float64(Float64(NN + MM) * Float64(Float64(NN - MM) * Float64(t_4 * t_0))))) * sqrt(t_3)) * t_5);
else
tmp = Float64(t_5 * sqrt(Float64(Float64(t_3 * Float64(Float64(Float64(NN / Float64(NN - MM)) / t_4) * Float64(4.0 / Float64(NN + MM)))) / t_0)));
end
return tmp
end
code[NN_, MM_, N1_, N_] := N[(N[Sqrt[N[(N[(N[(4.0 * NN), $MachinePrecision] * NN), $MachinePrecision] / N[(N[(N[(N[(N[(N[(NN + MM), $MachinePrecision] * N[(NN - MM), $MachinePrecision]), $MachinePrecision] * N[(N[(N1 - N), $MachinePrecision] + NN), $MachinePrecision]), $MachinePrecision] * N[(N[(N - N1), $MachinePrecision] + NN), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N + N1), $MachinePrecision] - NN), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N + N1), $MachinePrecision] + NN), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(N[(2.0 * NN), $MachinePrecision] + 1.0), $MachinePrecision] * N[(N[(2.0 * NN), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[NN_, MM_, N1_, N_] := Block[{t$95$0 = N[(N[(NN + N[(N - N1), $MachinePrecision]), $MachinePrecision] * N[(N1 - N[(-1.0 + N[(NN - N), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N1 + N), $MachinePrecision] - NN), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[Sqrt[N[(N[(NN * N[(4.0 * NN), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N1 - N), $MachinePrecision] - NN), $MachinePrecision] * N[(N[(N[(NN + MM), $MachinePrecision] * N[(NN - MM), $MachinePrecision]), $MachinePrecision] * N[(N[(N - N1), $MachinePrecision] - NN), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] * N[(1.0 + N[(NN + N[(N1 + N), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(1.0 + N[(NN * 2.0), $MachinePrecision]), $MachinePrecision] * N[(N[(NN * 2.0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(NN / N[(NN + N[(N[(N1 + N), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(NN + N[(N1 - N), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(N[(2.0 * NN + 1.0), $MachinePrecision] * N[(2.0 * NN + -1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -2e-62], N[(N[Sqrt[N[(N[(N[(NN / N[(t$95$4 / N[(NN / N[(N[(NN + N), $MachinePrecision] - N1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(NN + MM), $MachinePrecision]), $MachinePrecision] * N[(N[(4.0 / N[(t$95$1 * N[(NN - N[(N[(-1.0 - N), $MachinePrecision] - N1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(NN - MM), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(NN * 2.0 + 1.0), $MachinePrecision] * N[(NN * 2.0 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 0.0], N[(N[(N[Sqrt[N[(N[(4.0 * NN), $MachinePrecision] / N[(N[(NN + MM), $MachinePrecision] * N[(N[(NN - MM), $MachinePrecision] * N[(t$95$4 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Sqrt[t$95$3], $MachinePrecision]), $MachinePrecision] * t$95$5), $MachinePrecision], N[(t$95$5 * N[Sqrt[N[(N[(t$95$3 * N[(N[(N[(NN / N[(NN - MM), $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision] * N[(4.0 / N[(NN + MM), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]]]]]]