function [ U ] = rr( A )
U=A;
[rA cA]=size(A);
for j=1:rA-1
    for k=j+1:rA
        U(k,:)=U(k,:)-U(k,j)/U(j,j)*U(j,:);
    end
end
end

