\left(\left(-2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)}\right)}^{2}}\mathsf{hypot}\left(\sqrt{1}, {\left(\frac{U}{\left(J \cdot 2\right) \cdot \cos \left(\frac{K}{2}\right)}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(\cos \left(\frac{K}{2}\right) \cdot \left(-2 \cdot J\right)\right)double f(double J, double K, double U) {
double r137984 = -2.0;
double r137985 = J;
double r137986 = r137984 * r137985;
double r137987 = K;
double r137988 = 2.0;
double r137989 = r137987 / r137988;
double r137990 = cos(r137989);
double r137991 = r137986 * r137990;
double r137992 = 1.0;
double r137993 = U;
double r137994 = r137988 * r137985;
double r137995 = r137994 * r137990;
double r137996 = r137993 / r137995;
double r137997 = pow(r137996, r137988);
double r137998 = r137992 + r137997;
double r137999 = sqrt(r137998);
double r138000 = r137991 * r137999;
return r138000;
}
double f(double J, double K, double U) {
double r138001 = 1.0;
double r138002 = sqrt(r138001);
double r138003 = U;
double r138004 = J;
double r138005 = 2.0;
double r138006 = r138004 * r138005;
double r138007 = K;
double r138008 = r138007 / r138005;
double r138009 = cos(r138008);
double r138010 = r138006 * r138009;
double r138011 = r138003 / r138010;
double r138012 = 2.0;
double r138013 = r138005 / r138012;
double r138014 = pow(r138011, r138013);
double r138015 = hypot(r138002, r138014);
double r138016 = -2.0;
double r138017 = r138016 * r138004;
double r138018 = r138009 * r138017;
double r138019 = r138015 * r138018;
return r138019;
}



Bits error versus J



Bits error versus K



Bits error versus U
Results
Initial program 18.2
rmApplied sqr-pow18.2
Applied add-sqr-sqrt18.2
Applied hypot-def8.1
Final simplification8.1
herbie shell --seed 2019179 +o rules:numerics
(FPCore (J K U)
:name "Maksimov and Kolovsky, Equation (3)"
(* (* (* -2.0 J) (cos (/ K 2.0))) (sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) (cos (/ K 2.0)))) 2.0)))))