I'm writing a code that mostly consists of matrices manipulation and, in particular, diagonalizati