How can I make a program that solves for the reduced echelon form of a matrix?

I dont know how to show it but it involves a lot of row operations between each row of the matrix

confused :|

here's the definition of a echelon form

Row echelon form - Wikipedia, the free encyclopedia