qdiag function

Joint Approximate Diagonalization of a set of square, symmetric and real-valued matrices