Introduction
Guessing and substituting works for two equations, but not for ten. Gaussian elimination is a systematic recipe that works for any linear system: simplify the equations step by step, without changing their solutions, until the answer can be read off. It is the algorithm behind np.linalg.solve.