\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\left(\frac{\sqrt{2}}{\sqrt{k}} \cdot {\pi}^{\left(\frac{1 - k}{2}\right)}\right) \cdot \frac{{n}^{\left(\frac{1 - k}{2}\right)}}{{2}^{\left(\frac{k}{2}\right)}}double f(double k, double n) {
double r27720307 = 1.0;
double r27720308 = k;
double r27720309 = sqrt(r27720308);
double r27720310 = r27720307 / r27720309;
double r27720311 = 2.0;
double r27720312 = atan2(1.0, 0.0);
double r27720313 = r27720311 * r27720312;
double r27720314 = n;
double r27720315 = r27720313 * r27720314;
double r27720316 = r27720307 - r27720308;
double r27720317 = r27720316 / r27720311;
double r27720318 = pow(r27720315, r27720317);
double r27720319 = r27720310 * r27720318;
return r27720319;
}
double f(double k, double n) {
double r27720320 = 2.0;
double r27720321 = sqrt(r27720320);
double r27720322 = k;
double r27720323 = sqrt(r27720322);
double r27720324 = r27720321 / r27720323;
double r27720325 = atan2(1.0, 0.0);
double r27720326 = 1.0;
double r27720327 = r27720326 - r27720322;
double r27720328 = r27720327 / r27720320;
double r27720329 = pow(r27720325, r27720328);
double r27720330 = r27720324 * r27720329;
double r27720331 = n;
double r27720332 = pow(r27720331, r27720328);
double r27720333 = r27720322 / r27720320;
double r27720334 = pow(r27720320, r27720333);
double r27720335 = r27720332 / r27720334;
double r27720336 = r27720330 * r27720335;
return r27720336;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
Simplified0.3
rmApplied unpow-prod-down0.5
rmApplied unpow-prod-down0.5
rmApplied associate-/l*0.5
rmApplied div-sub0.5
Applied pow-sub0.5
Applied associate-*l/0.5
Applied associate-/r/0.5
Applied *-un-lft-identity0.5
Applied times-frac0.5
Simplified0.4
Final simplification0.4
herbie shell --seed 2019120
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))