//reduce the matrix to row echelon form
clear;
A=[2,-1,3,2;1,4,0,-1;2,6,-1,5]
disp(A,"A:")
disp("Applying row transformations")
disp("R1=R1-2*R2");
A(1,:)=A(1,:)-2*A(2,:);