Any positive definite symmetric matrix $M$ can be written as $M = A^tA$ with some invertible matrix $A$ with positive eigenvalues, and then the inequality becomes $\| Ax \| < 1$, and the set of solutions is the image of the unit ball $\| y \|<1$ under the inverse of the matrix $x = A^{-1} y$. The change of volume under linear maps is given by the determinant, so the volume of the ellipsoid is $c_n\det A^{-1} = \frac{c_n}{\det A} = \frac{c_n}{\sqrt{\det M}}$, where $c_n$ is the volume of the $n$-dimensional unit ball.
Oops, you don't assume that the matrix is symmetric. In that case, the ellipsoid is described by the equation $x^t N x$ with $N = \frac12 (M+M^t)$ symmetric and positive definite, so you get $ \frac{c_n}{\sqrt{\det N}}$