read MultiAlmkvistZeilberger:

print(MAZ(1,exp(-x^2/2-y^2/2-z^2/2),(x-y)^2*(x-z)^2*(y-z)^2,[x,y,z],n,N,{}));

time();

quit: