
The most difficult part is that how can we use the row operations to converse the "augmented matrix" into the "solution matrix" that contains "identity matrix" and our "solution matrix").
For example we have the system of equation is : 2x + y = 5 and 5x + 3y = 13. How can we change to solution matrix ?