No método Gauss Seidel realizamos uma decomposição A=M-N, onde M é uma matriz triangular inferior de A. O comando em Python no módulo import numpy ...
No método Gauss Seidel realizamos uma decomposição A=M-N, onde M é uma matriz triangular inferior de A. O comando em Python no módulo import numpy as np responsável por realizar esse procedimento é:
O método Gauss-Seidel não realiza uma decomposição A=M-N, mas sim uma decomposição A=L+D+U, onde L é uma matriz triangular inferior, D é uma matriz diagonal e U é uma matriz triangular superior.
Para realizar a decomposição A=L+D+U em Python utilizando o módulo numpy, podemos utilizar a função "tril" para obter a matriz triangular inferior, a função "triu" para obter a matriz triangular superior e a função "diag" para obter a matriz diagonal. O código ficaria assim:
```
import numpy as np
A = np.array([[4, 1, 1], [2, 7, 1], [1, 2, 6]])
L = np.tril(A, k=-1)
D = np.diag(np.diag(A))
U = np.triu(A, k=1)
print("Matriz A:")
print(A)
print("Matriz L:")
print(L)
print("Matriz D:")
print(D)
print("Matriz U:")
print(U)
```
Lembrando que esse código realiza a decomposição A=L+D+U e não A=M-N.
0
2
Faça como milhares de estudantes: teste grátis o Passei Direto
Compartilhar