In Werner Greub's book Linear Algebra, 4th ed. on p. 230, he gives this proof of the normal form for a skew transformation on a finite-dimensional real inner product space. (Note Greub's convention for the matrix of a transformation is the transpose of that normally used with left-hand notation.)
I believe this proof is incorrect because it is not true in general that the $a_n$ defined form an orthonormal basis of the space. For example in $\mathbb{R}^4$, if we define the transformation $\psi$ by $$e_1\mapsto e_2\qquad e_2\mapsto -e_1\qquad e_3\mapsto e_4\qquad e_4\mapsto -e_3$$ where $e_i$ is the $i$-th standard basis vector, then $\psi$ is skew and $\varphi=\psi^2=-\iota$ is diagonalized by the standard basis. If we follow the proof for this example, we get $a_1=e_1$, $a_2=\psi e_1=e_2$, $a_3=e_2$, and $a_4=\psi e_2=-e_1$, so the $a_n$ do not form a basis of $\mathbb{R}^4$.
Does anyone see a way to salvage this proof while still retaining its spirit (in particular, avoiding use of complex numbers)?