How do I compute the volume of the intersection of two $n$-dimensional ellipsoids?
Given an $n$-vector $c$ and a symmetric positive-definite $n\times n$ matrix $A$, define the ellipsoid $$E(c,A)=\{x|(A(x-c),x-c)<1\}$$ where $(\cdot,\cdot)$ if the dot product.
Then $$\mathrm{vol}(E(c,A))=\frac{u}{\det A}$$ where $u=\mathrm{vol}(E(0,1))$ is the volume of the unit sphere.
The question is: how do I compute the volume of the intersection $$\mathrm{vol}(E(c_1,A_1)\cap E(c_2,A_2))$$
I am more interested in being able to compute something relevant reasonably fast than in the exact correctness of the value. E.g., I would be happy to use the volumes of parallelepipeds (and their intersections) instead of ellipsoids.
EDIT: another acceptable alternative would be to define normal densities $f_1$ and $f_2$ (with mean $c_i$ and covariance $A_i$). What is $\int_{\Bbb{R}^n}f_1 f_2$? Something ugly, alas.