\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{1 \cdot {\left(\left(\left(2 \cdot \pi\right) \cdot \sqrt{n}\right) \cdot \sqrt{n}\right)}^{\left(\frac{1}{2}\right)}}{\sqrt{k} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{k}{2}\right)}}double f(double k, double n) {
double r130819 = 1.0;
double r130820 = k;
double r130821 = sqrt(r130820);
double r130822 = r130819 / r130821;
double r130823 = 2.0;
double r130824 = atan2(1.0, 0.0);
double r130825 = r130823 * r130824;
double r130826 = n;
double r130827 = r130825 * r130826;
double r130828 = r130819 - r130820;
double r130829 = r130828 / r130823;
double r130830 = pow(r130827, r130829);
double r130831 = r130822 * r130830;
return r130831;
}
double f(double k, double n) {
double r130832 = 1.0;
double r130833 = 2.0;
double r130834 = atan2(1.0, 0.0);
double r130835 = r130833 * r130834;
double r130836 = n;
double r130837 = sqrt(r130836);
double r130838 = r130835 * r130837;
double r130839 = r130838 * r130837;
double r130840 = r130832 / r130833;
double r130841 = pow(r130839, r130840);
double r130842 = r130832 * r130841;
double r130843 = k;
double r130844 = sqrt(r130843);
double r130845 = r130835 * r130836;
double r130846 = r130843 / r130833;
double r130847 = pow(r130845, r130846);
double r130848 = r130844 * r130847;
double r130849 = r130842 / r130848;
return r130849;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.5
rmApplied div-sub0.5
Applied pow-sub0.5
Applied frac-times0.4
rmApplied add-sqr-sqrt0.5
Applied associate-*r*0.5
Final simplification0.5
herbie shell --seed 2020100
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))