3

R積分の境界極限を解くことを含むアルゴリズムを実装しようとしています。たとえばa、次の積分が与えられた場合に見つけたいと思います。

integral_0^a exp(x) = 1/2

でそれを行う方法の大まかな考えがありますmatlab。しかし、これをどのように解決するのでしょうRか?

ご提案いただきありがとうございます。

4

1 に答える 1

5

を使用integrateして、積分を (数値的に) 計算しuniroot、方程式を (数値的に) 解くことができます。

f <- function(a) integrate( exp, 0, a )$value - 1/2
uniroot( f, c(-1, 1) ) # Look for a solution in [-1,1]
log(3/2) # Compare with the exact solution
于 2012-04-16T03:59:06.267 に答える