2012-04-16 3 views
3

Rのアルゴリズムを実装しようとしていますが、これは積分の境界限界を解くことです。例えば、私はaは、以下の積分与え探したい:未知の積分上限を解く

integral_0^a exp(x) = 1/2

私はmatlabでそれを行うにはどのように大まかなアイデアを持っています。しかし、これを解決するにはどうすればR

ご意見ありがとうございます。

答えて

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 
+0

ありがとうございました!これは私が必要とするものです。 – notrockstar

関連する問題