Solution (source code)

= Solution

Suppose $C=S\Lambda_CS^{-1}$ and $D=S\Lambda_DS^{-1}$ for an invertible $S$ and diagonal matrices $\Lambda_C,\Lambda_D$. Since diagonal matrices commute,
$$
CD=S\Lambda_C\Lambda_DS^{-1}=S\Lambda_D\Lambda_CS^{-1}=DC.
$$
Therefore \b[simultaneous diagonalization implies commutation]. This argument works over either $\mathbb R$ or $\mathbb C$ and explains why <simultaneous diagonalization> produces <commuting operators>.