Firstly, I will note that I have looked at many sources on Google, and also on Stack exchange, but have not been able to find a hint of how to go about proving this the way I am doing- or at least not one which I could connect with the symbols/notation I am using here. I think that the first answer in this post is quite similar, but I am not sure. I am not familiar with manipulating the symbols/notation that was used there.
So far I have:
$$(AB)_{ij}=p_{ij}=\sum _k (A)_{ik}(B)_{kj}$$
Then
$$\det(AB)=\sum_{j_1j_2...j_n}\epsilon _{j_1j_2...j_n}p_{1j_1}p_{2j_2}...p_{nj_n}$$
Where $\epsilon _{j_1j_2...j_n}$ is the Levi-Civita or epsilon tensor
$$\det(AB)=\sum_{j_1j_2...j_n}\epsilon _{j_1j_2...j_n}\sum_{k_1} (A)_{1k_1}(B)_{k_1j_1}...\sum_{k_n} (A)_{nk_n}(B)_{k_nj_n}$$
$$=\sum_{j_i,k_i}\epsilon _{j_1j_2...j_n}(A)_{1k_1}...(A)_{nk_n}(B)_{k_1j_1}...(B)_{k_nj_n}$$
Now I am stuck. But what is frustrating me is that I see exactly what needs to be done for me o get the result. What I am looking for is
$$\sum_{j_i,k_i}\epsilon _{j_1j_2...j_n}\epsilon _{k_1k_2...k_n}(A)_{1k_1}...(A)_{nk_n}(B)_{1j_1}...(B)_{nj_n}$$
So I somehow need to get another epsilon oiver the $k$ terms out of $(B)_{k_1j_1}...(B)_{k_nj_n}$, but I don't know how to motivate this?
Therefore, $\det(AB)$ is the volume of the set that results from acting on the unit cube, first, by $B$, and then by $A$. Multiplicativity follows from this alone.
For a rigorous proof without coordinates, see the definition of determinant in Halmos's "Finite-dimensional vector spaces."
– avs May 08 '17 at 16:51