\frac{p - {p}^{N}}{p - 1}\frac{p - {p}^{N}}{p - 1}double code(double p, double N) {
return ((double) (((double) (p - ((double) pow(p, N)))) / ((double) (p - 1.0))));
}
double code(double p, double N) {
return ((double) (((double) (p - ((double) pow(p, N)))) / ((double) (p - 1.0))));
}



Bits error versus p



Bits error versus N
Results
Initial program 0.0
Final simplification0.0
herbie shell --seed 2020152
(FPCore (p N)
:name "(/ (- p (pow p N)) (- p 1))"
:precision binary64
(/ (- p (pow p N)) (- p 1.0)))