Хорошо, я напишу свой вопрос, так как в латентной форме, я надеюсь, вы не против.Предположим, я хочу оценить интеграл в Rstudio 3.5.1. Я смотрел на пример
Rcode начинается здесь:
f <- function (x) {2/3 * (x [1] + x [2] + x [3])} </p>
adaptIntegrate (f, lowerLimit = c (0, 0, 0), upperLimit = c (0.5, 0.5, 0.5))
R-код заканчивается здесь
Что отлично работает, но если я хочу оценить следующий интеграл
$ \ int_ {0} ^ {1} \ int_ {0} ^ {y_ {1}} (y_ {1} - y_ {2}) \ cdot \ frac {1} {y_ {1}} dy_ {2} dy {1} $
Я пытался ввести его в Rstudioкак
Rcode начинается здесь:
f <- function (y) {(y [1] -y [2] * (1 / y [1]))} </p>
adaptIntegrate (f, lowerLimit = c (0,0), upperLimit = c (0, y [1]))
R-код заканчивается здесь
Но каждый раз, когда я пытаюсь оценить его,это дает мне ошибку
"Ошибка в adaptIntegrate (f, lowerLimit = c (0, 0), upperLimit = c (0, y [1])): объект 'y' не найден" *
Теперь мне интересно, нельзя ли вычислить интеграл с не вещественным числом как бунапример,в этом примере переменная.
Буду признателен за любую помощь или ссылки на то, как оценить эти типы интегралов в Rstudio