\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{1 \cdot {\left(\frac{1}{k}\right)}^{\frac{1}{4}}}{\sqrt{\sqrt{k}}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}double f(double k, double n) {
double r125782 = 1.0;
double r125783 = k;
double r125784 = sqrt(r125783);
double r125785 = r125782 / r125784;
double r125786 = 2.0;
double r125787 = atan2(1.0, 0.0);
double r125788 = r125786 * r125787;
double r125789 = n;
double r125790 = r125788 * r125789;
double r125791 = r125782 - r125783;
double r125792 = r125791 / r125786;
double r125793 = pow(r125790, r125792);
double r125794 = r125785 * r125793;
return r125794;
}
double f(double k, double n) {
double r125795 = 1.0;
double r125796 = 1.0;
double r125797 = k;
double r125798 = r125796 / r125797;
double r125799 = 0.25;
double r125800 = pow(r125798, r125799);
double r125801 = r125795 * r125800;
double r125802 = sqrt(r125797);
double r125803 = sqrt(r125802);
double r125804 = r125801 / r125803;
double r125805 = 2.0;
double r125806 = atan2(1.0, 0.0);
double r125807 = r125805 * r125806;
double r125808 = n;
double r125809 = r125807 * r125808;
double r125810 = r125795 - r125797;
double r125811 = r125810 / r125805;
double r125812 = pow(r125809, r125811);
double r125813 = r125804 * r125812;
return r125813;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
rmApplied add-sqr-sqrt0.4
Applied sqrt-prod0.5
Applied associate-/r*0.5
rmApplied add-sqr-sqrt0.5
Applied sqrt-prod0.5
Applied sqrt-prod0.6
Applied *-un-lft-identity0.6
Applied times-frac0.6
Taylor expanded around 0 0.4
Final simplification0.4
herbie shell --seed 2020045
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))