r - 求解未知的积分上限

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