a) É possível verificar a existência e unicidade da solução? Se sim, mostre. Se não, justifique.
b) Resolva-o pelo método da Eliminação de Gauss (sem pivotação).
c) A solução é exata? Mostre através do vetor de resíduos.
d) Resolva-o pelo método da Eliminação de Gauss com pivotação parcial.
e) Resolva-o pelo método da decomposição LU.
f) Se formos resolver o sistema por métodos iterativos (de Jacobi ou Gauss-Seidel), temos garantia da convergência do método? Justifique/explique.
pede (a) existência e unicidade da solução, (b) resolução por Eliminação de Gauss sem pivotação, (c) verificação exata via vetor de resíduos, (d) resolução por Eliminação de Gauss com pivotação parcial, (e) decomposição LU, (f) condição de convergência de métodos iterativos. Vamos usar teoria de determinante, eliminação de Gauss, decomposição LU e critérios de convergência (dominância diagonal, SPD).
Para (c) calculamos o vetor de resíduos r=b-Ax. Com x=(39/19,42/19,-47/19), verifica-se
Ax=(-5,-7,4)^T=b,
logo r=(0,0,0)^T. A solução é exata.
Passo 5
Para (d) com pivotação parcial trocamos a linha 1 pela linha 2 (maior pivô em módulo na coluna 1), e repetimos a eliminação. O resultado final da solução é o mesmo:
Resolvendo Ly=b e depois Ux=y dá novamente \bigl(x_{1},x_{2},x_{3}\bigr)=(39/19,42/19,-47/19).
Passo 7
Para (f) a condição suficiente de convergência de Jacobi ou Gauss–Seidel é dominância diagonal estrita ou SPD em A. Aqui as linhas 1 e 3 não satisfazem \lvert a_{ii}\rvert\ge\sum_{j\ne i}\lvert a_{ij}\rvert,
e A não é simétrica positiva definida. Portanto não temos garantia geral de convergência.
Resposta
a) Sim. \det(A)=19\neq0, logo existe solução única.
b) x_{1}=39/19,\;x_{2}=42/19,\;x_{3}=-47/19.
c) O vetor de resíduos é r=b-Ax=(0,0,0)^T, solução exata.
d) Com pivotação parcial chega-se ao mesmo (39/19,42/19,-47/19).
e) L=\begin{pmatrix}1&0&0\\\ 2&1&0\\tfrac12&-\tfrac74&1\end{pmatrix}, U=\begin{pmatrix}2&-3&1\\\ 0&-2&-3\\\ 0&0&-\tfrac{19}{4}\end{pmatrix},
solução (39/19,42/19,-47/19).
f) Não há garantia de convergência de Jacobi ou Gauss–Seidel pois A não é estritamente diagonal dominante nem SPD.