Adding sequential implementation of a multi-shift Hessenberg solve based upon that of Alg. 2 in Greg Henry's 'The shifted Hessenberg system solve computation'.