Двойной интеграл с противоречивыми результатами по сравнению с вольфрамом

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

Возьмем этот пример:

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, и она такая же, как здесь. У кого-нибудь есть решение для решения этой проблемы? большое спасибо!


person Valegard234    schedule 19.09.2014    source источник
comment
Когда я нажимаю на вашу ссылку W|A, я получаю сообщение о превышении стандартного времени вычисления. Звучит подозрительно.   -  person Hong Ooi    schedule 19.09.2014
comment
PS. Анонимный downvoter, объясните свои причины, пожалуйста.   -  person Hong Ooi    schedule 19.09.2014
comment
Ссылка у меня работает. скопируйте это в вольфрам: 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
comment
Да, это дает числовой результат. Но вы также видели коробку внизу с этим предупреждающим сообщением? Возможно, вычисление было остановлено преждевременно.   -  person Hong Ooi    schedule 19.09.2014
comment
боже мой, я не видел этого сообщения, и у меня действительно есть про-аккаунт. Если я увеличу время вычисления, это даст мне тот же результат. спасибо за подсказку.   -  person Valegard234    schedule 19.09.2014
comment
Возможно, вы захотите указать это как ответ или отредактировать свой вопрос, чтобы отметить, что вы его решили.   -  person Hong Ooi    schedule 19.09.2014


Ответы (1)


Увеличивая время вычисления, результаты согласуются. спасибо @Hong Ooi

person Valegard234    schedule 19.09.2014