double f(double i, double n) {
double r40090984 = 100.0;
double r40090985 = 1.0;
double r40090986 = i;
double r40090987 = n;
double r40090988 = r40090986 / r40090987;
double r40090989 = r40090985 + r40090988;
double r40090990 = pow(r40090989, r40090987);
double r40090991 = r40090990 - r40090985;
double r40090992 = r40090991 / r40090988;
double r40090993 = r40090984 * r40090992;
return r40090993;
}
double f(double i, double n) {
double r40090994 = i;
double r40090995 = -30308.320454546014;
bool r40090996 = r40090994 <= r40090995;
double r40090997 = n;
double r40090998 = r40090997 / r40090994;
double r40090999 = 1.0;
double r40091000 = r40090999 / r40090998;
double r40091001 = pow(r40091000, r40090997);
double r40091002 = r40091001 - r40090999;
double r40091003 = r40090998 * r40091002;
double r40091004 = 100.0;
double r40091005 = r40091003 * r40091004;
double r40091006 = 1.5640068539209153;
bool r40091007 = r40090994 <= r40091006;
double r40091008 = 16.666666666666668;
double r40091009 = r40090994 * r40091008;
double r40091010 = 50.0;
double r40091011 = r40091009 + r40091010;
double r40091012 = r40090994 * r40091011;
double r40091013 = r40091012 + r40091004;
double r40091014 = r40090997 * r40091013;
double r40091015 = r40090994 / r40090997;
double r40091016 = r40090999 + r40091015;
double r40091017 = pow(r40091016, r40090997);
double r40091018 = r40091017 - r40090999;
double r40091019 = exp(r40091018);
double r40091020 = log(r40091019);
double r40091021 = r40091020 / r40091015;
double r40091022 = r40091004 * r40091021;
double r40091023 = r40091007 ? r40091014 : r40091022;
double r40091024 = r40090996 ? r40091005 : r40091023;
return r40091024;
}
100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\begin{array}{l}
\mathbf{if}\;i \le -30308.320454546014:\\
\;\;\;\;\left(\frac{n}{i} \cdot \left({\left(\frac{1}{\frac{n}{i}}\right)}^{n} - 1\right)\right) \cdot 100\\
\mathbf{elif}\;i \le 1.5640068539209153:\\
\;\;\;\;n \cdot \left(i \cdot \left(i \cdot \frac{50}{3} + 50\right) + 100\right)\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{\log \left(e^{{\left(1 + \frac{i}{n}\right)}^{n} - 1}\right)}{\frac{i}{n}}\\
\end{array}



Bits error versus i




Bits error versus n
| Original | 42.6 |
|---|---|
| Target | 42.2 |
| Herbie | 19.3 |
if i < -30308.320454546014Initial program 27.7
Taylor expanded around inf 62.9
Simplified19.1
if -30308.320454546014 < i < 1.5640068539209153Initial program 50.0
Taylor expanded around 0 33.2
Simplified33.2
Taylor expanded around 0 17.2
Simplified17.2
if 1.5640068539209153 < i Initial program 30.8
rmApplied add-log-exp30.8
Applied add-log-exp31.0
Applied diff-log31.0
Simplified31.0
Final simplification19.3
herbie shell --seed 2019102
(FPCore (i n)
:name "Compound Interest"
:herbie-target
(* 100 (/ (- (exp (* n (if (== (+ 1 (/ i n)) 1) (/ i n) (/ (* (/ i n) (log (+ 1 (/ i n)))) (- (+ (/ i n) 1) 1))))) 1) (/ i n)))
(* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))))