I am trying to write my own rref and have the following code. It appears to have worked for the first column but not for all other columns for a test augmented matrix I entered. Can you please assist? Thanks!

%Loop through columns, then rows

%j=1:n-1 because matrix is augmented and I only want to check the square matrix values; however, row operations should be performed on the entire row of the matrix

for j=1:n-1

    for i=1:m

%if i=j, the code already performed an operation to transform the value to 1, so I want to skip matrix(i,j) if i=j

if i~=j

        if matrix(i,j)~=0
            if matrix(i,j)>0
            matrix(i,:)=matrix(i,:)-matrix(i,j)*matrix(j,:);
            elseif matrix(i,j)<0
                matrix(i,:)=matrix(i,:)+(abs(matrix(i,j)))*matrix(j,:);
            end
        end
    end
    end
end

Answer 1

[Sanjana Ramakrishnan]

Refer MATLAB source code for reference by typing the below command in MATLAB command window:

>>edit rref

You might also find the below file exchange link useful:

https://www.mathworks.com/matlabcentral/fileexchange/21583-fast-reduced-row-echelon-form