TaylorSeries[ Exp[x^2 + y^2] Exp[ - (x^2 - 2 u x y + y^2)/(1-u^2)]/ Sqrt[(1-u^2)], {u, 0, 5}]