function A = rdiag(n)
A = diag(fix(10*rand(1,n)));