Version2.1:

- update in MATRIX.f90: the matrix_diagonalize routines had default setup with
overwrite_matrix=FALSE. In this way overwriting with the matrix of the eigenvectors
that was often assumed to be true was never done in reality, leading to possible problems
in many circumstances. default is now 'V' that is TRUE.
TODO add a logical variable as in scipy.