I’m trying to implement a householder tridiagonalization routine to tridiagonalize small (n < 10) hermitian matrices in matlab. This is what I have so far: