https://github.com/cran/pracma
Tip revision: 26e049d70b4a1c237987e260cba68f6a9413736c authored by Hans W. Borchers on 09 April 2019, 04:10:07 UTC
version 2.2.5
version 2.2.5
Tip revision: 26e049d
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)
}