\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 r81669 = 1.0;
double r81670 = k;
double r81671 = sqrt(r81670);
double r81672 = r81669 / r81671;
double r81673 = 2.0;
double r81674 = atan2(1.0, 0.0);
double r81675 = r81673 * r81674;
double r81676 = n;
double r81677 = r81675 * r81676;
double r81678 = r81669 - r81670;
double r81679 = r81678 / r81673;
double r81680 = pow(r81677, r81679);
double r81681 = r81672 * r81680;
return r81681;
}
double f(double k, double n) {
double r81682 = 1.0;
double r81683 = k;
double r81684 = sqrt(r81683);
double r81685 = r81682 / r81684;
double r81686 = 2.0;
double r81687 = r81682 - r81683;
double r81688 = r81687 / r81686;
double r81689 = pow(r81686, r81688);
double r81690 = atan2(1.0, 0.0);
double r81691 = pow(r81690, r81688);
double r81692 = r81689 * r81691;
double r81693 = r81685 * r81692;
double r81694 = n;
double r81695 = pow(r81694, r81688);
double r81696 = r81693 * r81695;
return r81696;
}



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