После решения многих проблем с приведенным ниже интегралом мне все еще нужна помощь, чтобы найти правильное решение. Старался максимально упростить, спасибо за помощь.
Возьмем этот пример:
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, и она такая же, как здесь. У кого-нибудь есть решение для решения этой проблемы? большое спасибо!
integration [//math:exp(75*x)*exp(-35*y-((-y + 0.001283464)^2)/0.0001468898)//] [//math:dx dy//] for x from [//math:-2.5//] to [//math:min(0,y)//] for y from [//math:-2//] to [//math:2//]
- person Valegard234   schedule 19.09.2014