\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 r137440 = 1.0;
double r137441 = k;
double r137442 = sqrt(r137441);
double r137443 = r137440 / r137442;
double r137444 = 2.0;
double r137445 = atan2(1.0, 0.0);
double r137446 = r137444 * r137445;
double r137447 = n;
double r137448 = r137446 * r137447;
double r137449 = r137440 - r137441;
double r137450 = r137449 / r137444;
double r137451 = pow(r137448, r137450);
double r137452 = r137443 * r137451;
return r137452;
}
double f(double k, double n) {
double r137453 = 1.0;
double r137454 = k;
double r137455 = sqrt(r137454);
double r137456 = r137453 / r137455;
double r137457 = 2.0;
double r137458 = r137453 - r137454;
double r137459 = r137458 / r137457;
double r137460 = pow(r137457, r137459);
double r137461 = atan2(1.0, 0.0);
double r137462 = pow(r137461, r137459);
double r137463 = r137460 * r137462;
double r137464 = r137456 * r137463;
double r137465 = n;
double r137466 = pow(r137465, r137459);
double r137467 = r137464 * r137466;
return r137467;
}



Bits error versus k



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