What is the motivation to use the elimination method as opposed to the substitution method? I've always found the latter much easier. Are there any examples where the substitution method fails or when it is inferior to the elimination method?
For an example of the substitution method see:
Note that the elimination method is equivalent to the substitution method, where the choice of the variable for which we solve is constrained.
Example: Consider the $3 \times 3$ system $$ \begin{cases} x + y + z = 2\\ x + 2y + 3z = 5\\ x + 4y + 9z = 15 \end{cases} $$ With elimination, we would proceed as follows: $$ \pmatrix{1&1&1&2\\ 1&2&3&5\\ 1&4&9&15} \leadsto \pmatrix{1&1&1&2\\ 0&1&2&3\\ 0&3&8&13} \leadsto \pmatrix{1&1&1&2\\ 0&1&2&3\\ 0&0&2&4} \leadsto \pmatrix{1&1&1&2\\ 0&1&2&3\\ 0&0&1&2}. $$ At this point, we see that the original system of equations can be rewritten as $$ \begin{cases} x + y + z = 2\\ \quad \ \ y + 2z = 3\\ \quad \ \ \quad \ \ z = 2. \end{cases} $$ If we solve the first equation for $x$, substitute, solve the first of the resulting two equations for $y$, then the result is exactly the same. In particular, $$ x+y+z = 2 \implies x = 2-y-z. $$ Substituting this into the second and third equations yields $$ \begin{cases} (2-y-z) + 2y + 3z = 5\\ (2-y-z) + 4y + 9z = 15 \end{cases} \implies \begin{cases} y + 2z = 3\\ 3y + 8z = 13. \end{cases} $$ This is exactly what we have after the first step of elimination. Solving the second equation gives us $$ y + 2z = 3 \implies y = 3 - 2z,\\ 3(3-2z) + 8z = 13 \implies 2z = 4. $$ This is exactly what we have after the second step of elimination.