[問題] 積分計算值

[slow0201 (海星)]

大家好,有個問題想請教

描述如下:
函數為兩個常態的CDF和一個常態的PDF的乘積
積分下限為負無窮大,上限為無窮大

請問直接使用下面程式積分是否容易產生錯誤結果?
我跑出來的結果比論文提供的值低很多
請問是否有別的寫法可以使得我得到的結果較接近
論文提供的結果,謝謝!
--------------------------------------------
程式編寫為
x=1
y=1.486607
z=1
w=1.486607

fun1<-function(u)
{
  fun_ev=pnorm(x*u+y)*pnorm(-z*u+w)*dnorm(u)
  return(fun_ev)
}
AS=integrate(fun1,lower=-Inf,upper=Inf)$value
----------------------------------------------

出來的結果AS是0.7096956
----------------------------------------------







--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.109.73.41
→ celestialgod:我幫你用MATLAB算也是0.7097
→ celestialgod:你應該去看一下你的FUNCTION也沒有寫錯
→ celestialgod:                          有
→ celestialgod:其他FUNCTION如quadinf (package:pracma)
→ celestialgod:AS2=quadinf(fun, -Inf, Inf, tol=1e-10);AS2=0.7097
→ slow0201:感謝樓上的幫忙,實在是不知道文章中的0.87是怎麼算出來的
→ psinqoo:sum(AS) [1] 0.7096957    summary(AS)   0.7097