\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\left(\frac{1}{\sqrt{k}} \cdot \left({2}^{\left(\frac{1 - k}{2}\right)} \cdot {\pi}^{\left(\frac{1 - k}{2}\right)}\right)\right) \cdot {n}^{\left(\frac{1 - k}{2}\right)}double f(double k, double n) {
double r96915 = 1.0;
double r96916 = k;
double r96917 = sqrt(r96916);
double r96918 = r96915 / r96917;
double r96919 = 2.0;
double r96920 = atan2(1.0, 0.0);
double r96921 = r96919 * r96920;
double r96922 = n;
double r96923 = r96921 * r96922;
double r96924 = r96915 - r96916;
double r96925 = r96924 / r96919;
double r96926 = pow(r96923, r96925);
double r96927 = r96918 * r96926;
return r96927;
}
double f(double k, double n) {
double r96928 = 1.0;
double r96929 = k;
double r96930 = sqrt(r96929);
double r96931 = r96928 / r96930;
double r96932 = 2.0;
double r96933 = r96928 - r96929;
double r96934 = r96933 / r96932;
double r96935 = pow(r96932, r96934);
double r96936 = atan2(1.0, 0.0);
double r96937 = pow(r96936, r96934);
double r96938 = r96935 * r96937;
double r96939 = r96931 * r96938;
double r96940 = n;
double r96941 = pow(r96940, r96934);
double r96942 = r96939 * r96941;
return r96942;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
rmApplied unpow-prod-down0.5
Applied associate-*r*0.5
rmApplied unpow-prod-down0.5
Final simplification0.5
herbie shell --seed 2019325 +o rules:numerics
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))