\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\left({\left(n \cdot \pi\right)}^{\left(\frac{1 - k}{2}\right)} \cdot {2}^{\left(\frac{1 - k}{2}\right)}\right) \cdot \frac{1}{\sqrt{k}}double f(double k, double n) {
double r4628070 = 1.0;
double r4628071 = k;
double r4628072 = sqrt(r4628071);
double r4628073 = r4628070 / r4628072;
double r4628074 = 2.0;
double r4628075 = atan2(1.0, 0.0);
double r4628076 = r4628074 * r4628075;
double r4628077 = n;
double r4628078 = r4628076 * r4628077;
double r4628079 = r4628070 - r4628071;
double r4628080 = r4628079 / r4628074;
double r4628081 = pow(r4628078, r4628080);
double r4628082 = r4628073 * r4628081;
return r4628082;
}
double f(double k, double n) {
double r4628083 = n;
double r4628084 = atan2(1.0, 0.0);
double r4628085 = r4628083 * r4628084;
double r4628086 = 1.0;
double r4628087 = k;
double r4628088 = r4628086 - r4628087;
double r4628089 = 2.0;
double r4628090 = r4628088 / r4628089;
double r4628091 = pow(r4628085, r4628090);
double r4628092 = pow(r4628089, r4628090);
double r4628093 = r4628091 * r4628092;
double r4628094 = sqrt(r4628087);
double r4628095 = r4628086 / r4628094;
double r4628096 = r4628093 * r4628095;
return r4628096;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
Taylor expanded around 0 0.4
rmApplied add-sqr-sqrt0.5
Applied associate-*r*0.5
rmApplied unpow-prod-down0.4
Simplified0.5
Final simplification0.5
herbie shell --seed 2019172
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1.0 (sqrt k)) (pow (* (* 2.0 PI) n) (/ (- 1.0 k) 2.0))))