I am looking for an elementary proof (if such exists) of the following: $$ AB=BA \quad\Longrightarrow\quad AB^T=B^TA, $$ where $A$ and $B$ are $n\times n$ real matrices, and $A$ is a normal matrix, i.e., $AA^T=A^TA$ - it is true for complex matrices as well, with $A^T$ replaced by $A^*$.
There is a non-elementary proof of this using the exponential of a matrix and properties of entire functions.
Update. In the first version of the question, both $A$ and $B$ were supposed to be normal, but as Shlomi correctly pointed out, this is true even in the case when only one of them is normal.
I'm not sure whether this qualifies as "elementary", but here's one way:
Note that $A$ and $B$ form a commuting family of matrices, so that $A$ and $B$ are simultaneously upper-triangularizable. That is, there is some unitary $U$ such that $U^*AU$ and $U^*BU$ are both upper triangular.
Note that $U^*AU$ and $U^*B U$ are both normal and upper-triangular, which means that they are diagonal. That is, we have $$ A = U\pmatrix{ \lambda_1 & &\\ &\ddots&\\ &&\lambda_n}U^* \quad B = U\pmatrix{ \mu_1 & &\\ &\ddots&\\ &&\mu_n}U^* $$ From here, it is easy to verify that $AB^* = B^*A$.
Alternatively, there's a nice inductive proof here if we use the following two facts:
With that, we may use the fact that $A$ and $B$ have a common eigenvector to write $$ A = U\pmatrix{ \lambda_1 & c^T\\ 0&\tilde A} U^* \quad B = U\pmatrix{ \mu_1 & d^T\\ 0&\tilde B}U^* $$ for a unitary $U$ (whose first column is a common eigenvector to $A$ and $B$). Because $A$ and $B$ are normal, we see that $c = d = 0$. Because $A$ and $B$ are commuting normal matrices, $\tilde A$ and $\tilde B$ are (smaller) commuting normal matrices.
By induction, we may conclude (as before) that for some unitary $U$, we have $$ A = U\pmatrix{ \lambda_1 & &\\ &\ddots&\\ &&\lambda_n}U^* \quad B = U\pmatrix{ \mu_1 & &\\ &\ddots&\\ &&\mu_n}U^* $$ and we may reach the desired conclusion, just like last time.