function A = qrstep0(A) % one step of the QR iteration % with NO shifting [m,n]=size(A); rho = 0; [Q,R]=qr(A-rho*eye(n,n)); A = R*Q + rho*eye(n,n);