Why does Gaussian elimination work?

It works by preserving solution set through row operations.

Gaussian elimination uses elementary row operations that do not change the solution space of Ax=b. It converts system into upper triangular form for easy back substitution.