Factorisation PLU

Bonjour

Si $A$ est une matrice inversible, alors on peut trouver un triplet $(P,L, U)$ tel que $A=PLU$, avec $P$ une matrice de permutation, $L$ triangulaire inférieure unitaire et $U$ triangulaire supérieure inversible. Donc on a une factorisation $PLU$.
Comme par exemple : $$ \begin{pmatrix} 1& 2 \\ 3 & 4 \end{pmatrix}=\begin{pmatrix} 1& 0 \\ 0 & 1 \end{pmatrix}\begin{pmatrix} 1&0 \\ 3 & 1 \end{pmatrix}\begin{pmatrix} 1& 2 \\ 0 & -2 \end{pmatrix}.$$ Je veux écrire un algorithme en matlab de cette factorisation.

function [  C ] = my_matrix_mult( A,P,L,U )
[n,n]=size(A)=size(P)=size(L)=size(U);
G := Dom::Matrix(Dom::Real)(A);
[L, U, pivlist] := linalg::factorLU(G);

Mais la factorisation de cet algorithme nous donne deux matrices au lieu de trois, le cas ici.
Merci de vos aides.