Trying to prove that for any two matrices $A,B$ representing the same symmetric bilinear form $f:V \times V \to \mathbb{R}$, there is an invertible matrix $P$, such that $B = P^{T}AP$ ?
I have a method for constructing such $P$ by repeated change of bases (when it is a quadratic form) but can't see a way to prove this generally?
This is just a change of basis. Let $A$ and $B$ be the matrices of the same bilinear form, but with respect to two ordered bases $\{x_1,\ldots,x_n\}$ and $\{y_1,\ldots,y_n\}$ respectively. Let $P$ be the matrix such that $y_j=\sum_ip_{ij}x_i$, i.e. the matrix such that $\sum_jc_jy_j=\sum_i(Pc)_ix_i$ for any $c\in\mathbb R^n$. Then $B=P^TAP$.