\left(\frac{n - 1}{2} + \frac{x \cdot \left(n \cdot n - 1\right)}{12}\right) + \frac{\left(\left(x \cdot x\right) \cdot x\right) \cdot \left(1 - \left(\left(n \cdot n\right) \cdot n\right) \cdot n\right)}{720}\left(\frac{n - 1}{2} + \frac{x \cdot \left(n \cdot n - 1\right)}{12}\right) + \frac{\left(\left(x \cdot x\right) \cdot x\right) \cdot \left(1 - \left(\left(n \cdot n\right) \cdot n\right) \cdot n\right)}{720}double code(double n, double x) {
return ((double) (((double) (((double) (((double) (n - 1.0)) / 2.0)) + ((double) (((double) (x * ((double) (((double) (n * n)) - 1.0)))) / 12.0)))) + ((double) (((double) (((double) (((double) (x * x)) * x)) * ((double) (1.0 - ((double) (((double) (((double) (n * n)) * n)) * n)))))) / 720.0))));
}
double code(double n, double x) {
return ((double) (((double) (((double) (((double) (n - 1.0)) / 2.0)) + ((double) (((double) (x * ((double) (((double) (n * n)) - 1.0)))) / 12.0)))) + ((double) (((double) (((double) (((double) (x * x)) * x)) * ((double) (1.0 - ((double) (((double) (((double) (n * n)) * n)) * n)))))) / 720.0))));
}



Bits error versus n



Bits error versus x
Results
Initial program 12.3
Final simplification12.3
herbie shell --seed 2020153
(FPCore (n x)
:name "(+ (+ (/ (- n 1) 2) (/ (* x (- (* n n) 1)) 12)) (/ (* (* (* x x) x) (- 1 (* (* (* n n) n) n))) 720))"
:precision binary64
(+ (+ (/ (- n 1.0) 2.0) (/ (* x (- (* n n) 1.0)) 12.0)) (/ (* (* (* x x) x) (- 1.0 (* (* (* n n) n) n))) 720.0)))