I have the following matrix: $$ A = \begin{bmatrix} 2 & 1 & 1 & 0 \\ 4 & 3 & 2 & 1 \\ 8 & 7 & 8 & 5 \\ 6 & 7 & 9 & 8 \\ \end{bmatrix} $$
and I want to apply Gauss elemination through a series of elementary matrix multiplications.
To cancel out the numbers under the first diagonal element of $A$, $$E_1 = \begin{bmatrix} 1 & 0 & 0 & 0 \\ -2 & 1 & 0 & 0 \\ -4 & 0 & 1 & 0 \\ -3 & 0 & 0 & 1 \\ \end{bmatrix} $$
Applying $E_1$, it works as expected
$$A^1 = E_1A = \begin{bmatrix} 2 & 1 & 1 & 0 \\ 0 & 1 & 0 & 1 \\ 0 & 3 & 4 & 5 \\ 0 & 4 & 6 & 8 \\ \end{bmatrix}$$
For the second column, however, starting with $E_2 = \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{bmatrix}$, how do I eliminate numbers under the second diagonal element of $A$?
If I want to apply on $E_2$ the row operations that would make the numbers under the second diagonal element of $A$ zero, then I would get for the third row of $E_2$, assuming it starts as an identity matrix, $I_n$ wouldn't I get:
$$E_2 = \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ -7 & 0 & 1 & 0 \\ -7 & 0 & 0 & 1 \\ \end{bmatrix}$$
The row operations in the context of A are apply $-7R_1 + R_3$ and $-7R_1 + R_4$ to cancel out the numbers under the second diagonal element of $A$ - so, the above shows I applied those to $I_n$
Doing the multiplication though, I get
$$A^2 = E_2E_1A = \begin{bmatrix} 2 & 1 & 1 & 0 \\ 0 & 1 & 0 & 1 \\ -14 & -4 & -3 & 5 \\ -14 & -3 & -2 & 8 \\ \end{bmatrix}$$
What am I doing wrong?
Edit: Calculation error in $A^1$ at position $(4,3)$. It should be $9 + (-3)\times 1 = 6 \neq 5$