I would like to explain to my students that the determinant of a matrix $A$ is in fact the measure of how the volume changes under said matrix. To do that, I would like to start from the $n$-dim. cube (represented by the standard basis $e_1,\ldots,e_n$ of $\mathbb{R}^n$) and derive the formula for the volume of the parallelepiped determined by $Ae_1,\ldots,Ae_n$. Of course, I could go to the exterior algebra $\Lambda^n(\mathbb{R}^n)$ and work there, but I don't think this would be very pedagogical...
Thus, my question is: is there a down-to-earth way to recover the determinant of $A$ as the volume of the parallelepiped defined above? And if not, how can one make a compelling argument of why taking the exterior product $Ae_1\wedge\ldots\wedge Ae_n$ gives the volume?
My preferred pedagogical approach is to start with examples in $2$ and $3$ dimensions. In $2D$ we can easily see that the columns of the matrix $A$ are the vectors $\vec u= A\vec i$ and $\vec v=A \vec j$ ( where $\{\vec i, \vec j \}$ is the standard basis) and it is easy to show that the ''signed'' area of the parallelogram with sides $\vec u$ and $\vec v$ is the cross product $ \vec u \times \vec v = \det A$.
In a similar way we can show that the mixed product of the vectors $\vec u=A\vec i$,$\vec v=A \vec j$ and $\vec w=A\vec k$ is the volume of the of the transformed parallelepiped,and these vectors are the columns of the matrix $A$ so that $\vec u\cdot(\vec v \times \vec w)=\det A$.
Now the generalization to $n-$dimensional space can be done noting the properties of the determinant, especially the multi-linearity, and noting that we can go from the volume of a $n-$cube to an $n-$polyhedron by linear tranformations of the sides.