Problem 235


に対し、

とする。
s(5000) = -600,000,000,000
となるrを求めよ。答えは小数点以下12桁に丸めること。

あんまり考えてない。

u <- function(k,r) (900-3*k)*r^(k-1)
s <- function(n,r) sum(u(1:n,r))
S <- function(r) s(5000,r) + 6e11
sprintf("%.12f", uniroot(S, interval=c(1,1.005), tol=1e-15)$root)