Given a set of n<d commuting operators, either degenerate or non-degenerate, in a d-dimensional Hilbert space, is there an effective analytical method of finding an orthonormal basis composed of d eigenvectors common to all the operators in the set?
The operators are dxd complex square matrices...