n=5;
A=rand(n,n);
b=rand(n,1);
for i=1:n 
    for j=i+1:n
        A(j,i)=0;
    end
end
A
b
x=rwe(A,b)
b-A*x