Adding a multi-shift Hessenberg approach for computing pseudospectra.