\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{{\left(2 \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}} \cdot e^{\log \pi \cdot \frac{1 - k}{2}}double f(double k, double n) {
double r2285961 = 1.0;
double r2285962 = k;
double r2285963 = sqrt(r2285962);
double r2285964 = r2285961 / r2285963;
double r2285965 = 2.0;
double r2285966 = atan2(1.0, 0.0);
double r2285967 = r2285965 * r2285966;
double r2285968 = n;
double r2285969 = r2285967 * r2285968;
double r2285970 = r2285961 - r2285962;
double r2285971 = r2285970 / r2285965;
double r2285972 = pow(r2285969, r2285971);
double r2285973 = r2285964 * r2285972;
return r2285973;
}
double f(double k, double n) {
double r2285974 = 2.0;
double r2285975 = n;
double r2285976 = r2285974 * r2285975;
double r2285977 = 1.0;
double r2285978 = k;
double r2285979 = r2285977 - r2285978;
double r2285980 = r2285979 / r2285974;
double r2285981 = pow(r2285976, r2285980);
double r2285982 = sqrt(r2285978);
double r2285983 = r2285981 / r2285982;
double r2285984 = atan2(1.0, 0.0);
double r2285985 = log(r2285984);
double r2285986 = r2285985 * r2285980;
double r2285987 = exp(r2285986);
double r2285988 = r2285983 * r2285987;
return r2285988;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
Simplified0.3
rmApplied *-un-lft-identity0.3
Applied unpow-prod-down0.4
Applied times-frac0.4
Simplified0.4
rmApplied add-exp-log0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019128
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))