The result for any positive integer n
is the sum of two terms: (1) an Exp[x y]
term multiplied by a polynomial in x
and y
and (2) a polynomial in x
and `y'. The first term is
Exp[x y] LaguerreL[n, (x - y)^2]/n!
but I haven't figured out a general form for the second polynomial (assuming there is a closed form for that polynomial). So far for that second polynomial the sum of the x
and y
exponents is always even and the coefficient for x^i y^j
is the same as the coefficient for x^j y^i
. That's all I've determined. Maybe someone else can see the pattern.