Question

解决了以下积分的许多问题后，我仍需要一些帮助才能找到合适的解决方案。我试图尽可能地降低复杂性，谢谢你的帮助。

让我们举个例子：

b <- function(t) { # variable upper bound
  if (t<0) t
  else 0
}

Doubleintegral：

f<-function(g,t){
    exp(75*g)*
    exp(-35*t-((-t + 0.001283464)^2)/0.0001468898)
}

result <- (integrate(Vectorize(function(t) integrate(function(g) f(g,t), -2.5,b(t))$value), -2, 2, stop.on.error=FALSE)$value)

这给了我：

> (integrate(Vectorize(function(t) integrate(function(g) f(g,t), -2.5,b(t))$value), -2, 2, stop.on.error=FALSE)$value)
[1] 0.0002256195

虽然Wolfram Alpha只给我0.00012 .. Link to Wolfram

我检查了Wolfram中的函数，它与此处相同。有人有解决方案来解决这个问题吗？非常感谢你！

Answer 1

通过延长计算时间，结果是一致的。谢谢@Hong Ooi

与wolfram相比，双重积分与不一致的结果

1 个答案: