function A = RandomStoch(n)
% A is a random nxn matrix with positive entries and column sums equal
% to one..
A = rand(n,n);
for j=1:n
A(:,j) = A(:,j)/sum(A(:,j));
end