与 Wolfram 相比结果不一致的二重积分

Question

在解决了以下积分的许多问题之后，我仍然需要一些帮助才能找到正确的解决方案。我尝试尽可能降低复杂性，感谢您的帮助。

让我们举个例子：

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

二重积分：

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.. Wolfram 链接

我检查了Wolfram中的函数，它和这里的一样。有人有办法解决这个问题吗？非常感谢！

Answer 1

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