Let $\theta$ and $\psi$ be symmetric bilinear forms on a finite-dimensional real vector space $V$, and assume $\theta$ is positive definite. Show that there exists a basis $\{v_1,\ldots,v_n\}$ for $V$ and $\lambda_1,\ldots,\lambda_n\in\mathbb{R}$ such that $$\theta(v_i,v_j)=\delta_{i,j}\quad\text{and}\quad\psi(v_i,v_j)=\delta_{ij}\lambda_i$$ where $\delta_{ij}$ is the Kronecker delta function.
I think it's enough to choose a basis $\{w_1,\ldots,w_n\}$ for which the matrix representations of $\theta$ and $\psi$ are both diagonal. Then $\left\{\frac{w_1}{\sqrt{\theta(w_1,w_1)}},\ldots,\frac{w_n}{\sqrt{\theta(w_n,w_n)}}\right\}$ is the required basis.
The catch is finding a basis which simultaneously diagonalizes $\theta$ and $\psi$.
Take $V = \mathbb{R}^n$ for simplicity, and let $\theta (x,y) = \langle x, A y \rangle$, and $\psi (x,y) = \langle x, B y \rangle$, with $A,B \in \mathbb{R}^{n\times n}$. By assumption $A=A^T >0$, and $B= B^T$.
Since $A>0$, it has a Cholesky decomposition $A = Q^T Q$ (with Q invertible, of course). Now consider the matrix $Q^{-T} B Q^{-1}$. Since this is symmetric, it can be diagonalized by an orthogonal $U$ to give $U^T Q^{-T} B Q^{-1} U = \Lambda$, where $\Lambda$ is diagonal. Also notice that $U^T Q^{-T} A Q^{-1} U = I$ (since $Q^{-T} A Q^{-1} = I$).
Now choose a basis $v_i = Q^{-1} U e_i$, where $e_i$ is the standard basis for $\mathbb{R}^n$. Then we have $$ \theta(v_i, v_j) = \langle Q^{-1} U e_i , A Q^{-1} U e_j \rangle = \langle e_i , U^T Q^{-T} A Q^{-1} U e_j \rangle = \langle e_i , e_j \rangle, $$ and similarly $$ \psi(v_i, v_j) = \langle Q^{-1} U e_i , B Q^{-1} U e_j \rangle = \langle e_i , U^T Q^{-T} B Q^{-1} U e_j \rangle = \langle e_i , \Lambda e_j \rangle, $$ from which the desired result follows.