Could you please let Diag() accept lists (systems) as arguments? The function Maxima(eigens_by_jacobi()) returns a list of eigenvalues and the spectral matrix. Currently, you have to apply mat2sys() to the eigenvector list if you want to evaluate the spectral representation:
It is perhaps just cosmetics but also one source of potential frustration removed.