Does the integral of PDF of multi-normal distribution over quarter planes have a closed form?

I am interested in finding a closed form solution (wich I suspect does not exist) to the following integral

$$\displaystyle \int _a^{\infty }\int _b^{\infty } \frac{\exp \left(-\frac{x^2+y^2-2 c x y}{2 \left(1-c^2\right)}\right)}{2 \pi \sqrt{1-c^2}} dy dx$$

which corresponds to the integral of the PDF$(x,y)$ of a multiNormalDistribution (of covariance coefficient $c$) over the quarter plane $x>a$ and $y>b$. Here $a$ and $b$ are positive and $0<c<1$ (and I know a solution exists for $a=b=0$, but this is not sufficient for my purpose).

More generally I would be interested in the $3$D generalization of this problem.

I have tried in Mathematica to no avail.

Regards,

Generically, cumulative distribution function of multivariate Gaussian vector is not expressible in terms of cdf of standard normal random variable $\Phi(x)$. The book by Alan Genz and Frank Brentz, "Computation of multivariate Normal and t Probabilities" is good reference on the subject.
For a standard 2D Gaussian vector $(X,Y)$ with correlation coefficient $-1 < \rho <1$, the probability $\mathbb{P}(X>a,Y>b)$ can be expressed in terms of Owen's T-function.
• @AdelBibi Yes, because $F_{X,Y}\left(x,y)\right) = \Pr\left( X \leqslant x, Y \leqslant y\right) = 1 - \Pr\left(X > x\right) - \Pr\left(Y>y\right) + \Pr\left(X>x, Y>y\right) = F_X\left(x\right) +F_Y\left(y\right) - \left(1- \Pr\left(X>x, Y>y\right)\right)$. The latter bivariate probability can be expressed in terms of Owen's T-function as I claimed in the post. Jul 4 '17 at 17:33