用二重积分命令:

int(int('f(x,y)','x',-inf, 3),'y',-inf, 3)

其中f(x,y)是密度函数

[此贴子已经被作者于2009-3-28 9:55:21编辑过]

int(-1/2*erf((-b+w*z2*a)/a/(2-2*w)^(1/2)/b)*pi^(1/2)*a*(2-2*w)^(1/2)*exp(-1/2*z2^2/b^2*w)*exp(-1/2*z2^2/b^2)+1/2*erf(1/(2-2*w)^(1/2)*w*z2/b)*pi^(1/2)*a*(2-2*w)^(1/2)*exp(-1/2*z2^2/b^2*w)*exp(-1/2*z2^2/b^2),z2 = 1 .. 2)，这个好像是积不出来对吗？请问有没有其他数值算法可以把它积分出来，结果是具体一个近似值之类的？

应该是可以的,求出一个式子(ans) 来之后, 你用强制转换 double(ans)就可以了

谢谢啦！以后多交流