Gaussian elimination
The above procedure can be repeatedly applied to solve the equation multiple times for different b. In this case it is faster (and more convenient) to do an LU decomposition of the matrix "A" once and then solve the triangular matrices for the different b, rather than using Gaussian elimination each time. The matrices "L" and "U" could be thought to have "encoded" the Gaussian elimination process.