求解积分的未知上界

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