无限积分不一致？ （adaptIntegrate）

Question

我用过

library('cubature')
adaptIntegrate(doubleintegral, lowerLimit = c(-2.5, -2), upperLimit = c(0, 2), x=x,r=r,m=m,n=n)$integral

集成以下功能：

doubleintegral <- function(y,x,r,n,m){
  max(((-y[1])^(n-1)*(y[2]-y[1])^(m-1)*exp((x[7]+x[4])*y[1]))*
  exp(-x[4]*y[2] - (r[v]-y[2]-x[5]*x[3]+0.5*x[6]^2*x[3])^2/(2*x[6]^2*x[3])),0)}

我使用了这个示例参数

x <- numeric()
x[1] = 42
x[2] = 21
x[3] = 1
x[4] = 72.9
x[5] = 0.0332
x[6] = 0.0311
x[7] = 16.8
r <- numeric()
r = 0.0006
v = 1
r[v] = -0.036

w = 1
n = 2
i = 1
m = 2

我不明白以下输出是如何可能的：

> adaptIntegrate(doubleintegral, lowerLimit = c(-130, -4), upperLimit = c(0, 4), x=x,r=r,m=m,n=n)$integral
[1] 1.12184e-07
> adaptIntegrate(doubleintegral, lowerLimit = c(-2.5, -2), upperLimit = c(0, 2), x=x,r=r,m=m,n=n)$integral
[1] 2.516489e-07

通过扩大限制，我得到一个更小的值。我遇到了与功能相同的问题＆＃34;整合＆＃34;我申请了另一个功能。这怎么可能？

有人对这种双重整合有另一个想法吗？

感谢您的帮助！