为什么我的双积分是R不工作

Question

我在这里看到这篇文章double integral in R并切换功能和限制以下匹配，但它不工作。

InnerFunc = function(x) { x + (y^2) } 
InnerIntegral = function(z) { sapply(y, 
     function(z) { integrate(InnerFunc, 0, 2)$value }) } 
integrate(InnerIntegral, 0, 1) 

我得到这个错误：

错误集成（InnerFunc，0,2）：功能 评价给出了错误类型的结果

Answer 1

你的变量都是出于怪事。这应该做你想做的

InnerFunc <- function(x, y) { x + (y^2) } 
InnerIntegral <- function(y) { sapply(y, 
    function(z) { integrate(InnerFunc, 0, 2, y=z)$value }) } 
integrate(InnerIntegral, 0, 1) 
# 2.666667 with absolute error < 3e-14 

通知InnerFunc现在怎么样既x和y的功能。另请注意，我们如何应用通过y参数传递给InnerIntegral的向量（而不是忽略传入的z值）。我们还将当前值y传递给InnerFunc。

虽然您键入

InnerFunc <- function(x, y) { x + (y^2) } 

您绘制的数学看起来确实应该引起

所以我不知道你真正追求的是什么。