\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{{\left(n \cdot \pi\right)}^{\left(\frac{1}{2} - k \cdot \frac{1}{2}\right)} \cdot {2}^{\left(\frac{1}{2} - k \cdot \frac{1}{2}\right)}}{\sqrt{k}}double f(double k, double n) {
double r4511535 = 1.0;
double r4511536 = k;
double r4511537 = sqrt(r4511536);
double r4511538 = r4511535 / r4511537;
double r4511539 = 2.0;
double r4511540 = atan2(1.0, 0.0);
double r4511541 = r4511539 * r4511540;
double r4511542 = n;
double r4511543 = r4511541 * r4511542;
double r4511544 = r4511535 - r4511536;
double r4511545 = r4511544 / r4511539;
double r4511546 = pow(r4511543, r4511545);
double r4511547 = r4511538 * r4511546;
return r4511547;
}
double f(double k, double n) {
double r4511548 = n;
double r4511549 = atan2(1.0, 0.0);
double r4511550 = r4511548 * r4511549;
double r4511551 = 0.5;
double r4511552 = k;
double r4511553 = r4511552 * r4511551;
double r4511554 = r4511551 - r4511553;
double r4511555 = pow(r4511550, r4511554);
double r4511556 = 2.0;
double r4511557 = pow(r4511556, r4511554);
double r4511558 = r4511555 * r4511557;
double r4511559 = sqrt(r4511552);
double r4511560 = r4511558 / r4511559;
return r4511560;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
Simplified0.3
rmApplied *-un-lft-identity0.3
Applied sqrt-prod0.3
Applied associate-/r*0.3
Simplified0.3
rmApplied unpow-prod-down0.4
Final simplification0.4
herbie shell --seed 2019158 +o rules:numerics
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))