A eliminação gaussiana para solução de sistemas de equações lineares está entre os métodos mais eficientes. Escrevemos um script no Scinotes do Scilab (www.scilab.org) para efetuar o passo da eliminação progressiva.
Antes de eliminar os termos de uma coluna, a matriz sofre a operação de pivotamento, que consiste em colocar como pivot o elemento de maior valor absoluto entre os termos da coluna a sofrer a eliminação, excetuando os termos acima do termo da diagonal.
Segue-se, daí, a eliminação progressiva.
Depois de proceder a eliminação em cada coluna, a matriz passa por novo pivotamento.
A linha 51 é apenas para mostrar a matriz final, isto é, depois de pivotar e eliminar.
Antes de eliminar os termos de uma coluna, a matriz sofre a operação de pivotamento, que consiste em colocar como pivot o elemento de maior valor absoluto entre os termos da coluna a sofrer a eliminação, excetuando os termos acima do termo da diagonal.
No próprio script preparamos as matrizes de entrada para as operações a serem aplicadas.
Segue-se, daí, a eliminação progressiva.
Depois de proceder a eliminação em cada coluna, a matriz passa por novo pivotamento.
A linha 51 é apenas para mostrar a matriz final, isto é, depois de pivotar e eliminar.
Uma referência importante em métodos numéricos são os livros de Chapra-Canale.
show
ResponderExcluir