Divergente Integral en R es solucionable en Wolfram
Sé que hice la misma pregunta antes, pero como soy bastante nuevo aquí, la pregunta fue pobre y no reproducible. Por lo tanto, trato de hacerlo mejor aquí. (Si solo edito el viejo, probablemente nadie lo leerá)
Tengo esta doble integral que me gustaría integrar:Aquí hay una 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)
Ejecutar esto en R me da el error:
the integral is probably divergent
Cuando intento ejecutar la función sam en Wolfram, me da un valor adecuado: (tuve que cambiar g = xyt = y)
Enlazar:
Como puede ver, obtiene un resultado finito, ¿alguien puede ayudarme aquí?
Tracé la función en el área definida y no pude encontrar un problema de singularidad. ver:
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)
Gracias por su ayuda y espero que la pregunta esté mejor estructurada que la anterior.