https://github.com/cran/pracma
Tip revision: c1688b374d201c13fb40b4dda2d2a89e34b94ec6 authored by Hans W. Borchers on 23 January 2021, 09:10:02 UTC
version 2.3.3
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)
}