N=10;
A=eye(N);
q=rand(N,1);
A*q
