Let $\mathcal{M}_n(\mathbb{R})$ denote the vector space of matrices of size $n\times n$ on $\mathbb{R.}$.
For any $n\in\mathbb{N}$, one can find a basis of $\mathcal{M}_n(\mathbb{R})$ where all the matrices are non diagonalizable. However, can one find a basis of $\mathcal{M}_n(\mathbb{R})$ composed of diagonalizable matrices only ?
(Of course, those matrices don't have to be diagonalizable in the same basis)
For $n=2$, such a basis can be found; however for greater $n$s, I haven't found a general answer.
Let $E_{i,j}=((\delta_{i,j}))$ be an elementary matrix in $\mathcal{M}_n(\mathbb{R})$.
Some thoughts : one could easily gather $n(n+1)/$2 symmetric matrices to start a basis, namely the $E_{i,j}+E_{j,i}$ for $i\neq j$ and the $E_{i,j}$.
Had we been working in $\mathcal{M}_n(\mathbb{C})$, we could have added $n(n-1)/2$ antisymmetric $E_{i,j}-E_{j,i},i\neq j$ to get our basis, but that doesn't hold on $\mathcal{M}_n(\mathbb{R})$.