\frac{a - b}{{\left(\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}\right)}^{3}} - \frac{a + b}{{\left(\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}\right)}^{3}}\frac{a - b}{{\left(\sqrt{\left(a - b\right) \cdot \left(a - b\right) + r \cdot r}\right)}^{3}} - \frac{a + b}{{\left(\sqrt{\left(a + b\right) \cdot \left(a + b\right) + r \cdot r}\right)}^{3}}double code(double a, double b, double r) {
return ((double) (((double) (((double) (a - b)) / ((double) pow(((double) sqrt(((double) (((double) (((double) (a - b)) * ((double) (a - b)))) + ((double) (r * r)))))), 3.0)))) - ((double) (((double) (a + b)) / ((double) pow(((double) sqrt(((double) (((double) (((double) (a + b)) * ((double) (a + b)))) + ((double) (r * r)))))), 3.0))))));
}
double code(double a, double b, double r) {
return ((double) (((double) (((double) (a - b)) / ((double) pow(((double) sqrt(((double) (((double) (((double) (a - b)) * ((double) (a - b)))) + ((double) (r * r)))))), 3.0)))) - ((double) (((double) (a + b)) / ((double) pow(((double) sqrt(((double) (((double) (((double) (a + b)) * ((double) (a + b)))) + ((double) (r * r)))))), 3.0))))));
}



Bits error versus a



Bits error versus b



Bits error versus r
Results
Initial program 10.2
Final simplification10.2
herbie shell --seed 2020152
(FPCore (a b r)
:name "(- (/ (- a b) (pow (sqrt (+ (* (- a b) (- a b)) (* r r))) 3)) (/ (+ a b) (pow (sqrt (+ (* (+ a b) (+ a b)) (* r r))) 3)))"
:precision binary64
(- (/ (- a b) (pow (sqrt (+ (* (- a b) (- a b)) (* r r))) 3.0)) (/ (+ a b) (pow (sqrt (+ (* (+ a b) (+ a b)) (* r r))) 3.0))))