求解積分的未知上界

Question

我試圖在R中實現一個算法，該算法涉及求解積分的邊界限制。例如，我想找到a給出以下組成：

integral_0^a exp(x) = 1/2

我有一個粗略的想法如何做到這一點的matlab。但是如何解決這個問題呢？R？

感謝您的建議。

Answer 1

您可以使用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