I am trying to show that if $B_1$ and $B_2$ are orthonormal bases for $\mathbb{R}^n$, then the change of basis matrix $P$ from $B_1$ to $B_2$ is an orthogonal matrix.
I'm a bit stuck. I started with this: Let $x,y \in \mathbb{R}^n$. Then $[x]_{B_2}^t[y]_{B_2} = (P[x]_{B_1})^tP[y]_{B_1} = [x]_{B_1}^tP^tP[y]_{B_1}$. I want to show that $P^tP = I_n$. I know I need to apply the fact that $B_1$ and $B_2$ are orthonormal, but I don't see how to apply that to this expression. Am I approaching this properly or is there a better way to think about it?