Fonction llu_unit(A)
Calcule la matrice L (triangulaire inférieure) de la décomposition LU :
A = [[4,3],[6,3]]
L = llu_unit(A)
# L = [[1,0],[1.5,1]]
La diagonale de L est fixée à 1, ce qui simplifie le calcul de U.
Fonction ulu_unit(A)
Calcule la matrice U (triangulaire supérieure) correspondante :
U = ulu_unit(A)
# U = [[4,3],[0,-1.5]]
Ces deux fonctions permettent de reconstruire la matrice originale :
A = L * U