Решение для неизвестной верхней границы интеграла

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