exp(n sum_j (-1)^(j+1) a^j / j), the product expansion of (1+a)^n.
Usage
PowLogSer(a, n, terms = 12)
Arguments
- a
perturbation with abs(a) < 1, so the log series converges.
- n
exponent.
- terms
number of log-series terms, >= 1.
Value
list(exact, product_form, rel_error).
References
Morin, D. J. (2016). Probability: For the Enthusiastic
Beginner. Createspace. Eq (7.21).
Examples
PowLogSer(0.1, 20)$rel_error
#> [1] 1.421878e-13
