\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{{\left(\pi \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}}double f(double k, double n) {
double r38152792 = 1.0;
double r38152793 = k;
double r38152794 = sqrt(r38152793);
double r38152795 = r38152792 / r38152794;
double r38152796 = 2.0;
double r38152797 = atan2(1.0, 0.0);
double r38152798 = r38152796 * r38152797;
double r38152799 = n;
double r38152800 = r38152798 * r38152799;
double r38152801 = r38152792 - r38152793;
double r38152802 = r38152801 / r38152796;
double r38152803 = pow(r38152800, r38152802);
double r38152804 = r38152795 * r38152803;
return r38152804;
}
double f(double k, double n) {
double r38152805 = atan2(1.0, 0.0);
double r38152806 = n;
double r38152807 = 2.0;
double r38152808 = r38152806 * r38152807;
double r38152809 = r38152805 * r38152808;
double r38152810 = 1.0;
double r38152811 = k;
double r38152812 = r38152810 - r38152811;
double r38152813 = r38152812 / r38152807;
double r38152814 = pow(r38152809, r38152813);
double r38152815 = sqrt(r38152811);
double r38152816 = r38152814 / r38152815;
return r38152816;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
Simplified0.3
rmApplied clear-num0.4
rmApplied *-un-lft-identity0.4
Applied *-un-lft-identity0.4
Applied sqrt-prod0.4
Applied times-frac0.4
Applied *-un-lft-identity0.4
Applied times-frac0.4
Simplified0.4
Simplified0.3
Final simplification0.3
herbie shell --seed 2019128
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))