\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 r134422 = 1.0;
double r134423 = k;
double r134424 = sqrt(r134423);
double r134425 = r134422 / r134424;
double r134426 = 2.0;
double r134427 = atan2(1.0, 0.0);
double r134428 = r134426 * r134427;
double r134429 = n;
double r134430 = r134428 * r134429;
double r134431 = r134422 - r134423;
double r134432 = r134431 / r134426;
double r134433 = pow(r134430, r134432);
double r134434 = r134425 * r134433;
return r134434;
}
double f(double k, double n) {
double r134435 = 1.0;
double r134436 = k;
double r134437 = sqrt(r134436);
double r134438 = r134435 / r134437;
double r134439 = 2.0;
double r134440 = r134435 - r134436;
double r134441 = r134440 / r134439;
double r134442 = pow(r134439, r134441);
double r134443 = atan2(1.0, 0.0);
double r134444 = pow(r134443, r134441);
double r134445 = r134442 * r134444;
double r134446 = r134438 * r134445;
double r134447 = n;
double r134448 = pow(r134447, r134441);
double r134449 = r134446 * r134448;
return r134449;
}



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 2019356 +o rules:numerics
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))