Solution
= Solution
Using the unit length of $q_i$,
$$
P_i^2=q_iq_i^Tq_iq_i^T=q_i(q_i^Tq_i)q_i^T=q_iq_i^T=P_i.
$$
Then ordinary <matrix> expansion gives $(I_n-P_i)^2=I_n-2P_i+P_i^2=I_n-P_i$. \b[Both $P_i$ and its complementary <matrix> are <projectors>.] The second calculation in fact holds for every <idempotent matrix>, without requiring the <matrix> to equal its <transpose>.