O objetivo é aplicar iterações de Gauss–Seidel até que a diferença entre duas iterações consecutivas seja menor que \varepsilon (iremos usar a norma infinito: \|X^{(k+1)}-X^{(k)}\|_{\infty}<\varepsilon).
Passo 2
Vamos montar as fórmulas de iteração do Gauss–Seidel isolando cada variável em sua equação. A ideia do Gauss–Seidel é atualizar as componentes em sequência e, sempre que possível, usar os valores já atualizados na mesma iteração (isso acelera a convergência).