Write Gauss Elimination Code

I need to write a function code in matlab that can reduce any 3x3 matrix into reduced or echelon form. This code should perform the steps of “rref” without using it. The answer should show all the steps used to get the matrix into reduced or echelon form.