Simultaneous triangularization of commuting matrices for the solution of polynomial equations
We present an extension of the QR method to simultaneously compute the joint eigenvalues of a finite family of commuting matrices. The problem is motivated by the task of finding solutions of a polynomial system. Several examples are included.