e queremos decompor em uij = aij - sumk=1^i-1 lik,ukj, lij = 1ujjl(aij - sumk=1^j-1 lik,ukjr) , onde é triangular inferior com 1’s na diagonal e é triangular superior. Isso faz parte de “fatoração LU” ou “eliminação de Gauss sem pivô” e a gente vai usar o método de Doolittle (que já coloca 1’s na diagonal de ).
O que a gente vai precisar:
Montar com 1’s na diagonal e variáveis nos outros lugares.
Montar com variáveis na parte triangular superior.