La fonction renvoie un vecteur de trois matrices imbriquées, P, L et U, de sorte que P . M = L . U. La matrice P représente la matrice de pivotement, et L et U représentent les matrices triangulaires inférieures et supérieures, respectivement.