\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{1}{\sqrt{k}} \cdot {\left({\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{\sqrt{1} + \sqrt{k}}{1}\right)}\right)}^{\left(\frac{\sqrt{1} - \sqrt{k}}{2}\right)}double code(double k, double n) {
return ((double) (((double) (1.0 / ((double) sqrt(k)))) * ((double) pow(((double) (((double) (2.0 * ((double) M_PI))) * n)), ((double) (((double) (1.0 - k)) / 2.0))))));
}
double code(double k, double n) {
return ((double) (((double) (1.0 / ((double) sqrt(k)))) * ((double) pow(((double) pow(((double) (((double) (2.0 * ((double) M_PI))) * n)), ((double) (((double) (((double) sqrt(1.0)) + ((double) sqrt(k)))) / 1.0)))), ((double) (((double) (((double) sqrt(1.0)) - ((double) sqrt(k)))) / 2.0))))));
}



Bits error versus k



Bits error versus n
Results
Initial program 0.5
rmApplied *-un-lft-identity0.5
Applied add-sqr-sqrt0.6
Applied add-sqr-sqrt0.6
Applied difference-of-squares0.6
Applied times-frac0.6
Applied pow-unpow0.8
Final simplification0.8
herbie shell --seed 2020121
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))