Volume of $n$-dimensional parallelepiped as determinant
Solution 1:
This is essentially equivalent to showing that a linear operator represented by a matrix $A$ takes shapes to the domain to shapes in the range, and multiplies their volume by the determinant of $A$.
This is the way I convinced myself of this many years ago. First, it is obviously true for diagonal matrices. Next, it is also true for shear matrices, which are of the form identity matrix plus a matrix all of whose entries are zero except one entry. Finally, it is obviously true for permutation matrices. Now you merely need to prove that every matrix is a product of these kinds of matrices, and this can be done, for example, by going through the steps of the Gauss-Jordon method for solving systems of linear equations. And of course, you need to convince yourself that if a matrix $A$ multiplies volumes by a factor $a$, and another matrix multiplies volumes by a factor $b$, then $AB$ multiplies volumes by a factor $ab$. And you can see this by approximating any reasonable shape by lots of little parallelepipeds that fill in the shape.
Solution 2:
I am thinking aloud here, but maybe this will be helpful:
Gram schmidt on the vectors $ (v_1,\cdots, v_n) $ corresponds to gauss jordan on $ b_{i,j} $. For instance, let $ v_1 \rightarrow v_1 - c v_k $. Then: $$ b_{1,1} \rightarrow b_{1,1} - 2 c b_{1,k} + c b_{k,k} \\ b_{1,j} \rightarrow b_{1,j} - 2 c b_{k,j} \\ b_{j,1} \rightarrow b_{j,1} - 2 c b_{j,k} $$ One can arrive at this by row/column operations: $ b_{1,j}' = b_{1,j} - c b_{k,j} $ followed by $ b_{j,1}'' = b_{j,1}' - c b_{k,j}' $. So by gauss jordan you should be able to reduce to a diagonal matrix, corresponding to skewing your parallelepiped back into a hyper rectangle.