utilizando a fatoração LU sem pivoteamento parcial.
O conteúdo principal é decompor A em
A = LU,
onde L é triangular inferior com diagonal igual a 1 e U é triangular superior, e em seguida resolver por substituição direta (Ly=b) e regressiva (Ux=y).
Passo 2
Para obter L e U fazemos eliminações de Gauss sem trocar linhas. O primeiro multiplicador é