I am trying to evaluate a numerical integral and running into slow convergence issues. My problem is as follows: I have a 6×6 Hermitian matrix, H[x,y], whose elements are functions of two parameters x and y. I find the eigenvalues of this matrix, which are real eigs[x_,y_]:=Eigenvalues[H[x,y]] and expressed using the Root function as rootsRead more