A = randn(m, n); % for columns A = A ./ repmat(sqrt(sum(A.^2, 1)),m,1); % for rows A = A ./ repmat(sqrt(sum(A.^2, 2)),1,n);