read MultiAlmkvistZeilberger:

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

time();

quit: