R中几何级数的计算
我想计算一下:R中几何级数的计算,r,R,我想计算一下: $\sum_{j=1}^n r^j$ for $n=10, 20, 30, 40$, where $r=1.08$ 把答案看成一个向量。我和R一起工作;有人能给我解释一下怎么做吗?那么你正在玩power series。你可以: r <- 1.08 ## this will be a divergent series, toward `Inf` rr <- r^(1:40) cumsum(rr)[1:4 * 10] # [1] 15.64549 49.42292
$\sum_{j=1}^n r^j$ for $n=10, 20, 30, 40$, where $r=1.08$
把答案看成一个向量。我和R一起工作;有人能给我解释一下怎么做吗?那么你正在玩power series。你可以:
r <- 1.08 ## this will be a divergent series, toward `Inf`
rr <- r^(1:40)
cumsum(rr)[1:4 * 10]
# [1] 15.64549 49.42292 122.34587 279.78104
r也许(如果我正确解释了latex):
库(purrr)
rr的值未说明,我猜所要求的公式并不是复利问题的正确解决方案,但以下是对目前提出的问题的一个答案:
r = 5.3; for (n in c(10,20,30,40) ) print( sum( r^(1:n) ) )
[1] 21555896
[1] 3.769856e+14
[1] 6.593006e+21
[1] 1.153034e+29
仅用于循环并将输出作为向量:
out=c();
for(i in c(10,20,30,40))
{sum=0;
for(j in 1:i)
{sum=sum+1.08^j};
out=c(out,sum)};
out
输出:
[1] 15.64549 49.42292 122.34587 279.78104
或:
或者,以R的方式:
> sapply(seq(10,40,10), function(x) sum(1.08^(1:x)))
[1] 15.64549 49.42292 122.34587 279.78104
什么是r?常数?对于原子随机正态r
,sapply(seq(10,40,10),函数(x){sum(rnorm(1)^seq(x))}
Fyi,r有文档。例如,您可以键入?sum
或?sequence
或?sequence
。答案不是问题中所问的“作为向量”。如果这困扰您,请尝试cumsum(r^(1:n))[seq(10,40,by=10)]
[1] 15.64549 49.42292 122.34587 279.78104
> out=c();
> for(i in c(10,20,30,40))
+ out=c(out,sum(1.08^(1:i)))
> out
[1] 15.64549 49.42292 122.34587 279.78104
> sapply(seq(10,40,10), function(x) sum(1.08^(1:x)))
[1] 15.64549 49.42292 122.34587 279.78104