\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\left(\frac{{\pi}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}{\sqrt{k}} \cdot {2}^{\left(\frac{1}{2} - \frac{k}{2}\right)}\right) \cdot {n}^{\left(\frac{1}{2} - \frac{k}{2}\right)}double f(double k, double n) {
double r4165912 = 1.0;
double r4165913 = k;
double r4165914 = sqrt(r4165913);
double r4165915 = r4165912 / r4165914;
double r4165916 = 2.0;
double r4165917 = atan2(1.0, 0.0);
double r4165918 = r4165916 * r4165917;
double r4165919 = n;
double r4165920 = r4165918 * r4165919;
double r4165921 = r4165912 - r4165913;
double r4165922 = r4165921 / r4165916;
double r4165923 = pow(r4165920, r4165922);
double r4165924 = r4165915 * r4165923;
return r4165924;
}
double f(double k, double n) {
double r4165925 = atan2(1.0, 0.0);
double r4165926 = 0.5;
double r4165927 = k;
double r4165928 = 2.0;
double r4165929 = r4165927 / r4165928;
double r4165930 = r4165926 - r4165929;
double r4165931 = pow(r4165925, r4165930);
double r4165932 = sqrt(r4165927);
double r4165933 = r4165931 / r4165932;
double r4165934 = pow(r4165928, r4165930);
double r4165935 = r4165933 * r4165934;
double r4165936 = n;
double r4165937 = pow(r4165936, r4165930);
double r4165938 = r4165935 * r4165937;
return r4165938;
}



Bits error versus k



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