Я знаю, что задавал тот же вопрос раньше, но, поскольку я здесь довольно нов, вопрос был задан плохо и не воспроизводим. Поэтому я стараюсь сделать это лучше здесь. (Если я отредактирую только старую, вероятно, никто не будет читать это

У меня есть этот двойной интеграл, который я хотел бы интегрировать:Вот картинка

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)

Запуск этого в R дает мне ошибку:

  the integral is probably divergent

Когда я пытаюсь запустить функцию sam в Wolfram, она дает мне правильное значение: (мне пришлось переключить g = x и t = y)

Ссылка на сайт:

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 + ау% 2F% 2F] + для + х + от + [% 2F% 2Fmath% 3A-2,5% 2F% 2F] + до + [% 2F% 2Fmath% 3A0% 2F% 2F] + при + у + с + [% 2F% 2Fmath% 3A2% 2F% 2F] + до + [% 2F% 2Fmath% 3A2% 2F% 2F]

Как вы видите, это дает конечный результат, может кто-нибудь помочь мне здесь?

Я построил функцию на определенной области и не смог найти проблему с особенностями. увидеть:

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)

Спасибо за вашу помощь, и я надеюсь, что вопрос лучше структурирован, чем старый.