Hi, I'm a new user.
I have a question about the subroutine dhsein.
I have a 10x10 real matrix in Hessenberg form. I need to compute all the eigenvectors, so select is a true array of 10 elements. So I choose mm =10.
Sometimes the routine crashes because m exeed 10 (it is 11). I don't understand how this is possible.