【求助】请教一个费米多重积分的算法！

horizonhp

我是学物理的，经常碰到费米分布的积分。
费米分布就是行如f(x)/（1+exp(x))函数的积分，如果是一维的，这个积分很好截断。
如果象这样：f(x, y, z)*(1/(1+exp(x))*(1-1/(1+exp(y))*(1/(1+exp(z)))
而且x, y, z的积分限是从负无穷积到正无穷。
请问这样的积分应该用怎样的数值方法得到数值呢？
  
最好能给出几个源程序，谢谢了！