Une factorisation A=LU exprime une matrice carrée comme produit d'une matrice triangulaire inférieure L et d'une matrice triangulaire supérieure U. Avec une diagonale unitaire pour L, la substitution avant est particulièrement simple. Une matrice inversible n'admet pas toujours cette forme sans échanges de lignes : un pivot principal nul peut nécessiter une permutation. Pour PA=LU, résoudre Ly=Pb, puis Ux=y. La permutation agit aussi bien sur le côté droit que sur la matrice coefficient. Des pivaux principaux non nuls justifient l'élimination usuelle sans échange, et non uniquement l'inversibilité.
Dans Ly=b, calculer y de haut en bas, en soustrayant les termes déjà connus ; diviser par la diagonale actuelle sauf si celle-ci vaut un. Dans Ux=y, calculer x de bas en haut. Pour chaque équation, substituer à nouveau dans la ligne originale comme vérification. Réutiliser LU pour plusieurs seconds membres évite de répéter l'élimination : une factorisation dense nécessite un travail d'ordre cubique selon la dimension, tandis que chaque résolution triangulaire requiert un travail d'ordre quadratique. L'arithmétique exacte peut vérifier de petits exemples ; le pivotage numérique réduit certains problèmes d'arrondi mais ne peut éliminer la mauvaise condition inhérente.
Pour B:V→W et A:W→Z, ker B est inclus dans ker(A∘B), mais ce dernier peut être plus grand. Les vecteurs supplémentaires sont ceux que B envoie dans ker A. Restreindre B à ker(A∘B) : son image est im B∩ker A et son noyau est ker B. Le théorème rang-noyau appliqué à cette restriction donne dim ker(A∘B)=dim ker B+dim(im B∩ker A). C'est l'intersection, et non tout le noyau de A, qui détermine la nullité supplémentaire.
Pour les endomorphismes de R⁶ ayant pour noyaux A=2 et B=3, le rang de B est 3 et la dimension de l'intersection peut varier de 0 à 2. Ainsi, la nullité de A∘B varie de 3 à 5. Les bornes dépendent de l'espace intermédiaire commun : pour des sous-espaces de dimensions r et s dans un espace de dimension m, leur intersection a une dimension comprise entre max(0,r+s−m) et min(r,s). Inverser la composition peut modifier sa nullité, même si AB et BA sont tous deux définis. Choisissez des espaces compatibles avant d'appliquer ces formules.