双重积分给出错误对象 y not found

Question

我正在尝试进行这种集成

我正在尝试这段代码

integrand_inner <- function(x, a, b) {(a * x^2 + b * y *x)*y}

integrand_outer <- Vectorize(
  function(y, a, b){
    integrate(integrand_inner,
              (bo - b)*y/(a-ao),
              1,
              a=a,
              b=b)$value
  }
)

revenue <- function(a){
  b <- sqrt(1-a^2)
  
  return(integrate(integrand_outer, 0, 1, a = a, b=b))
}

但我得到了错误

Error in f(x, ...) : object 'y' not found

我不知道出了什么问题。

Answer 1

需要将 y 作为参数添加到函数并传递。

integrand_inner <- function(x,y, a, b) {(a * x^2 + b * y *x)*y}

integrand_outer <- Vectorize(
  function(y, a, b){
    integrate(integrand_inner,
              (bo - b)*y/(a-ao),
              1,
              a=a,
              b=b,
              y=y)$value
  }
)

revenue <- function(a){
  b <- sqrt(1-a^2)
  
  return(integrate(integrand_outer, 0, 1, a = a, b=b))
}