\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\sqrt{\left|\frac{a + b}{a} \cdot \frac{1}{\frac{a}{a - b}}\right|}double code(double a, double b) {
return sqrt(fabs((((a * a) - (b * b)) / (a * a))));
}
double code(double a, double b) {
return sqrt(fabs((((a + b) / a) * (1.0 / (a / (a - b))))));
}



Bits error versus a



Bits error versus b
Results
Initial program 14.1
rmApplied difference-of-squares14.1
Applied times-frac0.0
rmApplied clear-num0.0
Final simplification0.0
herbie shell --seed 2020058
(FPCore (a b)
:name "Eccentricity of an ellipse"
:precision binary64
:pre (<= 0.0 b a 1)
(sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))