factorisation LU
factorisation LU