Eu sei que fiz a mesma pergunta antes, mas como sou bastante nova aqui, a pergunta foi mal e não reproduzível. Portanto, eu tento fazer melhor aqui. (Se eu editar apenas o antigo, provavelmente ninguém o lerá)

Eu tenho essa integral dupla que eu gostaria de integrar:Aqui está uma foto

ff<-function(g,t) exp((16)*g)*exp(-8*t-(-t-0.01458757)^2/(0.0001126501))

integrate(Vectorize(function(t) integrate(function(g) 
                                          ff(g,t), -2.5,0)$value), -2, 2)

Executar isso em R me dá o erro:

  the integral is probably divergent

Quando tento executar a função sam no Wolfram, ele me fornece um valor adequado: (eu tive que alternar g = x et = y)

Ligação:

http://www.wolframalpha.com/input/?i=integration+[%2F%2Fmath%3Aexp%28%2816%29*x%29*exp%28-8*y-%28-y-0.01458757%29 ^ 2% 2F% 280.0001126501% 29% 29% 2F% 2F] + [% 2F% 2Fmath% 3Adx + dy% 2F% 2F] + para + x + de + [% 2F% 2Fmath% 3A-2,5% 2F% 2F] + a + [% 2F% 2Fmath% 3A0% 2F% 2F] + para + y + de + [% 2F% 2Fmath% 3A-2% 2F% 2F] + a + [% 2F% 2Fmath% 3A2% 2F% 2F]

Como você pode ver, obtém um resultado finito, alguém pode me ajudar aqui?

Plotei a função na área definida e não consegui encontrar um problema de singularidade. Vejo:

library('Plot3D')
x <- seq(-2.5,0, by = 0.01) #to see the peak change to: seq(-0.2,0, by = 0.001)
y <- seq(-2,2, by = 0.01) #"": seq(-0.1,0.1, by = 0.001)
grid <- mesh(x,y) 
z <- with(grid,exp((16)*x)*
  exp(-8*y-(-0.013615734-y-0.001+0.5*0.007505^2*1)^2/(2*0.007505^2)))
persp3D(z = z, x = x, y = y)

Obrigado por sua ajuda e espero que a questão esteja melhor estruturada do que a antiga.