18/06/29 10:52:22.56 zpGshx3p.net
>>281
理論値は
> collector <- function (n){
+ Hn=0
+ for (i in 1:n) Hn=Hn + 1/i
+ n*Hn
+ }
> collector(5)
[1] 11.41667
シミュレーションの値とほぼ一致しているのでどちらの計算も正しいのだろうな。
18/06/29 10:52:22.56 zpGshx3p.net
>>281
理論値は
> collector <- function (n){
+ Hn=0
+ for (i in 1:n) Hn=Hn + 1/i
+ n*Hn
+ }
> collector(5)
[1] 11.41667
シミュレーションの値とほぼ一致しているのでどちらの計算も正しいのだろうな。