\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{{\left(\left(\pi \cdot 2\right) \cdot n\right)}^{\left(\frac{\frac{1}{2} - \frac{k}{2}}{2}\right)}}{\frac{\sqrt{k}}{{\left({\left(\left(\pi \cdot 2\right) \cdot n\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}\right)}^{\frac{1}{2}}}}double f(double k, double n) {
double r5284117 = 1.0;
double r5284118 = k;
double r5284119 = sqrt(r5284118);
double r5284120 = r5284117 / r5284119;
double r5284121 = 2.0;
double r5284122 = atan2(1.0, 0.0);
double r5284123 = r5284121 * r5284122;
double r5284124 = n;
double r5284125 = r5284123 * r5284124;
double r5284126 = r5284117 - r5284118;
double r5284127 = r5284126 / r5284121;
double r5284128 = pow(r5284125, r5284127);
double r5284129 = r5284120 * r5284128;
return r5284129;
}
double f(double k, double n) {
double r5284130 = atan2(1.0, 0.0);
double r5284131 = 2.0;
double r5284132 = r5284130 * r5284131;
double r5284133 = n;
double r5284134 = r5284132 * r5284133;
double r5284135 = 0.5;
double r5284136 = k;
double r5284137 = r5284136 / r5284131;
double r5284138 = r5284135 - r5284137;
double r5284139 = r5284138 / r5284131;
double r5284140 = pow(r5284134, r5284139);
double r5284141 = sqrt(r5284136);
double r5284142 = pow(r5284134, r5284138);
double r5284143 = pow(r5284142, r5284135);
double r5284144 = r5284141 / r5284143;
double r5284145 = r5284140 / r5284144;
return r5284145;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
Simplified0.3
rmApplied sqr-pow0.4
Applied associate-/l*0.4
rmApplied div-inv0.4
Applied pow-unpow0.4
Final simplification0.4
herbie shell --seed 2019138
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))