For any real $n \times n$ skew-symmetric matrix $A,$ show there exists an orthogonal matrix $P$ such that $$PAP^T = \begin{pmatrix}{} \Lambda_1 & \\ & \Lambda_2 & \\ & & \ddots & \\ & & & \Lambda_\frac{n}{2} \end{pmatrix},$$ where $\Lambda_i = \begin{pmatrix}{} 0 & \lambda_i \\ -\lambda_i & 0 \\ \end{pmatrix}$, $\lambda_i \geq 0.$
I have found a decent proof of the above, but unfortunately it uses notions of linear (skew) transformations. I am having difficulty locating a proof that uses bilinear forms instead and I would like to try and understand that approach.