解决积分未知上限问题

Question

我正试图在R中实现一个涉及求解积分边界限制的算法。例如，我想找到给定以下积分的a：

integral_0^a exp(x) = 1/2

我大致知道如何在matlab中完成它。但是在R中如何解决呢？

感谢您的建议。

- notrockstar

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

- Vincent Zoonekynd