\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\left(\left(\frac{1}{\sqrt{k}} \cdot {2}^{\left(\frac{1 - k}{2}\right)}\right) \cdot {\pi}^{\left(\frac{1 - k}{2}\right)}\right) \cdot {n}^{\left(\frac{1 - k}{2}\right)}double f(double k, double n) {
double r4564505 = 1.0;
double r4564506 = k;
double r4564507 = sqrt(r4564506);
double r4564508 = r4564505 / r4564507;
double r4564509 = 2.0;
double r4564510 = atan2(1.0, 0.0);
double r4564511 = r4564509 * r4564510;
double r4564512 = n;
double r4564513 = r4564511 * r4564512;
double r4564514 = r4564505 - r4564506;
double r4564515 = r4564514 / r4564509;
double r4564516 = pow(r4564513, r4564515);
double r4564517 = r4564508 * r4564516;
return r4564517;
}
double f(double k, double n) {
double r4564518 = 1.0;
double r4564519 = k;
double r4564520 = sqrt(r4564519);
double r4564521 = r4564518 / r4564520;
double r4564522 = 2.0;
double r4564523 = r4564518 - r4564519;
double r4564524 = r4564523 / r4564522;
double r4564525 = pow(r4564522, r4564524);
double r4564526 = r4564521 * r4564525;
double r4564527 = atan2(1.0, 0.0);
double r4564528 = pow(r4564527, r4564524);
double r4564529 = r4564526 * r4564528;
double r4564530 = n;
double r4564531 = pow(r4564530, r4564524);
double r4564532 = r4564529 * r4564531;
return r4564532;
}



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.4
Applied associate-*r*0.4
Final simplification0.4
herbie shell --seed 2019174
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1.0 (sqrt k)) (pow (* (* 2.0 PI) n) (/ (- 1.0 k) 2.0))))