100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\begin{array}{l}
\mathbf{if}\;i \leq -20604699981894388:\\
\;\;\;\;100 \cdot \frac{{\left(\frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{elif}\;i \leq 1.0526062084788294 \cdot 10^{-06}:\\
\;\;\;\;100 \cdot \left(n + n \cdot \left(i \cdot \left(0.5 + i \cdot 0.16666666666666666\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{100 \cdot \left({\left(\frac{i}{n} + 1\right)}^{n} - 1\right)}{\frac{i}{n}}\\
\end{array}(FPCore (i n) :precision binary64 (* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
(FPCore (i n)
:precision binary64
(if (<= i -20604699981894388.0)
(* 100.0 (/ (- (pow (/ i n) n) 1.0) (/ i n)))
(if (<= i 1.0526062084788294e-06)
(* 100.0 (+ n (* n (* i (+ 0.5 (* i 0.16666666666666666))))))
(/ (* 100.0 (- (pow (+ (/ i n) 1.0) n) 1.0)) (/ i n)))))double code(double i, double n) {
return 100.0 * ((pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
double code(double i, double n) {
double tmp;
if (i <= -20604699981894388.0) {
tmp = 100.0 * ((pow((i / n), n) - 1.0) / (i / n));
} else if (i <= 1.0526062084788294e-06) {
tmp = 100.0 * (n + (n * (i * (0.5 + (i * 0.16666666666666666)))));
} else {
tmp = (100.0 * (pow(((i / n) + 1.0), n) - 1.0)) / (i / n);
}
return tmp;
}




Bits error versus i




Bits error versus n
Results
| Original | 47.9 |
|---|---|
| Target | 47.5 |
| Herbie | 15.0 |
if i < -20604699981894388Initial program 27.4
Taylor expanded around inf 64.0
Simplified17.9
if -20604699981894388 < i < 1.0526062084788294e-6Initial program 58.1
Taylor expanded around 0 26.4
Simplified26.4
Taylor expanded around 0 10.3
Simplified10.3
if 1.0526062084788294e-6 < i Initial program 32.0
rmApplied associate-*r/_binary64_274832.0
Simplified32.0
Final simplification15.0
herbie shell --seed 2020343
(FPCore (i n)
:name "Compound Interest"
:precision binary64
:herbie-target
(* 100.0 (/ (- (exp (* n (if (== (+ 1.0 (/ i n)) 1.0) (/ i n) (/ (* (/ i n) (log (+ 1.0 (/ i n)))) (- (+ (/ i n) 1.0) 1.0))))) 1.0) (/ i n)))
(* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))