I want to calculate such an integral
$$
I=\int_0^{\infty}\int_0^{\infty}dxdy\ x^ny^m\exp(-ax^2-by^2-gxy+c_1x+c_2y),
$$
where $a,b,g>0$ are positive and real, $n,m$ are positive integers and $c_1,c_2$ are arbitrary complex numbers. For sure I can firstly integrate over, say $y$, which is
$$
\int_{0}^{\infty}dy\ y^m\exp(-by^2+(c_2-gx)y)\sim F(1+\frac{m}{2},\frac{3}{2},\frac{(c_2-gx)^2}{4b})
$$
where $F(\alpha,\gamma,z)$ is the Kummer Hypergeometric function. Now I havo to deal with the integral like this
$$
\int_0^{\infty}dx\ x^n\exp(-ax^2+c_1x)F(\alpha,\gamma,\frac{(c_2-gx)^2}{4b}).
$$
How can I go on with it?