wolfram과 비교하여 일관성없는 결과가 두 배 적분

아래의 적분과 관련된 많은 문제를 해결 한 후에도 적절한 해결책을 얻으려면 도움이 필요합니다. 나는 당신의 도움 덕분에 가능한 복잡성을 줄이려고 노력했다.wolfram과 비교하여 일관성없는 결과가 두 배 적분

b <- function(t) { # variable upper bound 
    if (t<0) t 
    else 0 
}

Doubleintegral :

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

동안 만 0.00012 ... Link to Wolfram

Wolfram에서 함수를 확인한 결과는 여기와 같습니다. 누군가이 문제를 해결할 수있는 해결책이 있습니까? 대단히 감사합니다!

출처

2014-09-19 Valegard234

W | A 링크를 클릭하면 "표준 계산 시간 초과"라는 메시지가 나타납니다. 의심 스럽네. –

PS. 익명 downvoter, 제발 이유를 설명하십시오. –

링크가 적합합니다. 이것을 과거의 wolfram에 복사하십시오 :'integration [// math : exp (75 * x) * exp (-35 * y - ((- y + 0.001283464)^2) /0.0001468898 // // [ [//math:-2.5//]에서 [// math : min (0, y) //]까지의 x에 대한 [dx dy //] [// math : -2 //]에서 [// 수학 : 2 //]' – Valegard234

계산 시간을 늘리면 결과가 일관됩니다. thanks @Hong Ooi

출처

2014-09-19 13:26:17 Valegard234

wolfram과 비교하여 일관성없는 결과가 두 배 적분

답변

관련 문제