https://github.com/cran/pracma
Raw File
Tip revision: c1688b374d201c13fb40b4dda2d2a89e34b94ec6 authored by Hans W. Borchers on 23 January 2021, 09:10:02 UTC
version 2.3.3
Tip revision: c1688b3
sumalt.R
## 
##  s u m a l t . R  Summing alternating series
##


sumalt <- function(f_alt, n) {
    b <- 2^(2*n-1)
    c <- b  # ; s <- 0
    s <- 0.0
    for (k in (n-1):0) {
        t <- f_alt(k)
        s <- s + c*t
        b <- b * (2*k+1) * (k+1) / (2 * (n-k) * (n+k))
        c <- c + b
    }
    s <- s / c
    return(s)
}
back to top