How to find row-echelon matrix form (not reduced) in Python?
Bill M's answer is correct. When you find the LU decomposition the U matrix is a  correct way of writing M in REF (note that REF is not unique so there are multiple possible ways to write it). To see why, remember that the LU decomposition finds P,L,U such that PLU = M. When L is full rank we can write this as U = (PL)-1M. So (PL)-1 defines the row operations that you have to preform on M to change it into U.
I am working on a project for my Linear Algebra class. I am stuck with a small problem. I could not find any method for finding the row-echelon matrix form (not reduced) in Python (not MATLAB).
Could someone help me out?
Thank you.
(I use python3.x)