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)
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.
How to find row-echelon matrix form (not reduced) in Python?