function D=qriter(A,nmax)
n=max(size(A)); D=A;
for k=1:nmax
[Q,R]=qr(D);
D=Q'*D*Q;
end; 
