First time here? Check out the FAQ!
answered 2 years ago
Use echelon_form function with parameter transformation=True:
echelon_form
transformation=True
R, E = M.echelon_form(transformation=True) assert E*M == R