usando o método de Gauss–Seidel com chute inicial \mathbf{x}^{(0)}=(0,0,0,0)^{\mathrm{T}}, arredondando cada iteração para duas casas decimais e parando quando a variação máxima entre iterações for menor que 0,01.
Passo 2
Para aplicar Gauss–Seidel, isolamos cada variável em cada equação, usando sempre os valores mais recentes: